Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 136,535 G→C 6.9% intergenic (‑71/+35) speE ← / ← yacC spermidine synthase/putative lipoprotein YacC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913136,5350GC 6.9% 526.3 / 20.4 158intergenic (‑71/+35)speE/yacCspermidine synthase/putative lipoprotein YacC
Reads supporting (aligned to +/- strand):  ref base G (87/59);  new base C (4/7);  total (91/67)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 158 total aligned reads displayed.

ATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAGGGCT  >  NC_000913/136387‑136681
                                                                                                                                                    |                                                                                                                                                  
 tACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCa                                                                                                                                                  >  1:364956/1‑150 (MQ=255)
 tACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCa                                                                                                                                                  >  1:2114954/1‑150 (MQ=255)
 tACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCa                                                                                                                                                  >  1:1751977/1‑150 (MQ=255)
   cAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACg                                                                                                                                                >  1:2065345/1‑150 (MQ=255)
        cGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAGGGCGCAACATAATAGCTAACATTGACCGTGCATCCGCGGTCAat                                                                                                                                          <  1:925865/150‑1 (MQ=255)
          ttATCTACGGCAAAGTACTGCGCAAACTGGTCGTGTTGTGTTTCATGCCAATGTTTTTTTTCGGCCATTGGTTGATACCTCCTTTGTTAAGAACCGTAAAAAATGGGCGCAACTTAATAGCTAACATTGACCGCGGATGCACGGTCAata                                                                                                                                         <  2:92181/150‑1 (MQ=255)
            aTCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAATATTGACCGTGCATCCGCGGTCAATGtt                                                                                                                                       <  2:1817935/150‑1 (MQ=255)
             tCTACCGCAAAGTACTGCCCAAACTGTTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATAttt                                                                                                                                      <  2:2389622/150‑1 (MQ=255)
             tCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGTGCATCCGCGGTCAATGtta                                                                                                                                      <  2:799506/150‑2 (MQ=255)
              cTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGTGCATCCGCGGTCAATGttag                                                                                                                                     <  1:1020116/150‑3 (MQ=255)
              cTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGTGCATCCGCGGTCAATGttag                                                                                                                                     <  2:1799428/150‑3 (MQ=255)
              cTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGTGCATCCGCGGTCAATGttag                                                                                                                                     <  2:2057506/150‑3 (MQ=255)
              cTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATAtttg                                                                                                                                     >  1:521539/1‑149 (MQ=255)
              cTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTc                                                                                                                                     <  1:1608147/150‑1 (MQ=255)
              cTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTc                                                                                                                                     >  2:1100255/1‑150 (MQ=255)
                aCCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTa                                                                                                                                   <  2:1536120/150‑1 (MQ=255)
                aCCGCAAACTACTGCCCAAACTGCTCGTGTACCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGTGCATCCGCGGTCAATGttagct                                                                                                                                   <  1:2485568/150‑5 (MQ=255)
                 ccGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTAc                                                                                                                                  <  1:2403467/150‑1 (MQ=255)
                 ccGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTAc                                                                                                                                  >  1:1339662/1‑150 (MQ=255)
                 ccGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTAc                                                                                                                                  >  1:2032759/1‑150 (MQ=255)
                   gCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCCCTGTTTTTTTTCGGCCATGCGTTGATACCTCCTTTGTTAACACCCGGAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACaa                                                                                                                                >  2:2174884/1‑150 (MQ=255)
                   gCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACaa                                                                                                                                >  1:349093/1‑150 (MQ=255)
                      aaGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACaaaaa                                                                                                                             >  2:2434918/1‑150 (MQ=255)
                          aCTGACCAAACTGGTCGTTTAGCGTTTCATTCCTGTGTTTTTTTTCTGCGATGGGTTGATACGTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTa                                                                                                                         <  2:871200/150‑1 (MQ=255)
                              cccAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAACAGGTATCAg                                                                                                                     >  2:2179480/1‑150 (MQ=255)
                              cccAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAg                                                                                                                     >  1:829588/1‑150 (MQ=255)
                               ccAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAgg                                                                                                                    <  1:491090/150‑1 (MQ=255)
                               ccAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAgg                                                                                                                    <  1:267948/150‑1 (MQ=255)
                               ccAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAgg                                                                                                                    <  2:2255167/150‑1 (MQ=255)
                                cAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAggg                                                                                                                   >  2:2299833/1‑150 (MQ=255)
                                      ggTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGattatt                                                                                                             >  1:1280260/1‑150 (MQ=255)
                                       gTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTAttt                                                                                                            <  1:1432469/150‑1 (MQ=255)
                                          gtgtAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGAc                                                                                                         <  2:2065345/150‑1 (MQ=255)
                                           tgtAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACg                                                                                                        >  1:509520/1‑150 (MQ=255)
                                               gCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAgg                                                                                                    <  1:1018541/150‑1 (MQ=255)
                                               gCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAgg                                                                                                    >  2:362205/1‑150 (MQ=255)
                                               gCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAgg                                                                                                    <  2:973336/150‑1 (MQ=255)
                                                  tttCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGAGTATTTGACGTAGGCaa                                                                                                 >  2:1187166/1‑150 (MQ=255)
                                                     cATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCa                                                                                              <  1:1444443/150‑1 (MQ=255)
                                                       tGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAgg                                                                                            >  2:529695/1‑150 (MQ=255)
                                                       tGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGc                                                                                           >  2:1383688/1‑150 (MQ=255)
                                                        gCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGc                                                                                           >  2:502224/1‑150 (MQ=255)
                                                         ccACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCt                                                                                          <  2:364956/150‑1 (MQ=255)
                                                          cACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCtt                                                                                         <  2:1280260/150‑1 (MQ=255)
                                                            cTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTaa                                                                                       <  2:509520/150‑1 (MQ=255)
                                                             tGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAg                                                                                      <  2:614160/150‑1 (MQ=255)
                                                             tGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAg                                                                                      >  1:28893/1‑150 (MQ=255)
                                                             tGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAg                                                                                      <  2:2114954/150‑1 (MQ=255)
                                                                tttttttCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGaa                                                                                   >  2:2413005/1‑150 (MQ=255)
                                                                       cGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCggg                                                                            <  1:1111503/150‑1 (MQ=255)
                                                                        ggCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGc                                                                           <  1:168523/150‑1 (MQ=255)
                                                                         gCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGcc                                                                          >  2:1819470/1‑150 (MQ=255)
                                                                            aTGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAgg                                                                       >  1:1129328/1‑150 (MQ=255)
                                                                               ggTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGTTTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCt                                                                    <  2:327611/150‑1 (MQ=255)
                                                                                        ctcctTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACAGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGAAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAtttt                                                           >  1:1716768/1‑150 (MQ=255)
                                                                                        ctaattttttAACACCCGTAAAAAAAGTTCGCAACATAATATATAAAATTCACCGCCAATGCACTCTCAATATTTCTACAAAAAGCTAACAGGGATTTTTTTAAATAGGCAAGAACTATTAAGGAATCGCCGCCCATGGATTTGCAtttt                                                           >  2:74433/5‑150 (MQ=255)
                                                                                          cctTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAtttttt                                                         <  1:1662314/150‑1 (MQ=255)
                                                                                          cctTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCAtttttt                                                         <  1:436048/150‑1 (MQ=255)
                                                                                            tttGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAg                                                       <  2:348615/150‑1 (MQ=255)
                                                                                            tttGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAg                                                       <  1:2304159/150‑1 (MQ=255)
                                                                                               gTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAGTCGCGGGCCAGGGCTTTGCATTTTTTAGCgc                                                    >  2:1187268/1‑149 (MQ=255)
                                                                                                ttAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGAc                                                   >  2:2315374/1‑150 (MQ=255)
                                                                                                ttAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGAc                                                   <  1:2367041/150‑1 (MQ=255)
                                                                                                 tAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACg                                                  >  1:985285/1‑150 (MQ=255)
                                                                                                       ccGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAATCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatg                                            >  1:2480905/1‑150 (MQ=255)
                                                                                                       ccGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatg                                            >  2:1639723/1‑150 (MQ=255)
                                                                                                       ccGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatg                                            >  2:1929071/1‑150 (MQ=255)
                                                                                                        cGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTCTCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgc                                           <  1:1915003/150‑1 (MQ=255)
                                                                                                        cGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgc                                           <  2:2042043/150‑1 (MQ=255)
                                                                                                         gTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgcc                                          <  2:663236/150‑1 (MQ=255)
                                                                                                          tAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgcca                                         >  2:1421248/1‑150 (MQ=255)
                                                                                                           aaaaaaaGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgccaa                                        <  1:787885/150‑1 (MQ=255)
                                                                                                            aaaaaaGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgccaat                                       >  1:2148927/1‑150 (MQ=255)
                                                                                                               aaaGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGaatgccaatgcc                                    <  1:2065207/150‑1 (MQ=255)
                                                                                                                    gcgcAACATAATAGCTAACATTGACCGCGGATGCACGGTGAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTga                               >  2:2202092/1‑150 (MQ=255)
                                                                                                                     cgcAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTgag                              <  1:2042716/150‑1 (MQ=255)
                                                                                                                       cAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGAt                            >  2:2216700/1‑150 (MQ=255)
                                                                                                                           ataataGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATcgcg                        >  2:1766031/1‑150 (MQ=255)
                                                                                                                            taataGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAATCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCgg                       <  2:2480905/150‑1 (MQ=255)
                                                                                                                            taataGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCgg                       <  1:502224/150‑1 (MQ=255)
                                                                                                                              ataGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTa                     <  1:687210/150‑1 (MQ=255)
                                                                                                                              ataGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGATTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTa                     >  2:827235/1‑150 (MQ=255)
                                                                                                                                aGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGc                   >  1:57094/1‑150 (MQ=255)
                                                                                                                                 gCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCt                  >  1:329140/1‑150 (MQ=255)
                                                                                                                                  cTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTg                 >  1:1145626/1‑150 (MQ=255)
                                                                                                                                  cTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTg                 >  1:1100805/1‑150 (MQ=255)
                                                                                                                                  cTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGCATTATTTGACGTAGGCAAGCAAGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTg                 <  1:979986/150‑1 (MQ=255)
                                                                                                                                    aaCATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTc               >  2:1987343/1‑150 (MQ=255)
                                                                                                                                      taTTGACCGTGCATCCGCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCtt             >  2:1489860/2‑150 (MQ=255)
                                                                                                                                        ttGACCGTGCATCCGCGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCa           >  2:2159741/1‑150 (MQ=255)
                                                                                                                                         tGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCGTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCAc          <  1:2234808/150‑1 (MQ=255)
                                                                                                                                          gACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCAcc         >  1:1723289/1‑150 (MQ=255)
                                                                                                                                          gACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCAcc         <  2:529231/150‑1 (MQ=255)
                                                                                                                                           aCCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCg        >  1:1314302/1‑150 (MQ=255)
                                                                                                                                           aCCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCg        <  2:39603/150‑1 (MQ=255)
                                                                                                                                           aCCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCg        >  2:476256/1‑150 (MQ=255)
                                                                                                                                           aCCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCg        <  1:1819470/150‑1 (MQ=255)
                                                                                                                                             cgcgGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAg      >  2:2061845/1‑150 (MQ=255)
                                                                                                                                              gcgGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAgg     >  2:1268851/1‑150 (MQ=255)
                                                                                                                                              gcgGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAgg     <  1:2315374/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGACGTAGGCAAGCAGGCTTAAGGAATCGCGGGCCAGGGCTTTGCATTTTTTAGCGACGGGAATGCCAATGCCGCTGAGATCGCGGTAGCTGTCTTCACCGAGGGCT  >  NC_000913/136387‑136681

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