Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 144,812 T→G 5.5% L79* (TTA→TGA)  yadI → putative PTS Enzyme IIA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913144,8120TG 5.5% 287.5 / 7.6 110L79* (TTA→TGA) yadIputative PTS Enzyme IIA
Reads supporting (aligned to +/- strand):  ref base T (37/67);  new base G (4/2);  total (41/69)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.73e-01

Only 100 of 114 total aligned reads displayed.

GTGCCGGGCCGTGAATTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTT  >  NC_000913/144663‑144957
                                                                                                                                                     |                                                                                                                                                 
     gggCCGTGAATTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTAATAATGAAGcc                                                                                                                                              <  1:664148/150‑1 (MQ=255)
       gCCGTGAATTTCTGGCGCGGATTAAGCTATAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCcac                                                                                                                                            <  2:1691948/150‑1 (MQ=255)
       gCCGAGAATTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACAGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCcac                                                                                                                                            >  1:1406800/1‑150 (MQ=255)
         cGTGAATTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTAATAATGAAGCcacgc                                                                                                                                          <  1:669806/150‑3 (MQ=255)
         cGTGAATTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGAAATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACaa                                                                                                                                          >  2:311078/1‑150 (MQ=255)
           tGAATTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAac                                                                                                                                        <  1:459198/150‑1 (MQ=255)
              aTTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACt                                                                                                                                     >  1:1353296/1‑150 (MQ=255)
                ttCTGGCGCGGATTAAGCGCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTcc                                                                                                                                   <  1:1194368/150‑1 (MQ=255)
                      cgcgGATTAAGCGCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGc                                                                                                                             <  1:1548399/150‑1 (MQ=255)
                         ggATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTTAGCCAAAAACACTCCCGTTGCGaa                                                                                                                          >  2:666892/1‑150 (MQ=255)
                         ggATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGAGaa                                                                                                                          >  2:1445048/1‑150 (MQ=255)
                         ggATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATAGCGTGTCTTCATTATTAAGCAACAAACACTCCCGTTGCGaa                                                                                                                          >  2:672938/1‑150 (MQ=255)
                            ttAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTg                                                                                                                       <  1:1727875/150‑1 (MQ=255)
                             tAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGTGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACATACACTCCCGTTGCGAAGTGa                                                                                                                      <  2:735895/150‑1 (MQ=255)
                              aaGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGAt                                                                                                                     <  1:1480479/150‑1 (MQ=255)
                               aGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGAtt                                                                                                                    >  1:1335898/1‑150 (MQ=255)
                                gCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGAttt                                                                                                                   >  1:91102/1‑150 (MQ=255)
                                   ctAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCAGCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTg                                                                                                                <  2:1286267/150‑1 (MQ=255)
                                    tAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTgg                                                                                                               <  1:753422/150‑1 (MQ=255)
                                    tAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTgg                                                                                                               <  1:1606900/150‑1 (MQ=255)
                                       tatGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGTGGAATCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGgtgt                                                                                                            <  2:1061043/150‑1 (MQ=255)
                                       tatGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGgtgt                                                                                                            >  1:827429/1‑150 (MQ=255)
                                         tGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCa                                                                                                          >  2:928999/1‑150 (MQ=255)
                                           cgcAGCCGCATGATGGGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACg                                                                                                        <  1:406238/148‑1 (MQ=255)
                                               gCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTAc                                                                                                    >  2:6836/1‑150 (MQ=255)
                                               gCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTAc                                                                                                    <  1:1062333/150‑1 (MQ=255)
                                                ccGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTAcc                                                                                                   <  1:1257335/150‑1 (MQ=255)
                                                 cGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCg                                                                                                  <  2:1353296/150‑1 (MQ=255)
                                                 cGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCg                                                                                                  <  2:299425/150‑1 (MQ=255)
                                                  gCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACACCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGt                                                                                                 >  2:357738/1‑150 (MQ=255)
                                                   cATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGtt                                                                                                <  2:39682/150‑1 (MQ=255)
                                                   cATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGtt                                                                                                <  2:971689/150‑1 (MQ=255)
                                                     tgatgTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTaa                                                                                              >  2:1086991/1‑150 (MQ=255)
                                                      gatgTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGTGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAAt                                                                                             <  2:709506/150‑1 (MQ=255)
                                                         gtgCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGa                                                                                          <  1:331531/150‑1 (MQ=255)
                                                           gCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAAc                                                                                        <  1:274707/150‑1 (MQ=255)
                                                             gATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAg                                                                                      <  2:1592247/150‑1 (MQ=255)
                                                              aTGCTCTGCATGTAGCTGAGTCTGGTTAGGCTGTCATCTTCTTAACCGATATAGACGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCAACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAga                                                                                     <  1:258210/150‑1 (MQ=255)
                                                               tGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAgat                                                                                    <  1:1707009/150‑1 (MQ=255)
                                                                   cttcATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATATCCGGAGCGCCACCGTTTCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGatgatg                                                                                <  1:311078/147‑1 (MQ=255)
                                                                      cATGAAGTGGACTCTGGTGAGGGTGTTATCTTCTTTAGCGATATAGCCGGAGCGCATACGTATAGTGTTGCTTCATTTTTAAGCTACAAACACTCCCGTTGCGTAGTGATTTCTGGTGTCACGTTAACGTTAATTGAACAGATGATGGCt                                                                             <  2:1116976/150‑1 (MQ=255)
                                                                        tGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCTGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTg                                                                           <  2:669528/150‑1 (MQ=255)
                                                                         gAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGc                                                                          <  1:1049280/150‑1 (MQ=255)
                                                                         gAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGc                                                                          <  1:457796/150‑1 (MQ=255)
                                                                         gAAGCGGACTCTCGTGAGGGTGTCATCTTCTAAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATAGAACAGATGATGTCTTGc                                                                          >  2:1528052/1‑150 (MQ=255)
                                                                          aaGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGcc                                                                         <  2:1212622/150‑1 (MQ=255)
                                                                           aGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCg                                                                        <  1:1664254/150‑1 (MQ=255)
                                                                           aGCGGACTCTGGTGAGGGTCGCATCTTCTTAAACGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACTCTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCg                                                                        <  1:263088/150‑1 (MQ=255)
                                                                                 ctctGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAAc                                                                  >  2:790954/1‑150 (MQ=255)
                                                                                  tctGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAAcc                                                                 <  2:1340526/150‑1 (MQ=255)
                                                                                  tctGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAAcc                                                                 <  2:857037/150‑1 (MQ=255)
                                                                                   cttgTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTCTCGCGGGGCTTCATTATTACGCCACAAACACTCCCTTTTCGAATTGATTTCTGGCGTCACGCTACCGTTAATGGAACAGATGCGGGCCTGCCGTGAAACCa                                                                <  2:289298/147‑1 (MQ=255)
                                                                                    tGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCAt                                                               <  1:1586440/150‑1 (MQ=255)
                                                                                     ggTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATg                                                              <  2:749300/150‑1 (MQ=255)
                                                                                             tgtCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTc                                                      <  2:871658/150‑1 (MQ=255)
                                                                                               tCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACGATTTCat                                                    >  1:1075565/1‑149 (MQ=255)
                                                                                               gcATCTTCTTAACCGATCTAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGGGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCag                                                    <  2:297091/149‑1 (MQ=255)
                                                                                                 aTCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCagag                                                  <  1:1545428/150‑1 (MQ=255)
                                                                                                 aTCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCagag                                                  <  1:1098556/150‑1 (MQ=255)
                                                                                                  tcttctTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTGGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGt                                                 <  2:91102/150‑1 (MQ=255)
                                                                                                   cttcttAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGtt                                                <  1:176460/150‑1 (MQ=255)
                                                                                                     tcttAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTc                                              >  2:627288/1‑150 (MQ=255)
                                                                                                     tcttAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTc                                              <  2:1518901/150‑1 (MQ=255)
                                                                                                       ttAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTcgc                                            <  1:709965/150‑1 (MQ=255)
                                                                                                        tAACCGATATAGCCGGATCGCCACCGTATCGCGTGGCTGCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGGCACGTTACCGTTAATTGAACAGATGATGGCGTGCCGTGAAACCATGACCAGTTCAGAGTTTcgcg                                           <  1:1047881/150‑1 (MQ=255)
                                                                                                        tAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCATCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTcgcg                                           >  2:724734/1‑150 (MQ=255)
                                                                                                         aaCCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGa                                          >  2:1116927/1‑150 (MQ=255)
                                                                                                          aCCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAACCCATGACCAGTTCAGAGTTTCGCGAg                                         >  2:233813/1‑150 (MQ=255)
                                                                                                          aCCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAg                                         <  2:1472178/150‑1 (MQ=255)
                                                                                                          aCCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAg                                         <  2:1186169/150‑1 (MQ=255)
                                                                                                           ccGATATAGCCGGAGCGACACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATAGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGc                                        >  2:1724038/1‑150 (MQ=255)
                                                                                                              atatAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTa                                     <  1:811227/150‑1 (MQ=255)
                                                                                                               tataGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTAt                                    <  2:569659/150‑1 (MQ=255)
                                                                                                               tataGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTAt                                    >  1:476249/1‑150 (MQ=255)
                                                                                                                  aGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGt                                 >  2:1686494/1‑150 (MQ=255)
                                                                                                                  aGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGt                                 >  1:1437241/1‑150 (MQ=255)
                                                                                                                   gCCGGAGCGCCACCGTATCGCTTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTc                                <  2:1524237/150‑1 (MQ=255)
                                                                                                                    ccGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCg                               <  1:998822/150‑1 (MQ=255)
                                                                                                                    ccGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGAGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCg                               >  2:792295/1‑150 (MQ=255)
                                                                                                                     cGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTTCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGa                              <  1:1329265/150‑1 (MQ=255)
                                                                                                                          cgcCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGATCAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTgg                         >  2:166502/1‑150 (MQ=255)
                                                                                                                             cACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTg                      <  1:1086991/150‑1 (MQ=255)
                                                                                                                              aCCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCTAGCGTATTGTCGAACTGGGTgc                     >  2:1428756/1‑150 (MQ=255)
                                                                                                                              aCCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTgc                     <  1:159394/150‑1 (MQ=255)
                                                                                                                                 gTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGcc                  <  1:700464/150‑1 (MQ=255)
                                                                                                                                   aTCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCAGGTGTCACGTTACCGTTAATTGAACAGATGATGCCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCTAGCGTATTGTCGAACAGGGTGCGCCgg                >  1:1274660/1‑150 (MQ=255)
                                                                                                                                    tCGCGTGGCTTCATTATTAAGCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAg              <  2:959196/150‑1 (MQ=255)
                                                                                                                                    gtttgTGGCTTAATAATGAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAACCCAAGACCAGTTCAGACTTTCGCGAGCGTATTGTCGAACAGGGTGCGCCGGa               >  2:786260/5‑150 (MQ=255)
                                                                                                                                    gtttgTGGCTTAATAATGAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGa               >  1:1172775/5‑150 (MQ=255)
                                                                                                                                     tttgTGGCTTAATAATGAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAg              >  1:968515/4‑150 (MQ=255)
                                                                                                                                      gcgTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAgg             >  2:1465540/1‑150 (MQ=255)
                                                                                                                                      gcgTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAgg             >  1:144174/1‑150 (MQ=255)
                                                                                                                                       tgTGGCTTAATAATGAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGt            >  1:399691/2‑150 (MQ=255)
                                                                                                                                        gTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTg           <  1:365499/150‑1 (MQ=255)
                                                                                                                                            cTTCAATATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTAACGTTAATTGAACAGATGATGGATTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCTAGCGTATTGTCGAACAGGGTGCGCCGGAGGTCagta       >  1:47958/1‑150 (MQ=255)
                                                                                                                                             ttCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGagtag      <  1:851093/150‑1 (MQ=255)
                                                                                                                                              tCATTATTAAGCCACAAACACTCCCGTTTCGAACACATAACTGGAGTCACGTTACCGTTAATTGAAAAGATGATCGCTAGCCGTGAAACCATGACCAGTTCAGAGTATCGCGAGCCAATTGTCGAACTGGGTGCGCCGGAGGTGagtagt     >  2:1239997/1‑150 (MQ=255)
                                                                                                                                               cATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTc    >  1:955400/1‑150 (MQ=255)
                                                                                                                                                attattAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCt   <  2:1069905/150‑1 (MQ=255)
                                                                                                                                                 ttattaAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGGCACGGTACCTTTAATTGAACAGATGATGGCTTGCCGGGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCtt  <  2:861954/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GTGCCGGGCCGTGAATTTCTGGCGCGGATTAAGCTCTAATATGCTCAGCCGCATGATGTGCGATGCTCTGCATGAAGCGGACTCTGGTGAGGGTGTCATCTTCTTAACCGATATAGCCGGAGCGCCACCGTATCGCGTGGCTTCATTATTAAGCCACAAACACTCCCGTTGCGAAGTGATTTCTGGTGTCACGTTACCGTTAATTGAACAGATGATGGCTTGCCGTGAAACCATGACCAGTTCAGAGTTTCGCGAGCGTATTGTCGAACTGGGTGCGCCGGAGGTGAGTAGTCTT  >  NC_000913/144663‑144957

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