Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,472,437 T→C 8.0% N33D (AAT→GAT)  pyrL ← pyrBI operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,4370TC 8.0% 583.4 / 28.7 187N33D (AAT→GAT) pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (86/86);  new base C (8/7);  total (94/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 189 total aligned reads displayed.

CGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGG  >  NC_000913/4472291‑4472581
                                                                                                                                                   |                                                                                                                                                
  ggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAgg                                                                                                                                              <  4:955422/149‑1 (MQ=255)
     ttGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCAATCGAGGGGc                                                                                                                                           <  1:1057787/149‑1 (MQ=255)
        cTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTgg                                                                                                                                        <  3:495929/149‑1 (MQ=255)
        cTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTgg                                                                                                                                        <  4:497328/149‑1 (MQ=255)
        cTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTAATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCAATCGAGGGGCttc                                                                                                                                        <  3:582073/149‑3 (MQ=255)
         tttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                                                                                       >  2:1246252/1‑149 (MQ=255)
         tttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                                                                                       >  3:1277918/1‑149 (MQ=255)
         tttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                                                                                       >  4:1213786/1‑149 (MQ=255)
         tttCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCAATCGAGGGGc                                                                                                                                           >  4:570755/1‑145 (MQ=255)
         tttCAGTTTAGCCGCTGTCGCCAGCACCAGACTAGGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                                                                                       >  2:1121627/1‑149 (MQ=255)
            cAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAAt                                                                                                                                    <  1:407930/149‑1 (MQ=255)
              gTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGGTATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATgg                                                                                                                                  >  2:406753/1‑149 (MQ=255)
              gTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATgg                                                                                                                                  <  2:1283010/149‑1 (MQ=255)
              gTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATgg                                                                                                                                  >  4:419955/1‑149 (MQ=255)
               tttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATggg                                                                                                                                 <  2:93938/149‑1 (MQ=255)
                ttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGt                                                                                                                                >  1:224629/1‑149 (MQ=255)
                        tGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACg                                                                                                                        >  1:916541/1‑149 (MQ=255)
                           tgcaaGCACCAGATTAAGGACATCACAACTAAGTTCGTTTATGCTAATGAAATGTGTCTAAAATAACCGATTAGCTATTTTTTATCCTCTGTCGCATGGGA‑AAAAAAAAGCCCCTCGATTTAGGGGCTGGGAATGGGTGATGAACGGga                                                                                                                     <  2:800043/145‑1 (MQ=255)
                           cGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGga                                                                                                                     <  2:1075194/149‑1 (MQ=255)
                            gCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaa                                                                                                                    <  1:164170/149‑1 (MQ=255)
                             ccAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTCTCTGATATAACGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaag                                                                                                                   <  3:1295724/149‑1 (MQ=255)
                                  aCCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaaaa                                                                                                              >  1:654692/1‑149 (MQ=255)
                                   ccAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaa                                                                                                                    <  2:1212536/142‑1 (MQ=255)
                                      gATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGAt                                                                                                                             <  3:892284/130‑1 (MQ=255)
                                      gATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGAt                                                                                                                             >  4:892284/1‑130 (MQ=255)
                                       aTTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGt                                                                                                                                <  2:1250698/126‑1 (MQ=255)
                                       aTTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGt                                                                                                                                >  1:1250698/1‑126 (MQ=255)
                                       aTTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCa                                                                                                         <  3:992361/149‑1 (MQ=255)
                                        ttAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGa                                                                                                                              >  3:841279/1‑127 (MQ=255)
                                        ttAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGa                                                                                                                              <  4:841279/127‑1 (MQ=255)
                                            ggTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCa                                                                                                    >  3:409528/1‑149 (MQ=255)
                                              tCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGc                                                                                                  >  1:719735/1‑149 (MQ=255)
                                                 tCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGc                                                                                                          <  3:634350/138‑1 (MQ=255)
                                                    cgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCttt                                                                                            <  4:934673/149‑1 (MQ=255)
                                                     gACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtttt                                                                                           <  3:427818/149‑1 (MQ=255)
                                                      aCTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCttttt                                                                                          <  1:1068452/149‑1 (MQ=255)
                                                       cTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGAGCAACGGGAAGAAAAACGGCAGGCCAGCGTCtttttt                                                                                         <  3:330128/149‑1 (MQ=255)
                                                            gTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCt                                                                                                                                          >  3:33402/1‑95 (MQ=255)
                                                            gTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCt                                                                                                                                          <  4:33402/95‑1 (MQ=255)
                                                                tttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCgg                                                                                <  4:601351/149‑1 (MQ=255)
                                                                     ggAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAc                                                                                    <  3:553240/140‑1 (MQ=255)
                                                                     ggAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAc                                                                                    >  4:553240/1‑140 (MQ=255)
                                                                      gAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGt                                                                                                                                >  1:518507/1‑95 (MQ=255)
                                                                      gAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtttttt                                                                                         >  4:213438/1‑134 (MQ=255)
                                                                      gAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtttttt                                                                                         <  3:213438/134‑1 (MQ=255)
                                                                      gAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGAc                                                                          >  2:699301/1‑149 (MQ=255)
                                                                       aaaTGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATg                                                                                                                                   <  2:1258510/91‑1 (MQ=255)
                                                                          tGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCa                                                                                                    <  3:166691/119‑1 (MQ=255)
                                                                            atatGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATg                                                                    >  2:202373/1‑149 (MQ=255)
                                                                              atGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                                                                  <  4:157534/149‑1 (MQ=255)
                                                                              atGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                                                                  <  1:1272471/149‑1 (MQ=255)
                                                                                 tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGaa                                                               >  3:870602/1‑149 (MQ=255)
                                                                                 tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGaa                                                               >  4:1059769/1‑149 (MQ=255)
                                                                                  ttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac                                                              >  1:1477466/1‑149 (MQ=255)
                                                                                     tGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacaca                                                           >  4:1261287/1‑149 (MQ=255)
                                                                                      gATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTGTTTTTTCAGACGCGGTAAGACAAAATGTCGAacacac                                                          <  1:1214048/149‑1 (MQ=255)
                                                                                      gATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacacac                                                          <  4:1277918/149‑1 (MQ=255)
                                                                                        tataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTg                                                        <  3:1408844/149‑1 (MQ=255)
                                                                                           aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAAc                                                     <  2:1191317/149‑1 (MQ=255)
                                                                                            gCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtttt                                                                                           >  2:256769/1‑110 (MQ=255)
                                                                                            gCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAAAGTCGAACACACTGAAcc                                                    >  4:1185082/1‑149 (MQ=255)
                                                                                            gCGGATTAGCCATCTGTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtttt                                                                                           <  1:256769/110‑1 (MQ=255)
                                                                                               gATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATa                                                 <  1:347399/149‑1 (MQ=255)
                                                                                               gATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATa                                                 <  3:1090526/149‑1 (MQ=255)
                                                                                               gATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATa                                                 <  2:528229/149‑1 (MQ=255)
                                                                                                   aGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtc                                             >  1:629146/1‑149 (MQ=255)
                                                                                                    gCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcc                                            >  4:1159527/1‑149 (MQ=255)
                                                                                                      cATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACAtcctc                                          <  3:1416264/149‑1 (MQ=255)
                                                                                                          ttttATCTCCTGACGCCTGGGCAAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCg                                       >  4:399139/1‑149 (MQ=255)
                                                                                                          ttttATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGt                                                                   >  3:1172144/1‑120 (MQ=255)
                                                                                                          ttttATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCgg                                      >  2:371197/1‑149 (MQ=255)
                                                                                                              aTCTCCTGCCGCCTGGGC‑AAAAAAAAGCCCCTCGATTGGGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCaaa                                  <  1:688217/149‑1 (MQ=255)
                                                                                                                   cTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGa                                                                           >  4:1412451/1‑103 (MQ=255)
                                                                                                                   cTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                                                                  >  2:850457/1‑112 (MQ=255)
                                                                                                                     gACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCg                           >  4:1348975/1‑149 (MQ=255)
                                                                                                                          cTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac                                                              >  1:475554/1‑109 (MQ=255)
                                                                                                                          cTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac                                                              <  2:475554/109‑1 (MQ=255)
                                                                                                                              gC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATAc                  >  2:1268262/1‑149 (MQ=255)
                                                                                                                                aaaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACt                 >  1:1189409/2‑149 (MQ=255)
                                                                                                                                 aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACtt                >  1:1324660/1‑148 (MQ=255)
                                                                                                                                 aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc                >  1:988770/1‑149 (MQ=255)
                                                                                                                                 aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc                >  1:1287237/1‑149 (MQ=255)
                                                                                                                                 aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc                >  1:789872/1‑149 (MQ=255)
                                                                                                                                 aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACAGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc                >  4:1104546/1‑149 (MQ=255)
                                                                                                                                   aaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCAt              >  1:560605/1‑149 (MQ=255)
                                                                                                                                      ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGaa                                                               >  2:646955/3‑97 (MQ=255)
                                                                                                                                      ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATa                                                 >  2:346944/3‑111 (MQ=255)
                                                                                                                                      ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCTGCATTATACTCATCGt           >  3:931632/3‑149 (MQ=255)
                                                                                                                                      ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGt           >  1:926086/3‑149 (MQ=255)
                                                                                                                                      ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGt           >  4:441656/3‑149 (MQ=255)
                                                                                                                                       caGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTc          >  1:28308/2‑149 (MQ=255)
                                                                                                                                       aaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTc          >  1:749003/1‑149 (MQ=255)
                                                                                                                                           cccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGaa      <  1:1082168/149‑1 (MQ=255)
                                                                                                                                           cccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGaa      >  3:782608/1‑149 (MQ=255)
                                                                                                                                           cccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGaa      >  2:668194/1‑149 (MQ=255)
                                                                                                                                             cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc    >  1:881131/1‑149 (MQ=255)
                                                                                                                                             cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc    >  1:942664/1‑149 (MQ=255)
                                                                                                                                             cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc    >  2:1161338/1‑149 (MQ=255)
                                                                                                                                               cGATTGAGGGGCTGGGCATGGGTGATCAAAGGGAAGAAAAACGGCAGGGCAGCGTCTCTTTTCAGACTCGGTAAGACAAAGCGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgcgg  <  2:1049617/149‑1 (MQ=255)
                                                                                                                                                  ttGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCa                                                   <  2:993594/97‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGC‑AAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGG  >  NC_000913/4472291‑4472581

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑