Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA51.5% 38.8 / 103.3 69L267F (TTG→TTTcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (27/6);  new base A (0/35);  total (28/41)
Fisher's exact test for biased strand distribution p-value = 1.63e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.85e-09

GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCA  >  NC_000913/357742‑358034
                                                                                                                                                          |                                                                                                                                                    
gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACAGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                                     >  2:694653/1‑150 (MQ=255)
   aaGCGTAGCGCATCATGAAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                                 <  2:1141729/151‑1 (MQ=255)
    aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                                <  1:41165/151‑1 (MQ=255)
    aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                                >  2:1009346/1‑151 (MQ=255)
      cgcaGCCCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTCCAAAAATGCCTTCGCGGCGGCTGCCTTCCAGCAATTTTTCCTCCCCACCAAAAccg                                                                                                                                               >  2:938793/4‑150 (MQ=17)
       gTAGCGCACCAGGCAATTCCAGCCGCTGACCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTCTCCAACGCCATGTGCAAAAATGCCTTCGCGTCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                             <  2:1054173/151‑1 (MQ=255)
            gcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTCTCCGCAACACAAcggcaactc                                                                                                                                         >  2:274509/1‑145 (MQ=14)
             cATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTTCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTATTTTTTCGCCGGAATACAAccagaacg                                                                                                                                          >  2:41159/1‑142 (MQ=11)
              aTCAGGCAATTCCAGCCGCAGACCTGTGTTAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctct                                                                                                                                      >  1:632759/1‑151 (MQ=255)
              aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtctc                                                                                                                                       >  1:1075833/1‑150 (MQ=255)
                cAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                                    >  1:100978/1‑151 (MQ=255)
                 aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTcc                                                                                                                                    >  1:305618/1‑150 (MQ=255)
                     aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                                 >  1:478952/1‑149 (MQ=255)
                         ccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                           >  1:718399/1‑151 (MQ=255)
                          cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCg                                                                                                                         >  1:335326/1‑151 (MQ=255)
                             ccGCAGACTTGTGTCAGCGGCTCCCGTGTTTCTTTCCCCCCAACACCCGCGCTTTTACCCAACGCCATTTGCAAAAATCCCTTCGGGGCGGCTGTCTCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCggg                                                                                                                       <  1:140607/150‑1 (MQ=255)
                               gCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                      >  1:965456/1‑149 (MQ=255)
                                   cccTGTGTCAGCGGCTCCCGTGATTCTTTCCCCCCACCACCCGCGCTTTTACCCACCGCCATGTGAAAAAATGCCTCCGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAaga                                                                                                                 <  1:804919/149‑1 (MQ=255)
                                     gtgtgTCAGCGCTTCCGGGGATTCTTTCCCCCCAACCCCCAAGCATTTACCAAACCCTATGACCAAAAATCCTTCGGCGGCGTCTGTTTGCCGTTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                      <  2:349882/142‑1 (MQ=255)
                                      tgtgTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                                      >  1:598246/1‑142 (MQ=255)
                                       gtgtCGGCGGCTACCGTGATTCTTTCCCGCAAACACCCGCGCATTTACCCAACGCCATGTGAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGaaa                                                                                                           <  1:680265/150‑1 (MQ=255)
                                        tgccggCGGCTCCCGTGATTCATTCCCCCCAACACCCGCGATTTTACCCACCGCCTTGTAAAAAAATCCCTTCGCGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                                          <  1:702814/145‑1 (MQ=255)
                                        tgccAGCGCCTCCTGTGTTTCATTCCCGCCAACACCCGCGCATTTACCCAACGCCATTTGAAAAAATCCCTGCGCGGGGGTTGTCTGCCGGTTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                                             <  1:399420/144‑1 (MQ=255)
                                            ggcgctaccccgttttCATTCCCCCAAACACCCGCCATTTCTCCCACGCCTTGAGAAAAAAAGCCTTTCGCGGGGGCTGCCTGAAAGTTTTTTTTCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                                        <  2:576318/134‑1 (MQ=255)
                                             gCGGCTCCCGTGATTCATTCCCCCCAACACCCGCTCATTTACCCACCCCCATGTAAAAAAATGCTTCCGGGGGGGCTTTCGCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAt                                                                                                          <  1:431867/145‑1 (MQ=255)
                                               gcctcctgtgtTTCATTTCCGCAAACACCCGCGCATTTACCCACCGCCATGTGAAAAAATGCCTCCGGGGCGTTTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTaa                                                                                                  <  1:381073/144‑1 (MQ=255)
                                                   aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGGAAAAAACTTTAAccc                                                                                               >  1:427304/1‑151 (MQ=25)
                                                    ccttgtttCTTTCCCCCCAACACCCGCGCATTCCCCCACCGCCATGTGCAAAATGCCCTCCGGGGCGGTGGTCTGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGc                                                                                                <  2:1252020/142‑1 (MQ=255)
                                                    ccgtatttCATTCCCCCAAAAACCCGCGCTTTTATCCACCGCCATGTAAAAAAATGCCTTCGCGGCGGCTGTCGGCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAg                                                                                                 <  1:263594/141‑1 (MQ=255)
                                                    ccgtattcaaTTCCCGCCAACACCCGCCCTTTTACCCACCGCCTGGTGAAAAAAGCCTTCCGCGGCGGCGCTCTGCCAGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGccc                                                                                              <  2:174607/141‑1 (MQ=255)
                                                    ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGaaaaaaactttaaacc                                                                                              >  1:461433/1‑137 (MQ=25)
                                                    ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGACCAGAaaaagaagattcaccgc                                                                                              >  1:1027173/1‑137 (MQ=25)
                                                         aTTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGGAAAAAAGCTTTa                                                                                                   >  2:422978/1‑141 (MQ=25)
                                                            cATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACTTTaaaaccccccagtt                                                                                      >  1:715327/1‑138 (MQ=25)
                                                              ttcccGCAAACACCCCCGCTTTTCCCCACCGCATTGTAAAAAAATGCCTCCGCGGGGGCTTTCTCCCATTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                                                      <  1:617610/144‑1 (MQ=255)
                                                                cccccaaacaccCGCGCTTTCCCCCACCGCCATGGGAAAAAAGCCCTCCGGGGCGGCTGTCGGCAATTGTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgtg                                                                                <  2:333896/145‑1 (MQ=255)
                                                                cccGCCAACACCCGCCTTTTTACCCACCCCCTTGTGAAAAAATGCCTTCGCGGCGGTTGTCTGCCAGTTTTTTTCCCGCGGAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtg                                                                                  <  2:457935/148‑1 (MQ=255)
                                                                 cccccAACAACCCCGTATTTACCCAACGCCATTTGCAAAAATGCCTGCGCGGGGGCTTTCTGCCGGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTtgtgt                                                                                 <  1:1254925/146‑1 (MQ=255)
                                                                       acaccccctcttttACCCACCCCCATGTAAAAAAAGCCCTCCGGGGGGTCTGTCTGCCTTTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                          <  1:1089607/139‑1 (MQ=255)
                                                                          acCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACTTTTATCCCCCCCTGTTTTGTGGGGATCg                                                                         >  1:521794/1‑148 (MQ=25)
                                                                          acCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAAAAACTTTAAAACCGTCATGTTGTGTGGGAATc                                                                          >  1:636108/1‑147 (MQ=255)
                                                                            cgcgCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGGGGAAAAAATTTATACCCGGCATCTTGTGTggggagtggtgc                                                                     >  2:146945/1‑140 (MQ=25)
                                                                               gcttttACCCACCGCCATGTGAAAAAATGCCTCCGGGGCGGTTTTCTGCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGctgctg                                                                   <  1:975438/146‑1 (MQ=255)
                                                                                cATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAACTTTAAAACCCCCCTGTTTTGTGGGAAACg                                                                         >  1:1011957/1‑142 (MQ=37)
                                                                                  tttcccaacccccgtgtgAAAAATTCCCTGCGGGGCGGTGTTCGGCCAGTTATTTGCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGt                                                                 <  2:1231295/134‑1 (MQ=255)
                                                                                  tttacccaccccattgtgAAAAAATGCCTTCGGGGCGTCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                                <  1:880875/135‑1 (MQ=255)
                                                                                  tttATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGGGGAAAAACATTAAGCCCATCCTGGTTGTGGGGGATCCCTGCTgg                                                                  >  2:241619/1‑147 (MQ=25)
                                                                                  tttACCCACCGCCTTGTAAAAAATTCCTTTGGGGGCGGCTGTCTGCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTgg                                                                  <  1:458797/147‑1 (MQ=255)
                                                                                  tttACCCAACTCCTGGTGAAAAAAGCCCTTGGGGGCGCCTGCCTGCAAGTTTTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTaa                                                               <  1:691768/150‑1 (MQ=255)
                                                                                       caagcgcttgggcaaaaaaTCCCTCCGCGGCGTCTTTCTCCGGGTGTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGTGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                              <  2:20796/132‑1 (MQ=255)
                                                                                         accGCTATGTGAAAAAATGCCTGCGCGGGGGCTTTCTGCCAGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                                          <  2:13366/146‑1 (MQ=255)
                                                                                         accGCCATGTGAAAAAATCCCTCCGCGGCGGCTGTCTCCCATTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                         <  1:366642/147‑1 (MQ=255)
                                                                                           ccccATGTCAAAAATGGCCTGCGGGGCGGTTGTCTGCCAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                                       <  1:682148/147‑1 (MQ=255)
                                                                                             cccttggaaaaaaTGCCCTTCGGGCCGGCTGTCTCCCTGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                                       <  2:434569/140‑1 (MQ=255)
                                                                                                tgtgAAAAAATGCCTGCGGGGGGGTTGTCGCCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                                 <  1:725379/150‑1 (MQ=255)
                                                                                                tggcaaaaaaTGCCTCCGGGGCGGCTTTCTCCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAgg                                                  <  1:274649/144‑1 (MQ=255)
                                                                                                 gtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGAAAAAAAATTTAAACCGCCCTGTTGTGTGGGGATCCGTGCTggggaaaaagtg                                                        >  1:388109/1‑132 (MQ=255)
                                                                                                    cAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAAATTTTAAACCCTCCCGTTTTGTGGGAAGCGCCTCTGGGTAAAAAGTGGGAAggg                                                 >  1:339747/1‑146 (MQ=255)
                                                                                                    aaaaaaTCCTTGCGGGGGGGCTGTCGGCTTGTTTTTTCCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                              <  2:536214/132‑1 (MQ=255)
                                                                                                     aaaaaTCCCTCCGCGGCGGTTTTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                            <  1:186414/150‑1 (MQ=255)
                                                                                                     aaaaaGGCCTCGGGGGTGGCTGTTTGCGAGTTTTTTTGTCGCCAAAAAAAATCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                                 <  2:1052209/145‑1 (MQ=255)
                                                                                                            gCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGAGAGAAAAAACCTTTAAACCCCCCTGTTTTGTGGGGATCGCTGCTGGGGACAAAGTGGGAAGGGGAgggggggg                                       >  1:540473/1‑142 (MQ=255)
                                                                                                                        ggcggcTGTCTGCCAGCTATTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTcc                             <  1:1239325/151‑1 (MQ=255)
                                                                                                                           ggttGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAg                           <  1:695138/147‑1 (MQ=255)
                                                                                                                                tctgcctctttttttcCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                         <  1:646660/136‑1 (MQ=255)
                                                                                                                                tCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCGAGAGAAAAAAATTTTAAACCCCCCTGTTGTGTGGGGATAGCCGCTGGGGAAAATATGGGAAAGGGAGGGGGGGGGATCAGCTCCCGaa                         >  1:972720/1‑147 (MQ=255)
                                                                                                                                      gttttttttcCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                       <  1:1103892/138‑1 (MQ=255)
                                                                                                                                         tATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt               <  2:466298/148‑1 (MQ=255)
                                                                                                                                               tCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGAAGAGGAAAAAAATTTAACCCGGCATGTTGTGGGGGAAAAGCTGGTGGGAAAAAAATGGGAAAGGGAGAGGGGGGGAACAACCGCACAAACCCGCTAATTGAGtttg       >  1:601971/1‑148 (MQ=255)
                                                                                                                                                    ccgcAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCa  <  1:393123/151‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                                    
GATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCA  >  NC_000913/357742‑358034

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