Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,066,310 A→T 5.8% H288Q (CAT→CAA ygfI ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,066,3100AT 5.8% 254.1 / 4.5 86H288Q (CAT→CAAygfIputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (42/39);  new base T (3/2);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

CGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAA  >  NC_000913/3066164‑3066456
                                                                                                                                                   |                                                                                                                                                  
cGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGc                                                                                                                                                 <  2:1006571/150‑1 (MQ=255)
       aTCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCa                                                                                                                                          >  2:1040974/1‑150 (MQ=255)
          tGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCg                                                                                                                                       >  1:996457/1‑150 (MQ=255)
             aTCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTgg                                                                                                                                    >  1:511294/1‑150 (MQ=255)
               cGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACg                                                                                                                                  <  1:377106/150‑1 (MQ=38)
                gCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGACACGTCTTCGACTTAGACCGTAATTTAATTGCAACCGGCTCGATGCTGGGATAATGACTCGTCCCGCCACGATGCTTCAGCATATTGCTGAAGATCGTGACgg                                                                                                                                 <  1:229555/150‑1 (MQ=37)
                   aTCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGtat                                                                                                                              >  1:701511/1‑150 (MQ=255)
                    tCTGGAAAATGCGCCTGGCGGACATATTCACCACGAGCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGtata                                                                                                                             <  1:219959/150‑1 (MQ=255)
                      tGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACCTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAAc                                                                                                                           <  1:614599/150‑1 (MQ=255)
                       ggAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACg                                                                                                                          <  2:23049/150‑1 (MQ=255)
                       ggAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACg                                                                                                                          >  2:1199257/1‑150 (MQ=255)
                       ggAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCAAGATCATCAGCAATATGCTGAAGCATCGAGGCGGTATAACg                                                                                                                          >  1:166741/1‑150 (MQ=255)
                          aaaTGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACATACACCGTAATTAAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgatg                                                                                                                       <  1:1122598/150‑1 (MQ=255)
                              gcgcCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTTCACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGGGGTATAACGATGAttt                                                                                                                   >  2:187162/1‑150 (MQ=255)
                               cgcCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGACGCATCGTGGCGGTATAACGATGATTTc                                                                                                                  >  1:850957/1‑150 (MQ=255)
                                gcCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCt                                                                                                                 <  1:339620/150‑1 (MQ=255)
                                   tGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGc                                                                                                                                   <  1:1089987/129‑1 (MQ=255)
                                   tGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGc                                                                                                                                   >  2:1089987/1‑129 (MQ=255)
                                   tGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAt                                                                                                              <  1:521133/150‑1 (MQ=255)
                                    ggCGGACATATTCACGACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAtt                                                                                                             <  1:54241/150‑1 (MQ=255)
                                     gCGGACATATTCCCCACGAGCTTAGCCACGGCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGAGTCGGCGCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACTATGATTTCTGATTg                                                                                                            <  2:274165/150‑1 (MQ=255)
                                     gCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                                                            <  1:798585/150‑1 (MQ=255)
                                     gCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                                                            <  1:267913/150‑1 (MQ=255)
                                     gCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                                                                                                            <  2:745935/150‑1 (MQ=255)
                                       ggACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGaa                                                                                                          <  1:322065/150‑1 (MQ=255)
                                       ggACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGaa                                                                                                          >  2:813702/1‑150 (MQ=255)
                                       ggACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGaa                                                                                                          <  1:562393/150‑1 (MQ=255)
                                        gACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAAt                                                                                                         <  1:426417/150‑1 (MQ=255)
                                         aCATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATa                                                                                                        <  1:114065/150‑1 (MQ=255)
                                          cATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaa                                                                                                       >  2:858612/1‑150 (MQ=255)
                                          cATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaa                                                                                                       >  1:561218/1‑150 (MQ=255)
                                           atatTCACCACTATCTTACCCACGGCTTAGACTTAAACCGTAATTTAATTGCAACCGGCTCGCTGCTGGGGTAATTACTCGTCCCGTCACGATCTTCAGCAAAATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa                                                                                                      <  2:91049/150‑1 (MQ=255)
                                                caccacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGatat                                                                                                 <  2:1022783/150‑1 (MQ=255)
                                                caccacGATCTTAGCCACGTCTTCGACTGACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGatat                                                                                                 >  1:459392/1‑150 (MQ=255)
                                                  ccacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCCACCGGGTCGATGCTGGGTTAATGAGTCGTCACGTCACGCTCTTCAGCAATATGCTGAACCATCGTGGAGGTATAAGGAAGATTTATTATTTATTAAAGATAtta                                                                                               >  2:983701/1‑149 (MQ=255)
                                                   cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTAAAACGATGATTTCTGATTGAATAAAGATAttac                                                                                              >  1:44120/1‑148 (MQ=255)
                                                    acGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCt                                                                                             >  1:1205002/1‑150 (MQ=255)
                                                     cGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTg                                                                                            <  1:1300251/150‑1 (MQ=255)
                                                       aTCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATCATTTCTGATTGAATAAAGATATTCCTGCa                                                                                          >  1:749925/1‑150 (MQ=255)
                                                        tCTTAGCCACGTCTTCGACTTACACCGTAAGTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATACCGATGATTTCTGATTGAATAAAGCACTTCCTGcct                                                                                         >  1:708469/1‑148 (MQ=255)
                                                           tAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCAtata                                                                                      <  2:1290345/150‑1 (MQ=255)
                                                             gCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGTTATAACGATGATTTCTGATAGAATAAAGATATTCCTGCATATACa                                                                                    >  1:904704/1‑150 (MQ=255)
                                                               cACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATATCGATGATTACTGATTGAATAAAGATATTCCTGCATATACATc                                                                                  >  1:633718/1‑150 (MQ=255)
                                                               cACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGCATTAAGATATTCCTGCATATACCTc                                                                                  >  1:1137352/1‑150 (MQ=255)
                                                                   tCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                                                              <  1:1009343/150‑1 (MQ=255)
                                                                   gcTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGcc                                                                              <  2:994516/149‑1 (MQ=255)
                                                                       cGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAATATATTCCTGCATATACATCGGCCtgtg                                                                          >  1:124828/1‑150 (MQ=255)
                                                                       cGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCtgtg                                                                          <  1:880985/150‑1 (MQ=255)
                                                                       cGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCAtgtg                                                                          >  1:1179981/1‑150 (MQ=255)
                                                                         aCTTACACCGTAAATTAATTGCACCCGGCGCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCg                                                                        <  2:266090/150‑1 (MQ=255)
                                                                          cTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGGCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCgg                                                                       <  2:906372/150‑1 (MQ=255)
                                                                              cacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                                                                   >  2:372835/1‑150 (MQ=255)
                                                                              cacCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCa                                                                   >  2:1125335/1‑150 (MQ=255)
                                                                                 cGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCa                                                                <  2:466546/150‑1 (MQ=255)
                                                                                    aatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                                                             <  1:169237/150‑1 (MQ=255)
                                                                                    aatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                                                             <  1:173189/150‑1 (MQ=255)
                                                                                    aatttaattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCaaaa                                                             >  1:351050/1‑150 (MQ=255)
                                                                                        taattGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATAtt                                                         >  1:1001044/1‑150 (MQ=255)
                                                                                           atGCAACCGGCTCGTTGATGTGTTAATGACTAGTCCAGTCACGATCTTCATCAATATGCTGAAGCTTCGTGGCGGTATAACGATGATTTCTGATTGAATAAATATATTGCTGCATATAAATCTGCCTGAGCGGTGCATCAAAATATTTCa                                                      <  1:187597/149‑1 (MQ=255)
                                                                                              cAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCg                                                   >  1:1209940/1‑150 (MQ=255)
                                                                                                aCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTGCTGCATATACATCGGCCTGTGCGGTGCATCATACTTTTTCAGCGtt                                                 >  2:686824/1‑150 (MQ=255)
                                                                                                aCCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGtt                                                 >  2:13448/1‑150 (MQ=255)
                                                                                                  cGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGAGGTCCATCAAAATATTTCAGCGTTGc                                               >  1:414351/1‑150 (MQ=255)
                                                                                                       tgATGCTGGGTTAATGACTCGACCAGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGAt                                          <  1:544282/149‑1 (MQ=255)
                                                                                                       cGATGCTGGGTTAATGACTCGTCCCGTCACGCTCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCCTCAAAATATTTCAGCGTTGCCGGAt                                          >  1:212828/1‑150 (MQ=255)
                                                                                                          tGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATCCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAATGAATAAAGATTTTCGGGCATATACTTCGGACTGTGAGGTGCATCAAAATATTTCACCTTTGACGTATTCa                                       >  2:214803/1‑150 (MQ=255)
                                                                                                            cTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGATGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCAcc                                     <  2:338159/150‑1 (MQ=255)
                                                                                                                   aaTGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGcc                              >  2:639570/1‑150 (MQ=255)
                                                                                                                      gACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAgg                           <  1:813702/150‑1 (MQ=255)
                                                                                                                      gACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAgg                           <  2:351050/150‑1 (MQ=255)
                                                                                                                       aCTCGTCCCGTCACGATCTTCAGTAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATTAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGc                          <  2:909943/150‑1 (MQ=255)
                                                                                                                          cGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCaaa                       >  1:503404/1‑150 (MQ=255)
                                                                                                                          cGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCaaa                       >  2:759817/1‑150 (MQ=255)
                                                                                                                           gTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCaaaa                      >  1:174061/1‑150 (MQ=255)
                                                                                                                               cGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGCTTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCg                  >  1:788854/1‑150 (MQ=255)
                                                                                                                               cGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCg                  <  2:404863/150‑1 (MQ=255)
                                                                                                                                gTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCgg                 <  1:334312/150‑1 (MQ=255)
                                                                                                                                gTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCgg                 <  1:387246/150‑1 (MQ=255)
                                                                                                                                gTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCgg                 >  2:938248/1‑150 (MQ=255)
                                                                                                                                 ccACGATGCTTCAGCATATTGCTGAAGATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCggg                >  2:356797/2‑150 (MQ=255)
                                                                                                                                    cgatgcTTCAGCATATTGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGAt              >  1:77978/6‑150 (MQ=255)
                                                                                                                                      aTCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCAc           >  1:1229037/1‑150 (MQ=255)
                                                                                                                                      aTCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCAc           >  1:663967/1‑150 (MQ=255)
                                                                                                                                         cTTCAGCATATTGCTGAAGCATCGTGGCGGTAGAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGAGGTGCTTCAAAATTTTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCg         >  1:951135/1‑150 (MQ=255)
                                                                                                                                            cAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCa      <  1:387647/150‑1 (MQ=255)
                                                                                                                                            aaccattatGCTGAAGCATCGAGGCGGTATAACGATGATTTATGATTGAAAAAAGATATTCCTGCAGATACATCGGCGTGTGCGGCGCAGCAAACAATTTCAGCGTAGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCa      <  1:240077/144‑1 (MQ=255)
                                                                                                                                                aaTATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGaaa  >  2:841317/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAA  >  NC_000913/3066164‑3066456

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: