Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,066,309 T→A 6.7% I289F (ATT→TTT) ‡ 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,3090TA 6.7% 227.9 / 5.4 86I289F (ATT→TTT) ‡ygfIputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (42/38);  new base A (3/3);  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 = 5.80e-01

GTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAA  >  NC_000913/3066160‑3066456
                                                                                                                                                      |                                                                                                                                                   
gTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAat                                                                                                                                                     >  1:1069710/1‑150 (MQ=255)
    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)
                                                                                                                                                      |                                                                                                                                                   
GTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGCATATACATCGGCCTGTGCGGTGCATCAAAATATTTCAGCGTTGCCGGATTCACCAGTCGCCAGGCAAAATGCGGGATCACCGTCAGAAA  >  NC_000913/3066160‑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: