Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,052,257 C→T 6.0% intergenic (+17/+32) gnsA → / ← yccM putative phosphatidylethanolamine synthesis regulator/putative 4Fe‑4S membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,052,2570CT 6.0% 363.5 / 7.8 117intergenic (+17/+32)gnsA/yccMputative phosphatidylethanolamine synthesis regulator/putative 4Fe‑4S membrane protein
Reads supporting (aligned to +/- strand):  ref base C (50/59);  new base T (4/3);  total (55/62)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 120 total aligned reads displayed.

GCCGACTTTATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGA  >  NC_000913/1052109‑1052406
                                                                                                                                                    |                                                                                                                                                     
gCCGACTTTATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCAcccc                                                                                                                                                      >  2:1848386/1‑150 (MQ=255)
    aCTTTATCGCTCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGc                                                                                                                                                  <  1:2109133/150‑1 (MQ=255)
    aCTTTATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGc                                                                                                                                                  <  1:1189140/150‑1 (MQ=255)
          tCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTg                                                                                                                                            <  2:1164243/150‑1 (MQ=255)
           cgcgCAAAAAATCGCCGAGCTGAACAAGAATACAGGCAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATCACCGGGCTTGAAAGTAACGATCTCAAAAACTAAATAAGGTGTTTTTGTGAACATCACCACGTGCGAGgaga                                                                                                                                           >  2:613314/1‑150 (MQ=255)
             cgcAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATg                                                                                                                                         >  1:838146/1‑150 (MQ=255)
              tcAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGt                                                                                                                                        <  1:999309/149‑1 (MQ=255)
              gcAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGt                                                                                                                                        <  1:1541835/150‑1 (MQ=255)
              gcAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAAAAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGt                                                                                                                                        <  2:1913656/150‑1 (MQ=255)
                aaaaaaTCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCTTATTCAAACTAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGTTc                                                                                                                                      <  1:332640/150‑1 (MQ=255)
                 aaaaaTCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTcc                                                                                                                                     <  2:387548/150‑1 (MQ=255)
                  aaaaTCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGTTcac                                                                                                                                    <  2:715321/150‑3 (MQ=255)
                    aaTCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGACAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCt                                                                                                                                  >  2:1346289/1‑150 (MQ=255)
                       cGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCttgt                                                                                                                               <  2:672937/150‑1 (MQ=255)
                        gCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCTCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCttgtt                                                                                                                              <  1:716870/150‑1 (MQ=255)
                         ccGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCttgttg                                                                                                                             <  1:695853/150‑1 (MQ=255)
                          cGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGc                                                                                                                            >  2:1438939/1‑150 (MQ=255)
                          cGAGCTGAACAAGAATACAGGGAAAGAAGGCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGc                                                                                                                            <  2:1416610/150‑1 (MQ=255)
                           gAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCt                                                                                                                           <  2:2195019/150‑1 (MQ=255)
                               tGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAAtt                                                                                                                       <  1:460498/150‑1 (MQ=255)
                                gAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTTATTTAGTGAACATCACCCCGTGCGAGGTGAGGTTCCGCTTGTTGCTAAttt                                                                                                                      >  2:1445937/1‑150 (MQ=255)
                                gAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAAttt                                                                                                                      >  2:576486/1‑150 (MQ=255)
                                gAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAAttt                                                                                                                      <  2:1633999/150‑1 (MQ=255)
                                  aCAAGAATACAGTGAAAGAAGTATCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATAAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAg                                                                                                                    <  2:1412008/150‑1 (MQ=255)
                                  aCAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAg                                                                                                                    >  1:2052998/1‑150 (MQ=255)
                                  aCAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTGATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTAGCTAATTTat                                                                                                                    >  1:1969256/1‑149 (MQ=255)
                                   cAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATCACCGGGCTTGAAAGTAATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGt                                                                                                                   <  2:1117546/150‑1 (MQ=255)
                                   aaaGAATACATGGAAAGAAGTCTGTGAAATTCGCTTCACGGCACGAGAAAAAATGACCGGGCTTGAATGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGt                                                                                                                   <  2:1641352/149‑1 (MQ=255)
                                     aGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGa                                                                                                                 >  1:1714975/1‑150 (MQ=255)
                                       aaTACATGGAAAGAAGTCTCTGAAATTCACTTCACCGCACGAGAAAAAAAGACCGGCCTTGAAAGTTATGATGTCAAATTCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGAcc                                                                                                               <  2:381841/150‑1 (MQ=255)
                                          aCAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAAt                                                                                                            >  2:1986809/1‑150 (MQ=255)
                                          aCAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAAt                                                                                                            <  2:210200/150‑1 (MQ=255)
                                           cAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTACCGCTTGTTGCTAAATTAGTGACCAATc                                                                                                           >  2:793510/1‑150 (MQ=255)
                                           cAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACACCGTGCGAGGTGATGTTCCGCTTGTTGCTATTTTAGTGACCAATc                                                                                                           >  2:129132/1‑150 (MQ=255)
                                            aGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCa                                                                                                          >  1:1751174/1‑150 (MQ=255)
                                            aGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCa                                                                                                          <  2:2115854/150‑1 (MQ=255)
                                              ggAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATTATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCAtt                                                                                                        >  1:1631991/1‑150 (MQ=255)
                                               gAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTg                                                                                                       >  2:55392/1‑150 (MQ=255)
                                                aaaGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTgg                                                                                                      <  2:2090549/150‑1 (MQ=255)
                                                aaaGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTgg                                                                                                      <  1:1601828/150‑1 (MQ=255)
                                                      gTCTCTGAAATTCGCTTCTCCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTCATGTTCCGCTTGTTGGTAATTTAGTGACCAATCATTGGCGCTtg                                                                                                >  1:1428421/1‑150 (MQ=255)
                                                         tctGATATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTgg                                                                                             <  2:2126889/150‑1 (MQ=255)
                                                         tatGATATTCGCTTCACAGCACGAGAAAAATTGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTgg                                                                                             <  1:91676/148‑1 (MQ=255)
                                                               aTTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTaa                                                                                       >  2:1449061/1‑150 (MQ=255)
                                                                  cGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCg                                                                                    >  2:485271/1‑150 (MQ=255)
                                                                     ttCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCg                                                                                 <  2:1071006/150‑1 (MQ=255)
                                                                      tCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCgg                                                                                <  2:1395227/150‑1 (MQ=255)
                                                                      tCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCgg                                                                                <  1:1134155/150‑1 (MQ=255)
                                                                      tCACCGCACGAGAAAAAATGACCGGGCTGGAAAGGTATGATGTCATAATCAAAAGAATGTGATTGTGTGAACATCCCCCCGTGCTAGCTGATGTTCCGATTGTTGCTAATTTAGTGACCAATCAATGCCGCTTGTGGAATTAAGctgcgg                                                                                >  2:1055321/1‑145 (MQ=255)
                                                                      gcACCGCACGAGAAAAAAAGACCGGGCTTGGAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGcggcgg                                                                                <  2:2171729/149‑1 (MQ=255)
                                                                          cGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACa                                                                            >  1:282617/1‑150 (MQ=255)
                                                                             aCGATAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAAtt                                                                         <  1:1066275/150‑1 (MQ=255)
                                                                              cGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtc                                                                        <  1:28234/150‑1 (MQ=255)
                                                                                agaAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtcct                                                                      >  2:2155051/1‑150 (MQ=255)
                                                                                agaAAAAATGACCGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtcctc                                                                     >  2:1547750/1‑150 (MQ=255)
                                                                                 gaAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtcctc                                                                     <  1:1604043/150‑1 (MQ=255)
                                                                                   aaaaaTGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCGCCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCg                                                                   >  1:1517803/1‑150 (MQ=255)
                                                                                   aaaaaTGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAATCGTCGGTACAATTCCTCCg                                                                   <  1:145201/150‑1 (MQ=255)
                                                                                   aaaaaTGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCg                                                                   <  1:1411112/150‑1 (MQ=255)
                                                                                   aaaaaTGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCg                                                                   <  1:829988/150‑1 (MQ=255)
                                                                                    aaaaTGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGCCCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCgg                                                                  >  1:219122/1‑150 (MQ=255)
                                                                                         gACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCg                                                             >  1:1418024/1‑150 (MQ=255)
                                                                                            cGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTTAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGc                                                          <  2:443923/150‑1 (MQ=255)
                                                                                             gggCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCt                                                         >  2:1999857/1‑150 (MQ=255)
                                                                                             gggCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCt                                                         >  2:650373/1‑150 (MQ=255)
                                                                                               gCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGATACAATTCCTCCGGCACCGGGCttt                                                       >  2:2186515/1‑150 (MQ=255)
                                                                                                 ttGAAAGTTATGATGTCAAAATCAAAATACTGGGATTTTGTGAACATCACCCCGGGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGc                                                     <  1:1183255/150‑1 (MQ=255)
                                                                                                 ttGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGc                                                     <  1:505671/150‑1 (MQ=255)
                                                                                                 ttGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGc                                                     >  2:1063952/1‑150 (MQ=255)
                                                                                                  tGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGcc                                                    <  1:844913/150‑1 (MQ=255)
                                                                                                  tGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGcc                                                    <  2:2051165/150‑1 (MQ=255)
                                                                                                   gAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGGCCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCa                                                   <  2:768091/150‑1 (MQ=255)
                                                                                                   gAAAGTTATGAAGTCAAAATCAAAATAATGTGATTTTGTGAACATCACACCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTTACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCa                                                   >  1:1453919/1‑150 (MQ=255)
                                                                                                         tATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTcc                                             >  2:207237/1‑150 (MQ=255)
                                                                                                           tgatgTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCg                                           <  1:548478/150‑1 (MQ=255)
                                                                                                               gTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATa                                       <  1:1438939/150‑1 (MQ=255)
                                                                                                                tcaaaatcaaaaTAATGTGATTTTGTGAACATCACCCCGTGCGAGCTGATGTTCCGCTTGTTGCTAAATTAGTTAACAATCATTGGCGCTTGTGGAATTAAGCGTAGGTACAATTCCTCCGGCACCGGGCTTTGCCATAATCCCGCATAc                                      >  1:1094202/1‑150 (MQ=255)
                                                                                                                 caaaatcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACa                                     <  2:1302719/150‑1 (MQ=255)
                                                                                                                 caaaatcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACa                                     >  1:968573/1‑150 (MQ=255)
                                                                                                                        aaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCGTGTTGCTACTTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGAGTCTGTACAATTACTACGACACAGGCCTTTGCCATACTCCCGACTACATTGCTTa                              >  2:356861/1‑150 (MQ=255)
                                                                                                                           ataatGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAAcc                           <  1:2193513/150‑1 (MQ=255)
                                                                                                                           ataatGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAAcc                           <  1:1444778/150‑1 (MQ=255)
                                                                                                                              atGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAAt                        >  1:459250/1‑150 (MQ=255)
                                                                                                                               tgtgATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATc                       >  1:825300/1‑150 (MQ=255)
                                                                                                                                 tgATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCAc                     >  2:1043595/1‑150 (MQ=255)
                                                                                                                                 tgATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCAc                     >  2:1742471/1‑150 (MQ=255)
                                                                                                                                  gATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTTTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTATCCAATCAcc                    >  1:1472837/1‑150 (MQ=255)
                                                                                                                                  gATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCAcc                    <  2:1631991/150‑1 (MQ=255)
                                                                                                                                  gATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCAcc                    >  1:2097164/1‑150 (MQ=255)
                                                                                                                                    ttttGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGc                  <  2:581111/150‑1 (MQ=255)
                                                                                                                                       tgtgAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCaaa               >  1:987018/1‑150 (MQ=255)
                                                                                                                                        gtgAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATTCTCCCGCATACATTGCGTAACCAATCACCGCaatc              <  1:83731/150‑3 (MQ=255)
                                                                                                                                        gtgAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCGTTGCCATACTCCCGCATACATTGCGTAACCAATCCCCGCAAAc              >  2:1963552/1‑150 (MQ=255)
                                                                                                                                        cggAACATCACCTCGCACGGGGTGATTTTCCGCTTGTTGCTAATTTAGTGACCAATCATGGGCGCGTGGGGAATTAAGCGTCGGTACAATTGCTCCGGCACCGGTCGTTGCCATACTCCCGCATGCATTGCGTTACCAAGCACCGCaacc              >  2:872841/3‑148 (MQ=255)
                                                                                                                                             cATCACCTCGCACGGGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATaaa         >  1:948284/1‑149 (MQ=255)
                                                                                                                                             cATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAAt         <  1:1986809/150‑1 (MQ=255)
                                                                                                                                              aTCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATc        >  1:2202020/1‑150 (MQ=255)
                                                                                                                                                 aCCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATcccc     >  1:72076/1‑150 (MQ=255)
                                                                                                                                                  ccccGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCa    <  1:576486/150‑1 (MQ=255)
                                                                                                                                                    ccGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTTCCATACTCCCGCATACATTGCGTAACCTATCACCGCAAACATAATCACCAgc  >  1:389367/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCGACTTTATCGCGCAAAAAATCGCCGAGCTGAACAAGAATACAGGGAAAGAAGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGA  >  NC_000913/1052109‑1052406

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: