Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8750AT49.2% 37.0 / 37.5 64N272K (AAT→AAAcynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (24/8);  new base T (0/31);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 8.39e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG  >  NC_000913/357731‑358023
                                                                                                                                                       |                                                                                                                                                       
cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                                                                                                                  <  1:413633/151‑1 (MQ=255)
cTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                                                                                                                  <  2:1037507/151‑1 (MQ=255)
 tCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTc                                                                                                                                                  >  1:956012/1‑150 (MQ=255)
      ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                            >  1:1324960/1‑150 (MQ=255)
      ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTccgccg                                                                                                                                            >  2:735856/1‑150 (MQ=255)
      ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTcccccg                                                                                                                                            >  1:1065075/1‑150 (MQ=255)
           gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                      >  2:1553464/1‑151 (MQ=255)
           gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGcaacaa                                                                                                                                      >  1:569531/1‑151 (MQ=255)
           gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGAGATTCATTTCCGCCAACAACCGCGAATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCCGCTATTTTTCCGCCGcaacaa                                                                                                                                      >  2:578935/1‑151 (MQ=255)
               aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAc                                                                                                                                   >  1:415924/1‑150 (MQ=255)
                gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccg                                                                                                                                 >  1:775338/1‑151 (MQ=255)
                             ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCag                                                                                                                    >  1:505602/1‑151 (MQ=255)
                             ggCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCag                                                                                                                    >  1:195687/1‑151 (MQ=255)
                                attcccggCCGAAGCCCGGTGCAAGCGCCTCCCTTAATTCTTTCCCCCCAACACCCGCTCTTTTCCCCAACCCCATTTGAAAAATGCCCTCGGCGGGGGCTGTCGGCCATTTTTTTTCCCGCCCCAACAAAACCGCCGTTCTCTCCagtag                                                                                                                 <  1:1327964/144‑1 (MQ=255)
                                attccccgCCGAAGACCTGTGTCGGGGCCTCCCGTGATTCATTCCCCCCAACACCCGCGCTTTTACCCAACCCATTGTGAAAAAAGGCCTCCGGGGGGGTTTTCTGCCAGCATTTTTCCCGCCCAAACAAAACCGCCGTTCTCTCCagtag                                                                                                                 <  1:958434/144‑1 (MQ=255)
                                atccccACCCGCAGCCGTGTGCCGGCGCCTACGGTGTTTCATTCCCCCCAACACCCGCGATTTTACCCACCCCATTGTCCAAAAATCCCCTCGGGGCGGCTGTCTGCAGGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCagtag                                                                                                                 <  2:371554/147‑1 (MQ=255)
                                actcccggCCGCAGCCTTGTGCCGGCGGCTACGGTGATTCATCCCCCCCAACCCCCGCGCTTTCCCCCAACGCCATGTGAAAAAAGGCCTCCGCGGCGGCTGTCGGCAGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagtag                                                                                                                 <  2:44042/144‑1 (MQ=255)
                                aaccccAGCCGAACCCTTGTGACGGCGGCCACTGTGATACTTCCCCGCAAACAACCGCTCTTTTACCAAACGCATTGTAAAAAAATCCCTCCGCGGCGGCTGCCTGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCagtag                                                                                                                 <  2:1090162/147‑1 (MQ=255)
                                aaTCCCAGCCGCAGACTGGGTCCAGCGCCACCCGGAATTCATTCCCGCAAACACCCGCGCTTTTCCCAACCGCCTGGTGCAAAAATCCCTCCGCGGGGGCTTTCGGCAGGTTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCagtag                                                                                                                 <  2:368585/151‑1 (MQ=255)
                                 aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                >  1:740216/1‑151 (MQ=255)
                                 aTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagtagt                                                                                                                >  1:1623062/1‑151 (MQ=255)
                                        ccGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggg                                                                                                         >  2:281465/1‑151 (MQ=255)
                                             gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcagg                                                                                                    >  1:1106021/1‑146 (MQ=255)
                                             gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGcgaca                                                                                                    >  2:402903/1‑147 (MQ=255)
                                              aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggtcgag                                                                                                   >  2:802305/1‑145 (MQ=255)
                                               cttgtgTCGGCGGATCCCGTGATTCTTTCCCCCCAACACCCGCCCTTTTACCCACCGCCAGGTGAAAAAATGCCTTCGGGGGGGCTTTCTGCAAGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAgaga                                                                                                 <  1:509576/149‑1 (MQ=255)
                                               ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcaggga                                                                                                  >  1:21410/1‑144 (MQ=255)
                                                  gtgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCggggggcgaggaaa                                                                                               >  2:355973/1‑141 (MQ=255)
                                                  gtgtCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCgggggcaggggaaa                                                                                               >  1:452466/1‑141 (MQ=18)
                                                      ccggCGTCTACGGTATTTCATTCCCCCCAACACCCGGCCATTACCCCACCGCGATGGCAAAAAAGGCTTGGGGGGGGTGGGTCGACCATCTGTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAg                                                                                            <  1:337003/147‑1 (MQ=255)
                                                              tcccgtttttCTTTCCCCCCAACACCCGCGATTTTTCCAACCGCAAGGTAAAAAATGGCTTCCGGGGGGGCTGCCTCCCAGTTTTTTTCCCGCGGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                                                                   <  1:411282/143‑1 (MQ=255)
                                                              tACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGGAAAAACTTTaaaa                                                                                    >  1:1233587/1‑148 (MQ=25)
                                                              ccctgtgtTTCATTCCCCCAAACACCCGGCTTTTTACCAACCCCCGTGTGAAAAAATCCTTCCGGGGGGGCGGCCTCCCGGTTTTTTTCCCGCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                                                                   <  2:980553/147‑1 (MQ=255)
                                                              cccggtatttCATCCCCGCAAACACCCGCGTTTTTCCCCACCCCCTGGAGAAAAATTCCCTCCGCGGCGGTTGTCGCCCGGTTTTTTTCCCGCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGcc                                                                                   <  2:602424/143‑1 (MQ=255)
                                                               aCCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAaaaaaaacattttacc                                                                                   >  1:1366820/1‑138 (MQ=25)
                                                                     aTTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGTAAGAGAAACAACTTTaacaccccgcgt                                                                            >  2:1404033/1‑141 (MQ=25)
                                                                      ttcttcccccacaacaaCCGCTATTTAACCAAACCCCGTGGGCAAAACTGCGTGGGGGGGGGCTGTCGGCAGGTTTTTTTCCCGCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgt                                                                           <  2:863313/140‑1 (MQ=255)
                                                                       tctttccccccaccaCCCGCCCATTTACCCACCCCCTTGTGCAAAAATCCTTCCGCGGCGGCTGTCCCCCAGCTTTTTTCCCGCCCCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgtt                                                                          <  1:763080/142‑1 (MQ=255)
                                                                        ttttccccccAACCCCCGCGCTTTTCCCCAACCCCATGTAAAAAAATGCCTTCGGGGGGGCTGTCTGCCAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                         <  1:347240/143‑1 (MQ=255)
                                                                        ctttccccccaccaCCCGCGCTTTTACCAACCGCCATGTGAAAAAATCCCTCCGCGGCGGCTGTCTGCGAGTTATTTTCCCGCCGCAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                         <  2:191470/143‑1 (MQ=255)
                                                                        ctttccccccAACACCCGCGCTTTTCCCCAACGCCATGTAAAAAAATGCCTTCGCGGGGGCTTCCTCCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                         <  1:199571/143‑1 (MQ=255)
                                                                        catcccccccAACACCCGCGCTTTTACCCACCGCCAGGTAAAAAAATCCCTGCGGGGCGGCTGTCTCCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCAtgttg                                                                         <  2:78311/143‑1 (MQ=255)
                                                                             ccaccAACACCCGCGCTTTCACCCAACCCCATGTAAAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTCCCGCCACAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTgg                                                                    <  2:1529981/148‑1 (MQ=255)
                                                                              ccccaccaCCCGCGCTTTCCCCCACCCCCTTGTGAAAAAATGCCTTCGGGGCGGCTGTCTGCCAGCTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                                   <  2:499741/149‑1 (MQ=255)
                                                                                              tttACCCACCGCTATGAAAAAAATGGCTTTCGGGGCGGCTGTCGGCCTGCTTTTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                   <  1:176602/151‑1 (MQ=255)
                                                                                              tttACCCAACCCCTTGTAAAAAAAGCCTTCCGGGGCGCCTGTCTGCCAGTTTTTTTCCCGCCACAACAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                   <  2:597374/151‑1 (MQ=255)
                                                                                               tttaccAACCCCATTGGAAAAAAAGCCTTGCGCGGCGGTTGGATGCCAGTTTTTTCCCGCGCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                   <  2:1336393/147‑1 (MQ=255)
                                                                                               ttATCCAACCCCCTGGGAAAAAAGGCCCTTCGCGGGGCCTGTCGCCGAGTTTTTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAATGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAAc                                                   <  1:1366406/151‑1 (MQ=255)
                                                                                                  cccccaacccatggtgCAAAACTCCTTCCGCGGCGCCTTTCTGCTAGTTTTTTTCCCGCCAAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa                                                 <  2:727415/138‑1 (MQ=255)
                                                                                                    cccaccccctgggaaaaaaaGCCCTCCGGGGCGGTTTTCGCCCGCTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtg                                               <  1:706090/136‑1 (MQ=255)
                                                                                                             ttgtcaaaaaaCCCCTCCGCGGCGGCTTTCTGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAggg                                      <  2:1206674/145‑1 (MQ=255)
                                                                                                              tgtcaaaaagTGCTTCCGGGGGGGTTTTCGGCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                     <  1:1598891/146‑1 (MQ=255)
                                                                                                              tgcgaaaaagTCCCTCCGCGGCGGCTGTCTGCCAGCTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                     <  2:166928/146‑1 (MQ=255)
                                                                                                              tgaacAAAAGTCCCCTCGCGGCGGCTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGa                                     <  2:109731/147‑1 (MQ=255)
                                                                                                               gtgCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGAGAAAAAAACTTTAACCCCGCCTGTTTTGTGGGCGACGCCGCCGGTTAAAAAGTGGGAAGGGag                                    >  2:1399586/1‑150 (MQ=255)
                                                                                                                   aaaaaTCCCTTCGGGGGGGTTTTCGCCCAGTTTTTTTCCCGCCCAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                <  1:677849/151‑1 (MQ=255)
                                                                                                                                          ttgtctccaatttATTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  1:357415/139‑1 (MQ=255)
                                                                                                                                          ttgcctgcagtttATTTTCCCCCCAAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  2:441985/139‑1 (MQ=255)
                                                                                                                                          ttGTCTCCAAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  1:1270052/150‑1 (MQ=255)
                                                                                                                                          ctgcctgcGAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  2:1403422/151‑1 (MQ=255)
                                                                                                                                          cTGTCTGCAGGTTTTTTTCCCGCCGAAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  2:1105498/151‑1 (MQ=255)
                                                                                                                                          cTGTCTCCCATTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCGAGAAc             <  1:1053335/151‑1 (MQ=255)
                                                                                                                                          cTGTCGGCCAGTTTTTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc             <  1:844665/151‑1 (MQ=255)
                                                                                                                                                     cTATTTTCCCGCCGAAAAAAACCCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg  <  2:1077036/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                       
CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCGGCGGCTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGCAAGAGAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTG  >  NC_000913/357731‑358023

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: