Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8800TC50.0% 91.4 / 37.5 84K271E (AAA→GAA) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (35/7);  new base C (0/42);  total (35/49)
Fisher's exact test for biased strand distribution p-value = 1.02e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357732‑358028
                                                                                                                                                           |                                                                                                                                                          
tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                         >  2:849935/1‑151 (MQ=255)
 cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccg                                                                                                                                                        >  2:1078544/1‑151 (MQ=255)
 cGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcc                                                                                                                                                         <  2:957394/150‑1 (MQ=255)
    aGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                                     >  2:938952/1‑151 (MQ=255)
     ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCcgccgc                                                                                                                                                    >  1:336945/1‑151 (MQ=255)
         tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                                 >  2:1673857/1‑150 (MQ=255)
          gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                                <  1:485522/150‑1 (MQ=255)
          gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                                >  2:1141075/1‑150 (MQ=255)
           aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                               >  2:1496254/1‑150 (MQ=255)
             actcgcaGCGCACCAGGCAATTCCAGCCGCAGACCTGTTTCTTCGGCTACCGTGATTCATTTCCGCCAACACCCGCGCATTTATCCAACGCCATGTGCAAACATGCATTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACaaaa                                                                                                                                             <  2:444192/144‑1 (MQ=255)
               gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑ccg                                                                                                                                        <  2:686945/151‑1 (MQ=255)
                cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑ccgc                                                                                                                                       <  2:271784/151‑1 (MQ=255)
                      gcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTtct                                                                                                                                 >  1:300182/1‑151 (MQ=255)
                       cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTtctc                                                                                                                                >  2:1175717/1‑151 (MQ=255)
                        aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTtctct                                                                                                                               >  1:1761495/1‑151 (MQ=255)
                          cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTcc                                                                                                                             >  1:501054/1‑151 (MQ=255)
                              cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCagt                                                                                                                          >  1:221927/1‑150 (MQ=255)
                              cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCagt                                                                                                                          >  2:364711/1‑150 (MQ=255)
                              cAATTCCAGCCGCAGACCTGTGTCAGCGGCTACAGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAC‑‑CCGCCGAtctct                                                                                                                               >  2:886532/1‑145 (MQ=255)
                               aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCagt                                                                                                                          >  2:282525/1‑149 (MQ=255)
                                   ccgcccGCGGCCCTGTGTCAGCGGCTCCCGTGTTTCATTCCCGCCACCACCCGCGCATTTACCCACCGCCTTGTGAAAAAATGCCTCCGCGGCGGCTTTCTGCCAGTTTTTTTCCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTgg                                                                                                                     <  1:1360504/146‑1 (MQ=255)
                                   ccgcccGCAGCCCTGTGTCAGCGCCTCCCGTATTTCATTCCCCCCACCACCCGCGCTTTTACCCACCCCCATGTGAAAAAATCCTTCCGCGGGGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTgg                                                                                                                     <  1:621951/146‑1 (MQ=255)
                                   ccAACCGCAACCCTGTCTCGCCGCCTACGGTGATTCATTCCCTCCAACAGCCGCGCTTTTCCCCAACCCTATGTGCAAAAATCCCTCCGGGGGGGCTTCCTCCCAGTTTTTTTCCCCCCCCAAAAAAG‑‑CCGCCGTTCTCTCCAGTAGTGGc                                                                                                                    <  1:1208048/151‑1 (MQ=255)
                                       ccGCAGCCCTGTGTCAGGGGCTCCCGTGATTCTTTCCCCCCAACACCCGCTCTTTTACCCACCGCCTTGTGAAAAACTGCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                 <  1:1184444/150‑1 (MQ=255)
                                       ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                 >  1:1618270/1‑150 (MQ=255)
                                         gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                >  1:355433/1‑149 (MQ=255)
                                          cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                >  1:1358088/1‑148 (MQ=255)
                                           accccTGTGTCAGCGTCTCCCGTTTTTCTTTCCCCCCAACACCCGCGCTTTTCCCCACCCCCATGTAAAAAAATGCTTTCGGGGCGGTTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAag                                                                                                            <  1:391625/148‑1 (MQ=255)
                                            gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCggggg                                                                                                               >  2:1459979/1‑146 (MQ=255)
                                             aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                >  1:690899/1‑145 (MQ=255)
                                              ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCgggggggg                                                                                                            >  2:243700/1‑144 (MQ=255)
                                                  ttccAGCGCCTCCCGTATTTCATTCCCCCAAACACCCGCGCATTTTCCCAACGCCATTTGCAAAATTGCCTTCGCGGCGGCTGTCTCCCAGCTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                                      <  1:1132613/147‑1 (MQ=255)
                                                  tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAGGAGaaaa                                                                                                      >  1:258826/1‑149 (MQ=255)
                                                   gccggggcctcccGTTATTCTTTCCCCCAAACACCCGCCCTTTTACCCACCCCCAGGTGAAAAAAGGCCTCCGGGGCGGCTGTCGCCCAGTTTTTTCCCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                                         <  2:87153/135‑1 (MQ=255)
                                                     cggcgcCTCCCGTGTTTCTTTCCCCCAAACACCCCCGCTTTTACCCACCGCCGTGTAAAAAAATCCCTCCGCGGCGGCTTTCCGCCATTTTTTTTCCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCaa                                                                                                             <  2:1342987/138‑1 (MQ=255)
                                                     cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGgaaagaaaaa                                                                                                     >  1:150684/1‑146 (MQ=255)
                                                         ggCTACCGTGATTCTTTCCCCCCAACACCCCCGCTTTTACCCACCCCCATGGCAAAAATTCCCTCCGCGGCGGTTGTCGCCCAGCTTTTTTCCCCCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                             <  1:1674945/151‑1 (MQ=255)
                                                         gcttcccGTGTTTCTTTCCCCCAAACACCCGCGCTTTTACCCACCGCCATGTCAAAAATTGCTTTCGGGGGGGTTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTa                                                                                              <  1:121229/145‑1 (MQ=255)
                                                         gcctgccGGGATTTTTTTCCCCCAACACCCGCGATTTACTCCACCGCTATGTGAAAAAGTCCCTTCGCGGCGTCTGTTGGCAGGTTTTTTTCCCGCCACAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGc                                                                                                   <  1:1726483/141‑1 (MQ=255)
                                                         gcctcccGTGTTTCATTCCCCCAACCCCCCGCCCTTTTACCCAACCCCTTGTAAAAAAATCCCTGCGGGGCGGCTGTCGGCAACTTTTTTTCCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa                                                                                                       <  2:1551206/137‑1 (MQ=255)
                                                          ggcaccgttaattctttCCCCCAACCACCCGCCCTTTTCCCCACCCCCATGTGAAAAAAGCCCTCCGGGGCGGCTTCCGGCCGCTTTTTTTCCCCCCACAAAAAAG‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCttt                                                                                               <  1:145877/135‑1 (MQ=255)
                                                              ccgtttttCTTTCCCCCCAACACCCGCGCTTTTATCCACCCCATTGTGCAAAAATGCCTTCGCGGCGGTTTTCTGCCAGTTTTTTTCCCGCCCAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                                                        <  1:396902/145‑1 (MQ=255)
                                                              ccgtttttCATTCCCCCAAACACCCGCCATTTCCCCCACCCCATGGTAAAAAATTGCCTGCGGGGCGGTTTTCTCCCAGTTTTTTTCCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                                                                          <  2:1490047/143‑1 (MQ=255)
                                                              ccGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCACCCCCATGTGAAAAAATGCCTGCGGGGGGGTTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                                                        <  1:895558/151‑1 (MQ=255)
                                                              ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAGGCGaaaaaaactttaaaacc                                                                                         >  1:1112961/1‑137 (MQ=25)
                                                              ccGTGATTCATTCCCCCCACCACCCGCGCTTTTACCCACCGCCATGTGCAAAAATGCTTTCGGGGCGGCTGTCTCCCAGCTTTTTTCCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                                                        <  1:1570986/151‑1 (MQ=255)
                                                                gTGATTCTTTTCCCCAAACACCCGCGCTTTTACCCAACGCCATTTGAAAAAATCCCTTCGCGGCGGTTTTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCg                                                                                        <  1:1500348/149‑1 (MQ=255)
                                                                gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGGGAGaaaaaaattttaaccc                                                                                          >  1:833361/1‑135 (MQ=25)
                                                                   tttCATTCCCCCAAACACCCGCTCTTTTACCCACCCCCATGTAAAAAATTCCTTGCGCGGGGGCTGTCCGCCAGTTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa                                                                                     <  2:1278433/148‑1 (MQ=255)
                                                                    ttCTTTCCCCCCAACACCCGCCCTTTTACCCACCGCCATGTGAAAAAATCCCTCCGGGGCGGCTGTCTGCAAGTTTTTTTCCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCa                                                                                     <  1:1754055/148‑1 (MQ=255)
                                                                    ttCATTCCCCCAAACACCCGCGCTTTTACCCACCGCCTTGTAAAAAAATCCCTGCGGGGCGGCTTTCTCCCAGTTTTTTTCCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAt                                                                                    <  1:814042/149‑1 (MQ=255)
                                                                        tttCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCGCGAGGAAAAGGttttaatcccccccggt                                                                                 >  2:92223/1‑134 (MQ=25)
                                                                        tccccGCCAACACCCTCCAATTTAGCCAGCGCCATGTCAAAAACTGCCTCCGCGGTGGTTGTCTCCACGTTATTTCCCCGCGCAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttgt                                                                               <  1:1213426/147‑1 (MQ=255)
                                                                          ccccccAACACCCGCGCTTTTACCCAACGCCTTTTGAAAAAAGCCCTCCGCGGCGGCGTTCTCCCTGTTTTTTTCCCGCCACAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                           <  1:938787/147‑1 (MQ=255)
                                                                           ccccaaacaccCGCGTTTTTACCCACCGCCAGGTCAAAAAATCCCTCCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCAAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                          <  2:609523/146‑1 (MQ=255)
                                                                           ccccaaacaccCGCGCTTTTACCCACCGCCATGTAAAAAAAGGCCTTCGCGGGGGTTTTCTGCCGGTTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                          <  1:197563/146‑1 (MQ=255)
                                                                                      ccGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCCGAAGAAAAAACATTTAAACCCGCAtgttgttgtg                                                                           >  2:1530859/1‑137 (MQ=25)
                                                                                       cgcataTTTACCCACCGCCTTTGAAAAAATGCCCTTGGGGGGGGTTGTCGCCCAGTTTTTTTCCCGCCTAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                               <  2:840735/146‑1 (MQ=255)
                                                                                         cgcATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAAAGAAAAAAATTTAAAACCGCCCTGTTTTGTGGGGGTAGctgctg                                                             >  2:1356024/1‑150 (MQ=25)
                                                                                             tttACCCACCGCCATGTAAAAAATTGCTTCCGGGGGGGCTTTCCCCAAGTTTTTTTCCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                          <  1:1385519/149‑1 (MQ=255)
                                                                                                    caacccccttgccaaaaagTCCCTCCGGGGCGCTTGTCCCCCAGTTTTTTTCCCCCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                        <  2:833322/132‑1 (MQ=255)
                                                                                                     accGCCATGGGAAAAAAGGCCTTCGGGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                   <  1:849974/147‑1 (MQ=255)
                                                                                                         ccATGTGAAAAAGTCCCTCCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaa                                               <  1:1175796/149‑1 (MQ=255)
                                                                                                         caaTGTGAAAAAACGCTTCCGGGGGGGCTTTCCCCCAGTTTTTTTCCCCCCGAAAAAAAC‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                              <  2:1244408/148‑1 (MQ=255)
                                                                                                          cATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAATTTAAAACCCCCCTGTTGTGTGGGGATCGCTGCTGGGTACAATGTgg                                                 >  1:155441/1‑146 (MQ=37)
                                                                                                           aTGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAACTTTTAAACCCTCCTGTTGTGTGGGGAACGCTGCTGGGGACAAAGTGGGAAgg                                            >  1:151766/1‑150 (MQ=255)
                                                                                                            tgtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAAAAAAAAAATTTAAACCGCCCTGCTGTTTGGGAAACGCCGCTGGGAACAAtgatg                                                 >  2:1421831/1‑141 (MQ=255)
                                                                                                            tgtgAAAAAATGCCTCCGGGGCGGTTTCCTGCCAGTTATTTTCCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                           <  2:1037665/150‑1 (MQ=255)
                                                                                                               taaaaaaaTCCCTCCGGGGGGGCTGTCGCCAAGTTTTTTTCCCGCCAAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                              <  1:774948/142‑1 (MQ=255)
                                                                                                               gtaaaaaaTCGCCTTCGGGGCGGCTGTCTCCCATTTTTTTTCCCCCCAAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                         <  2:1626142/148‑1 (MQ=255)
                                                                                                                 aaaaatTGCTTCCGGGGGGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                      <  1:1460675/150‑1 (MQ=255)
                                                                                                                  aaaaaTGCCTCCGCGGCGGTTTTCTCCCAGTTTTTTTCCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                         <  2:524827/147‑1 (MQ=255)
                                                                                                                  aaaaaGCCCTGCGCGGCGGCTGTCGCCCAGTTTTTTTCCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                          <  2:1367030/146‑1 (MQ=255)
                                                                                                                      aaaagtccttccggggcggcTGTCGCCCAGTTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGGAATGGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                       <  1:1674015/135‑1 (MQ=255)
                                                                                                                              ctcggggggggtggTCTCCCAGTTTTTTTCCCCCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGa                                <  2:336928/135‑1 (MQ=255)
                                                                                                                                  gggggggcTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                              <  1:1360944/141‑1 (MQ=255)
                                                                                                                                    ggcggcTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAATTTAAACCCCCCCTGTTGTGTGGGGAAAGCCGCTTGGGACAAAGTGGGAAGGGGAGTGGGGGGGATCAACTcc                       >  1:925086/1‑151 (MQ=255)
                                                                                                                                    ggcggcTGTCCCCCATCTTTTTTCCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc                       <  1:174616/151‑1 (MQ=255)
                                                                                                                                           gTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAAAGAAAAAACATTTAAGACCCCATGTTGTGGGGGAACAGCTGCCGGGAACAAATTGGGAAGGGGAGGGGGGGGGAACAGCCCCCGaaac                 >  1:813919/1‑148 (MQ=255)
                                                                                                                                            tCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGAAAAAAATTTTAAACCCCCCTTTTTTGTGGGGAACGCTGCTGGGTAAAAAGGGGGAAGGGGAGGGGGGGGGGACAGCTCCAGaa                   >  1:1226859/1‑147 (MQ=255)
                                                                                                                                               gCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAACTTTTAAACCCTCAATTTGTGTGGCAAACCCCGCCGGGTACAATTTGGGGAAGGGAGGGGGGGGGATCAGCACCCGAAAcgcccc            >  1:1765955/1‑148 (MQ=255)
                                                                                                                                                ccattttttttcCCGCCACAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa           <  2:1504798/144‑1 (MQ=255)
                                                                                                                                                ccattttttttcCCCCCACAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTa           <  2:1761331/144‑1 (MQ=255)
                                                                                                                                                ccAGCTTTTTTCCCCCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc               <  1:243709/147‑1 (MQ=255)
                                                                                                                                                ccAGCTATTTTCCCCCCCAAAAAAAC‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                     <  1:1011785/141‑1 (MQ=255)
                                                                                                                                                 ccgCTTTTTTCCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa          <  1:481697/149‑1 (MQ=255)
                                                                                                                                                     ttttttcCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt         <  1:25767/146‑1 (MQ=255)
                                                                                                                                                     ttttttcCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg       <  1:504264/148‑1 (MQ=255)
                                                                                                                                                     ttttttcCCGCCACAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAtt        <  1:26326/147‑1 (MQ=255)
                                                                                                                                                      ttttttCCCCCCCAAAAAAAACCCGCCGTTCCCCCAGGTATGGGCGGGGAAAGAGAAATAGTTTTAAGCCCTTAAGTTGTGGGGGCATCCCCTGCTGGTACAATGTGGGAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg       <  2:211673/149‑1 (MQ=255)
                                                                                                                                                       ttttcCCGCCGCAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAg     <  1:830430/150‑1 (MQ=255)
                                                                                                                                                        tttcCCGCCGAAAAAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt    <  1:775460/150‑1 (MQ=255)
                                                                                                                                                        tttcCCCCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGt    <  2:1350731/150‑1 (MQ=255)
                                                                                                                                                         tttCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAACTTTTAGACCGCCCTGTTTTGTGGGGAACGCTGCTGGGAAAAAAATGGGAAGGGGAGGGCGGGGGATCAACTCCCCAACCGCGCCAATTTAGttt  >  1:693299/1‑150 (MQ=255)
                                                                                                                                                           |                                                                                                                                                          
TCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAA‑‑CCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTC  >  NC_000913/357732‑358028

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