Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,897,913 T→C 8.1% intergenic (‑43/+51) ygcS ← / ← ygcU putative transporter YgcS/putative FAD‑containing dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,897,9130TC 8.1% 316.5 / 13.4 110intergenic (‑43/+51)ygcS/ygcUputative transporter YgcS/putative FAD‑containing dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (50/51);  new base C (3/6);  total (53/57)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 110 total aligned reads displayed.

CGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCG  >  NC_000913/2897765‑2898057
                                                                                                                                                    |                                                                                                                                                
cGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTg                                                                                                                                                 <  2:398596/150‑1 (MQ=255)
cGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTg                                                                                                                                                 >  1:1484989/1‑150 (MQ=255)
   gaACATAACCGTCGGTCAGGTGTGCGCAGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTgaga                                                                                                                                              <  1:1446160/150‑1 (MQ=255)
     aCATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAg                                                                                                                                            >  2:471446/1‑150 (MQ=255)
        tAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCg                                                                                                                                         <  1:1341052/150‑1 (MQ=255)
         aaCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGa                                                                                                                                        <  1:68542/150‑1 (MQ=255)
         aaCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgt                                                                                                                                        <  1:483249/150‑2 (MQ=255)
          aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGaa                                                                                                                                       <  1:53024/150‑1 (MQ=255)
          aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGaa                                                                                                                                       <  1:1260252/150‑1 (MQ=255)
          aCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgtt                                                                                                                                       <  2:1387630/150‑3 (MQ=255)
           ccGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGaaa                                                                                                                                      <  2:1301682/150‑1 (MQ=255)
           ccGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGaaa                                                                                                                                      >  2:179467/1‑150 (MQ=255)
            cGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGTTGAATTTTTCATCTTTACACCTTCATTTTCCGGTAATAAGAACTGCTTCTCAGGCGAGAAGCCgttct                                                                                                                                     <  2:1214810/150‑5 (MQ=255)
               cGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGc                                                                                                                                  <  2:1079390/150‑1 (MQ=255)
                   cAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTTAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCACCg                                                                                                                              >  2:889639/1‑150 (MQ=255)
                    aGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGc                                                                                                                             <  2:849438/150‑1 (MQ=255)
                        gtgCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCg                                                                                                                         <  2:1337387/150‑1 (MQ=255)
                        gtgCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCg                                                                                                                         <  2:551454/150‑1 (MQ=255)
                        gtgCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCg                                                                                                                         <  2:1200487/150‑1 (MQ=255)
                          gcgcCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGGAGCCGAAAAGGCAACGCCCCGTc                                                                                                                       >  2:127502/1‑150 (MQ=255)
                           cgcCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCa                                                                                                                      >  1:1443502/1‑150 (MQ=255)
                           cgcCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCa                                                                                                                      >  2:704895/1‑150 (MQ=255)
                           cgcCGAAAATGACCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGACCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCa                                                                                                                      >  2:1462598/1‑150 (MQ=255)
                             ccGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATc                                                                                                                    >  2:851841/1‑150 (MQ=255)
                               gAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATcgc                                                                                                                  <  1:760405/150‑1 (MQ=255)
                                aaaGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATTTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATcgcg                                                                                                                 <  1:240007/150‑1 (MQ=255)
                                  aGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGa                                                                                                               <  2:801483/150‑1 (MQ=255)
                                   gTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGaa                                                                                                              >  1:1535315/1‑150 (MQ=255)
                                        cgcAGCAATGCGGCAGTGAAAACGGTAAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTg                                                                                                         >  2:328530/1‑150 (MQ=255)
                                         gcagcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGc                                                                                                        >  2:1479843/1‑150 (MQ=255)
                                                   ggCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACatt                                                                                              <  1:1013691/150‑1 (MQ=255)
                                                   ggCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACatt                                                                                              >  2:1223323/1‑150 (MQ=255)
                                                        tGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTAtttt                                                                                         >  1:1321560/1‑150 (MQ=255)
                                                         gAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTAttttt                                                                                        >  1:1186863/1‑150 (MQ=255)
                                                          aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTc                                                                                       <  1:704384/150‑1 (MQ=255)
                                                          aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTc                                                                                       >  2:900587/1‑150 (MQ=255)
                                                             aCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGCAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAAt                                                                                    <  1:127502/150‑1 (MQ=255)
                                                              cGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATc                                                                                   >  2:1268249/1‑150 (MQ=255)
                                                                  tAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgat                                                                               >  2:195150/1‑150 (MQ=255)
                                                                   aaGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgata                                                                              <  2:261285/150‑1 (MQ=255)
                                                                    atcGGTAAATCATCCATTCTCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgatag                                                                             >  2:1004348/3‑150 (MQ=255)
                                                                    agtggTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgatag                                                                             <  2:1062075/147‑1 (MQ=255)
                                                                    aGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgatag                                                                             <  2:1053153/150‑1 (MQ=255)
                                                                    aGCGGTAAATCATCCATTCGCACCGGGGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGATCGGCTTCTCGCCTGAGAAGCCTAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgatag                                                                             <  2:1424550/150‑1 (MQ=255)
                                                                      cGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgatagat                                                                           <  1:893650/150‑1 (MQ=255)
                                                                           aaTCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTAc                                                                      >  2:1528737/1‑150 (MQ=255)
                                                                            atcatcCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTAcc                                                                     <  1:253959/150‑1 (MQ=255)
                                                                              catcCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCg                                                                   >  2:1031087/1‑150 (MQ=255)
                                                                              catcCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCg                                                                   <  2:668818/150‑1 (MQ=255)
                                                                                 ccATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTAt                                                                <  2:538694/150‑1 (MQ=255)
                                                                                      cGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGac                                                                                                   >  2:412062/1‑110 (MQ=255)
                                                                                      cGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGac                                                                                                   <  1:412062/110‑1 (MQ=255)
                                                                                      cGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATa                                                           >  2:993336/1‑150 (MQ=255)
                                                                                      cGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATa                                                           >  2:137837/1‑150 (MQ=255)
                                                                                      cGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTCCAATCgg                                                                                 <  2:1184538/128‑1 (MQ=255)
                                                                                      cGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTCCAATCgg                                                                                 >  1:1184538/1‑128 (MQ=255)
                                                                                          ccGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGc                                                       <  2:1360987/150‑1 (MQ=255)
                                                                                          ccGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGc                                                       >  1:377690/1‑150 (MQ=255)
                                                                                               gAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTa                                                  <  1:1147925/150‑1 (MQ=255)
                                                                                                aaGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAg                                                 <  1:689973/150‑1 (MQ=255)
                                                                                                    gtTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATc                                             <  1:704895/150‑1 (MQ=255)
                                                                                                    gtTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATc                                             <  1:1054852/150‑1 (MQ=255)
                                                                                                    gtTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATc                                             <  1:195150/150‑1 (MQ=255)
                                                                                                     ttCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCg                                            <  1:1031087/150‑1 (MQ=255)
                                                                                                     ttCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCg                                            <  1:851841/150‑1 (MQ=255)
                                                                                                     ttCATCTTTACCACTTCATTGTACGGTAATAAGAACGGCTTCTCGACTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCg                                            >  2:54837/1‑150 (MQ=255)
                                                                                                          cTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCTCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTg                                       >  1:770811/1‑150 (MQ=255)
                                                                                                             tACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCtt                                    >  2:713382/1‑150 (MQ=255)
                                                                                                              aCCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCttt                                   >  1:1533468/1‑150 (MQ=255)
                                                                                                               ccccTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCtttt                                  >  1:140865/1‑150 (MQ=255)
                                                                                                               ccccTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCtttt                                  >  2:25876/1‑150 (MQ=255)
                                                                                                                 ccTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTc                                <  1:1150044/150‑1 (MQ=255)
                                                                                                                 ccTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAACCCGAAAACCCAACGCCCCGTCATAGCGCAACGCAGCCTGCCACATTATTTTTCAATCTCATAGATACTACCCGTATTCATAAAGCCATTAGGATCGAAATGCttatta                                >  2:431446/1‑149 (MQ=255)
                                                                                                                 aattCAATGTCAGGTAATAAGAACGGCTTCTCGCCTGAGAAGACGAAAAGGCAACGCACCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTc                                >  2:487132/3‑150 (MQ=255)
                                                                                                                  cTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCa                               <  2:568864/150‑1 (MQ=255)
                                                                                                                   ttCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAg                              <  2:1443502/150‑1 (MQ=255)
                                                                                                                    tCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCGACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGc                             >  1:399421/1‑150 (MQ=255)
                                                                                                                    tCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGc                             <  2:136034/150‑1 (MQ=255)
                                                                                                                     cATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGcc                            <  1:1223609/150‑1 (MQ=255)
                                                                                                                     cATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGcc                            >  2:1489022/1‑150 (MQ=255)
                                                                                                                      aTTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAGCTGCTTTTTCAGccc                           >  2:561792/1‑150 (MQ=255)
                                                                                                                          tCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTcc                       >  2:801719/1‑150 (MQ=255)
                                                                                                                           ccGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCa                      >  2:653181/1‑150 (MQ=255)
                                                                                                                            cGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAg                     >  1:295173/1‑150 (MQ=255)
                                                                                                                                aataaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAAc                 >  1:117745/1‑150 (MQ=255)
                                                                                                                                 ataaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACg                <  2:1179536/150‑1 (MQ=255)
                                                                                                                                 ataaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAAAGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACg                >  2:1231214/1‑150 (MQ=255)
                                                                                                                                  taaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGc               >  1:1057582/1‑150 (MQ=255)
                                                                                                                                  ccttttcGGCTTCTCAGGCGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGc               >  1:78246/7‑150 (MQ=255)
                                                                                                                                   cttttcGGCTTCTCAGGCGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGcc              >  2:1376638/6‑150 (MQ=255)
                                                                                                                                    aGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGccc             >  2:756380/1‑150 (MQ=255)
                                                                                                                                    aGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGccc             >  1:1432378/1‑150 (MQ=255)
                                                                                                                                     gAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCa            <  2:1184118/150‑1 (MQ=255)
                                                                                                                                      ttcGGCTTCTCAGGCGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAg                     <  1:1529630/138‑1 (MQ=255)
                                                                                                                                       aCGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCAcg          <  2:1527998/150‑1 (MQ=255)
                                                                                                                                       aCGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCAcg          <  1:989924/150‑1 (MQ=255)
                                                                                                                                        cGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCAcga         <  1:1004348/150‑2 (MQ=255)
                                                                                                                                          gCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACgcgc       >  2:858722/1‑150 (MQ=255)
                                                                                                                                           cTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCt      >  2:885738/1‑150 (MQ=255)
                                                                                                                                               tcGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCg  >  2:61847/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCG  >  NC_000913/2897765‑2898057

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