Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,985,112 C→G 6.8% intergenic (‑43/+27) araG ← / ← araF L‑arabinose ABC transporter ATPase/L‑arabinose ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,985,1120CG 6.8% 297.9 / 8.1 97intergenic (‑43/+27)araG/araFL‑arabinose ABC transporter ATPase/L‑arabinose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (46/43);  new base G (4/3);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTAAA  >  NC_000913/1984968‑1985260
                                                                                                                                                |                                                                                                                                                    
cTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCa                                                                                                                                                 >  1:1698913/1‑150 (MQ=255)
 tGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCag                                                                                                                                                <  2:269603/150‑1 (MQ=255)
     cGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggg                                                                                                                                            <  1:555339/150‑1 (MQ=255)
       gCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATGATGCAGAggggg                                                                                                                                          <  2:1329173/150‑2 (MQ=255)
          tAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAAtt                                                                                                                                       >  2:1409472/1‑150 (MQ=255)
            gCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAAtttt                                                                                                                                     >  2:1018358/1‑150 (MQ=255)
             cAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACTGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAAttttt                                                                                                                                    <  1:1403732/150‑1 (MQ=255)
              aGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTc                                                                                                                                   >  2:1368847/1‑150 (MQ=255)
               gTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGATACGTTTTACCGATGCCGCGAAATGAGAGAGACGGGGGAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTcc                                                                                                                                  <  2:1588281/150‑1 (MQ=255)
                tCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCg                                                                                                                                 <  1:1361415/150‑1 (MQ=255)
                 cAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAACCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACCCCCCCGCTGCATCATGCAGAGGGGATGTTTTCCgg                                                                                                                                >  2:328190/1‑150 (MQ=255)
                  aaaaCTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGCCTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGc                                                                                                                               >  1:614564/1‑150 (MQ=255)
                  aaaaCTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGc                                                                                                                               <  2:1369045/150‑1 (MQ=255)
                      cTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAt                                                                                                                           >  2:726169/1‑150 (MQ=255)
                          tatCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAttact                                                                                                                       >  1:1381600/1‑150 (MQ=255)
                          tatCCGTCAGCGCCTTAACGCCCGGAAACGTTTTACCGATTCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAAttact                                                                                                                       <  2:526494/150‑1 (MQ=255)
                                 cAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGcc                                                                                                                >  1:1492161/1‑150 (MQ=255)
                                   gcgcCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATAATGAAGAGGGGATTTTTTCCGGCAAATTAATTACCGCATa                                                                                                              >  2:938132/1‑150 (MQ=255)
                                      ccTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCCGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAAc                                                                                                           <  1:470229/150‑1 (MQ=255)
                                            cGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCttttt                                                                                                     >  2:1027868/1‑150 (MQ=255)
                                              ccGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtc                                                                                                   <  2:64652/150‑1 (MQ=255)
                                               cGGGTAACGTTTTAACGATGCAGCGAAATGAGAGATACGGTGTAGACTGTTGCTTAAAGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtct                                                                                                  <  2:1613669/150‑1 (MQ=255)
                                               cGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtct                                                                                                  <  1:34227/150‑1 (MQ=255)
                                                gggTAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTTTTGGATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtctc                                                                                                 <  1:938132/150‑1 (MQ=255)
                                                gggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtctc                                                                                                 >  2:1569528/1‑150 (MQ=255)
                                                gggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTtctc                                                                                                 >  1:1476217/1‑150 (MQ=255)
                                                 ggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTcc                                                                                                <  2:1523068/150‑1 (MQ=255)
                                                 ggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTcc                                                                                                <  2:1765822/150‑1 (MQ=255)
                                                 ggAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTcc                                                                                                <  2:1521889/150‑1 (MQ=255)
                                                  tataCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCa                                                                                               <  2:1812011/147‑1 (MQ=255)
                                                  gAAACGTTTTCCCGATTCCGCGTAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCa                                                                                               <  2:888864/150‑1 (MQ=255)
                                                    aaCGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGt                                                                                             >  2:190898/1‑150 (MQ=255)
                                                      cGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttc                                                                                           <  1:488596/150‑1 (MQ=255)
                                                       gTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttct                                                                                          >  2:645120/1‑150 (MQ=255)
                                                       gTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGACTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttct                                                                                          >  2:965725/1‑150 (MQ=255)
                                                         tttACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGttcttc                                                                                        >  1:603927/1‑150 (MQ=255)
                                                           tACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTcttctt                                                                                      <  2:58824/150‑1 (MQ=255)
                                                           tACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAAACTTTTTTCTCCAGTTcttctt                                                                                      >  2:1298446/1‑150 (MQ=255)
                                                            aCCGATGCCGCGAATTGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTCGTCGTTCACACCCCCTCTGCATCTTTCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTGTCCAGTTCTTCttt                                                                                     <  2:1826047/150‑1 (MQ=255)
                                                            aCCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCttt                                                                                     >  1:564501/1‑150 (MQ=255)
                                                            aCCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCttt                                                                                     >  2:559148/1‑150 (MQ=255)
                                                            aCCGATGCCGCGAAATGAGAGATACGGGGGAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCttt                                                                                     <  2:164160/150‑1 (MQ=255)
                                                             ccGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTcttctt                                                                                      <  2:1758355/148‑1 (MQ=255)
                                                             ccGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTcttctt                                                                                      >  1:1758355/1‑148 (MQ=255)
                                                             ccGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTa                                                                                    >  2:477925/1‑150 (MQ=255)
                                                                aTGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAg                                                                                 <  1:518893/150‑1 (MQ=255)
                                                                  gCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGtt                                                                               <  2:884397/150‑1 (MQ=255)
                                                                   ccGCGAAATGAGAGATACGGGGGAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCTTTCACACCCCCTCTGCATCATGCAGAGGTGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTg                                                                              <  1:43418/150‑1 (MQ=255)
                                                                    cgcgAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGt                                                                             >  2:725073/1‑150 (MQ=255)
                                                                          aTGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCAcgcg                                                                       <  2:1826634/150‑1 (MQ=255)
                                                                           tGAGAGATACGGGGGAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGt                                                                      <  2:1476704/150‑1 (MQ=255)
                                                                              gagaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGAt                                                                   >  1:24688/1‑150 (MQ=255)
                                                                               agaTACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATc                                                                  <  1:1018358/150‑1 (MQ=255)
                                                                                  tACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGt                                                               <  1:377030/150‑1 (MQ=255)
                                                                                     ggggTAGACTGTCGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTAcc                                                            <  2:817990/150‑1 (MQ=255)
                                                                                          aGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTc                                                       >  1:1801863/1‑150 (MQ=255)
                                                                                          aGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTc                                                       <  1:1811800/150‑1 (MQ=255)
                                                                                            aCTGTTGCATAACGTCTCCTTGAATCACGGGTCGTTAACACCCCCTCTGGATGAAGCTCAGGGGAATTTTTCCGGCAAATTACTTACCGCCTTAACCTTTTTTCTCCAGTTCTTCTTTAAACTTGTCACGCGTGATCAGTTCCACGTggg                                                     <  2:724204/150‑1 (MQ=255)
                                                                                            aCTGTTGCATAACGTCTCCGTGATTCAGTGGTCGTTCACACCCCCTCTGCTTCATGAAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCgg                                                     <  1:708125/150‑1 (MQ=255)
                                                                                            aCTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCAACATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCgg                                                     <  2:1607005/150‑1 (MQ=255)
                                                                                              tGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTa                                                   <  1:1494767/150‑1 (MQ=255)
                                                                                              tGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTa                                                   <  2:20966/150‑1 (MQ=255)
                                                                                               gTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCGTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTaa                                                  >  1:1071977/1‑150 (MQ=255)
                                                                                                 tGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACt                                                >  2:1700180/1‑150 (MQ=255)
                                                                                                   cATAACGTCTCCGTGAATCACTGGTCGTTGACCCCCCCTCTGCAGCATGCAGAGGGGCATTTTTCCGGCAAATTACTTACCGACTAATCCTTTTTTAGGCAGGTAGTCTTTTAAGTTGTCACGCGGGGTCAGTTCCCCGTCgggtacgtc                                              >  2:906670/1‑143 (MQ=255)
                                                                                                           ctcCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAAtt                                      >  2:1305733/1‑150 (MQ=255)
                                                                                                             tcGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAAtttt                                    >  2:23805/2‑150 (MQ=255)
                                                                                                             tcGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAAtttt                                    >  2:426809/2‑150 (MQ=255)
                                                                                                              cGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTg                                   >  1:689614/1‑150 (MQ=255)
                                                                                                                tGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGc                                 <  1:494897/150‑1 (MQ=255)
                                                                                                                 gAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCg                                >  2:308978/1‑150 (MQ=255)
                                                                                                                   aTCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGt                              >  2:1361900/1‑150 (MQ=255)
                                                                                                                    tCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGtt                             <  2:1596821/150‑1 (MQ=255)
                                                                                                                     cACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTc                            >  1:1156380/1‑150 (MQ=255)
                                                                                                                      aCTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCa                           <  1:1409472/150‑1 (MQ=255)
                                                                                                                        tGGTCGTTCACACCCCCTCTGCATCATGCAGTGGGGTATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAAc                         <  1:682409/150‑1 (MQ=255)
                                                                                                                        tGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGATCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAAc                         <  2:1013627/150‑1 (MQ=255)
                                                                                                                         ggTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACg                        <  2:736148/150‑1 (MQ=255)
                                                                                                                          gTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGt                       >  2:1590839/1‑150 (MQ=255)
                                                                                                                            cGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGGGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCt                     >  2:1426354/1‑150 (MQ=255)
                                                                                                                               tCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCtttt                  >  2:999433/1‑150 (MQ=255)
                                                                                                                               tCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCtttt                  <  1:1201164/150‑1 (MQ=255)
                                                                                                                               tCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCtttt                  >  1:192886/1‑150 (MQ=255)
                                                                                                                               tCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCtttt                  >  1:491977/1‑150 (MQ=255)
                                                                                                                                 acacCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGACGGTAACTTCGGTAAATTTTGGCGGGTCAACGTCTTTGGc                >  1:1052571/1‑150 (MQ=255)
                                                                                                                                 aattccccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCATTTgg                >  1:1779504/5‑149 (MQ=255)
                                                                                                                                    cccccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGGACCACGTCGGTATATGCGGTAAATTTTGGCGGTTCAACGTCTTGTGCTAc             >  1:1357566/1‑150 (MQ=255)
                                                                                                                                     ccccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAAAGTCTTTTGCTAcc            >  2:1321598/1‑150 (MQ=255)
                                                                                                                                     ccccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAcc            >  1:872989/1‑150 (MQ=255)
                                                                                                                                     ccccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAccc           >  2:642755/1‑150 (MQ=255)
                                                                                                                                      cccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAACCCTTTTTTATCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAAGTTTGGAGGTTCAAAGTCTTTTGCTAccc           >  1:859571/1‑150 (MQ=255)
                                                                                                                                      cccTCTGCATGATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAccc           >  1:859565/1‑150 (MQ=255)
                                                                                                                                      cccTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTAccc           >  2:754354/1‑150 (MQ=255)
                                                                                                                                        ctctGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAg         <  2:1698913/150‑1 (MQ=255)
                                                                                                                                           ttcataatGGAGAGGGGAATTTTTCCGGCAAATTTCTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTg      <  2:775343/148‑1 (MQ=255)
                                                                                                                                           tGCATCATGCAGAGGGGAATTTTTCCGGCAAATTTCTTACCGCCTAAACCTTTTTTCTCCAGTTCATCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAATTTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTg      <  2:1216486/150‑1 (MQ=255)
                                                                                                                                           tGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTg      >  2:1461465/1‑150 (MQ=255)
                                                                                                                                            gCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTACACGTTTTTTCTCCAGTTCTTCTTTAAAGTTGTAAAGCGTCATCTGTACCACGGCGGTATCTTCGGTAAATTTTTGCGGTTCAACGTCTTTTGCTACACAGTTGt     >  2:79279/1‑150 (MQ=255)
                                                                                                                                               tcatGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTaaa  <  1:190898/150‑1 (MQ=255)
                                                                                                                                               tcatGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTaaa  >  2:661226/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CTGACCGGCATAGCAGTCAAAACTAATATCCGTCAGCGCCTTAACGCCGGGAAACGTTTTACCGATGCCGCGAAATGAGAGATACGGGGTAGACTGTTGCATAACGTCTCCGTGAATCACTGGTCGTTCACACCCCCTCTGCATCATGCAGAGGGGAATTTTTCCGGCAAATTACTTACCGCCTAAACCTTTTTTCTCCAGTTCTTCTTTAAAGTTGTCACGCGTGATCAGTACCACGTCGGTAACTTCGGTAAATTTTGGCGGTTCAACGTCTTTTGCTACCCAGTTGTAAA  >  NC_000913/1984968‑1985260

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: