Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,8900CA48.7% 47.7 / 60.3 78L267F (TTG→TTT) cynRtranscriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (28/12);  new base A (1/37);  total (29/49)
Fisher's exact test for biased strand distribution p-value = 1.11e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCG  >  NC_000913/357740‑358029
                                                                                                                                                          |                                                                                                                                                
cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGcaac                                                                                                                                                  >  2:560089/1‑151 (MQ=255)
 tGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGcaac                                                                                                                                                  >  2:1489631/1‑150 (MQ=255)
   aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGcaac‑aa                                                                                                                                               >  1:1510675/1‑150 (MQ=255)
   aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑aaa                                                                                                                                              <  2:965819/151‑1 (MQ=255)
     aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑aaaa                                                                                                                                             <  1:548192/150‑1 (MQ=255)
     aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAAc                                                                                                                                            >  1:1090089/1‑151 (MQ=255)
     aaGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAAc                                                                                                                                            <  1:1302530/151‑1 (MQ=255)
      aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAAcc                                                                                                                                           <  2:570887/151‑1 (MQ=255)
       gCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAAccg                                                                                                                                          >  1:1742135/1‑151 (MQ=255)
         gtggcgcATCAGGCAATTCCAGCCGCAGACCTGTGCCAGCGGCTACCGTGATTCATTCCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAAccgcc                                                                                                                                        <  1:1019065/148‑1 (MQ=255)
               cATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAATTGCTTT‑CGGGCCGGTGGCCTCCA‑AGC‑AATTTTCCGGCGCAAA‑CAAAACGCCGGtc                                                                                                                                    >  2:516288/1‑148 (MQ=14)
                aTCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTtctct                                                                                                                                 >  2:1840328/1‑151 (MQ=255)
                 tCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTtctctc                                                                                                                                >  1:1741195/1‑151 (MQ=255)
                   aGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTcc                                                                                                                               >  2:164634/1‑150 (MQ=255)
                   aGGCAATTCCAACCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACTACCGCACATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAATCGCCGTTCTCTcc                                                                                                                               >  1:1420711/1‑150 (MQ=255)
                         tgccAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                        >  1:275211/3‑151 (MQ=255)
                           ccccccccAGACTTGTGCCAGCGCCTCCCGTATTTCTTTCCCCCC‑AACACCCGCGCTTTTCCCCAC‑CCCCTTGTGAAAAAAGCCCTT‑CGGGGCGGCTGTCTCCC‑AGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCag                                                                                                                             <  1:980264/137‑1 (MQ=255)
                           ccAGCCGCTGTCCTGTTCCAGCGGCTCCCTTATTTCTTTCCCCCA‑AACACCCGCGCATTTCCCCAC‑CCCATTGTAAAAAAATGCCTC‑CGCGGCGGCTGTCTCCC‑AGCTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGc                                                                                                                      <  1:167275/151‑1 (MQ=255)
                           ccAGCCGCAGCCCTGTGCCAGGGCCCCCCGTTATTCTTTCCCCCC‑AACACCCGCGCTTTTATCCAC‑CCCCTTGTGAAAAAAGGCCTC‑CGCGGCGGCTGTCTGCA‑AGTTATTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                        <  2:392517/149‑1 (MQ=255)
                           ccAGCCGAAGACTTGTGTCAGGGGCTCCGGGGATTCATTCCCGCA‑AACACCCGCGCTTTTACCCAA‑CGCCTTGTGCAAAAATGCCTC‑CGGGGGGGCTGTCTCCC‑AGTTATTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                        <  2:933101/149‑1 (MQ=255)
                            cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑g                                                                                                                   >  1:358664/1‑151 (MQ=255)
                            cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑g                                                                                                                   >  2:1054622/1‑151 (MQ=255)
                            cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑g                                                                                                                   >  2:289151/1‑151 (MQ=255)
                            cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑g                                                                                                                   >  2:620202/1‑151 (MQ=255)
                               ccacagccgtgtgCCAGCGGCTCCCGGGATTCTTTCCCCCC‑AACACCCGCGCTTTTACCCAC‑CGCCATGTGCAAAATTGCCTT‑CGCGGCGCCTGTCTGCC‑AGTTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑gggg                                                                                                                <  1:74405/142‑1 (MQ=255)
                               ccGCAGCCCTGTGTCAGCGGCTACCGTGATTCTTTCCCGCC‑AACACCCGCGATTTTACCCAC‑CGCCATGTGAAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGTTTTTTTTCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑ggg                                                                                                                 <  2:1664158/150‑1 (MQ=21)
                                cgcacccgtgttTCAGCGGCTCCCTTTTTTCTTTTCCCCC‑AACACCCCCCCATTTCCCCAA‑CGCCTTTTGAAAAACTGCCTC‑CGGGGCGGTTTTCTGCC‑AGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑gggg                                                                                                                <  2:490702/142‑1 (MQ=255)
                                cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑ggg                                                                                                                 >  1:1840195/1‑149 (MQ=255)
                                  cAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑ggcgg                                                                                                               >  2:1385238/1‑146 (MQ=255)
                                   aGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑ggggg                                                                                                               >  2:1828828/1‑147 (MQ=255)
                                    gACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑gggg                                                                                                                >  2:276482/1‑146 (MQ=255)
                                      ccTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGagag                                                                                                          >  1:679009/1‑150 (MQ=255)
                                       cTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGCAgaga                                                                                                         >  2:1090180/1‑150 (MQ=255)
                                        tgtgTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑gggg                                                                                                                >  1:1927349/1‑142 (MQ=255)
                                              aGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑gggggcgaggagaaaaa                                                                                                   >  1:1050189/1‑136 (MQ=255)
                                                cGGCTCCCGTGTTTCATCCCCCCC‑AACACCCGCGCTTTTCCCCAA‑CCCCTTGTGAAAAAATGCTTG‑CGGGGCGGCTGTCTCCC‑AGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGaa                                                                                                        <  1:1963481/142‑1 (MQ=255)
                                                 ggCTCCCGTAATTCATCCCCCCC‑AACACCCGCGCTTTTACCCAC‑CCCCATGTGAAAAAATCCCTC‑CGCGGCGGCTGTCGGCA‑AGCTTTTTTCCCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTa                                                                                              <  1:1714806/150‑1 (MQ=255)
                                                 ggCTACCGTGATTCTTTCCCGCC‑AACACCCGCGCATTTACCCAC‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTTTCTGCC‑AGTTATTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTaa                                                                                             <  2:90508/151‑1 (MQ=25)
                                                     ccctgttttCTTTCCCGCC‑AACACCCG‑GTATTTTACCAACCCCCATGTGAAAAAAGCCCTC‑CGCGGCGCCTTTCTGCA‑AGTTTTTTTCCCCCCCAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGccc                                                                                         <  1:1029627/144‑1 (MQ=255)
                                                      ccgtatttCATTCCCCCC‑AACACCCGCCCTTTTACCCAC‑CCCCTTGTGAAAAATGGCCTC‑CGCGGCGGCTGTCTGCC‑AGTTATTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGc                                                                                           <  1:885827/142‑1 (MQ=255)
                                                       ccggatttctttCCCCCAAACACCCGCGTTTTTACCCAC‑CCCCAGGTAAAAAATTCCTTT‑CGGGGCGGCTGTCGCCC‑AGTTTTTTTCCCGCCCAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGc                                                                                           <  1:1907464/139‑1 (MQ=255)
                                                           aTTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAGAGaaaaaaattttaaacc                                                                                          >  2:760795/1‑132 (MQ=25)
                                                             tCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAAAAAC‑CTTAAAACCCTCCtg‑tt                                                                                >  2:984516/1‑151 (MQ=25)
                                                             ccATTCCCCCA‑AACACCCCCGAATTTACCAAC‑CGCAATGT‑AAAAAAAGCCCT‑CGGGGGGGGCTTCCGGC‑AGGTTTTTTTCCCCCCAAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCAt                                                                                    <  1:1008859/147‑1 (MQ=255)
                                                                 ttCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAGAGAAAAAAC‑CTTTAACCCGTCCCT‑Ttgtgt                                                                            >  1:1181364/1‑151 (MQ=25)
                                                                   ccccc‑AACACCCGCGCATTTCCCAAC‑CGCCATGTAAAAAAATCCCTT‑CGGGGGGGCTGTCGGCC‑AGTTTTTTTCCCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCAtg‑t                                                                                 <  1:605189/141‑1 (MQ=255)
                                                                   ccccc‑AACACCCGCCCATTTACCCAC‑CGCCATGTGAAAAAATCCCTT‑CGCGGCGGCTTTCTGCC‑AGCTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑Ttgtgt                                                                            <  2:1397849/146‑1 (MQ=255)
                                                                      cc‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGGCAAG‑GAAAAAGA‑TTTAAGCCGGTc                                                                                      >  2:1384590/1‑138 (MQ=25)
                                                                         aacaccCGCGCTTTTACCCAC‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTTTTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATc                                                                     <  2:756120/151‑1 (MQ=38)
                                                                                gcccattaccccAA‑CGCCAGGTCAAAAAAGCCCTG‑CGCGGCGGCTGCCGGCC‑AGTTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCgctgct                                                               <  1:606108/140‑1 (MQ=255)
                                                                                       tATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAAAGAAAAAACATTTAAACCC‑CCATG‑TTTTGGGGGGATCGCTGCTGGtgaaa                                                        >  2:397757/1‑146 (MQ=25)
                                                                                          ccAC‑CGCCATGTGAAAAAATGCCTC‑CGCGGCGGCTGTCTGCC‑AGTTGTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                    <  2:1963078/151‑1 (MQ=255)
                                                                                           ccacccccgtgtaaaaaaaGCCCTC‑CGCGGCGGCTGTCTGCC‑AGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAAtg                                                     <  2:1085837/136‑1 (MQ=255)
                                                                                            aa‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAGCGAAAAAAA‑ATTAAACCCCCCCTG‑TTGTTGGGGCATCGCTGCTggggaca                                                        >  1:236927/1‑140 (MQ=21)
                                                                                             acccccATGTGC‑AAAATGCCTTCCGCGGCGTCTGTCGGCC‑AGTTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                                  <  1:1268688/147‑1 (MQ=255)
                                                                                               cccaaTGTGAAAAAATGCCTC‑CGCGGCGGTTGTCTGCC‑GGCTATTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGa                                                 <  2:1474553/146‑1 (MQ=255)
                                                                                                 ccATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGAAGAGAAAAAAC‑CTTAATCCCCGCCTGCTTTTGTGG‑AATCGgtgctgtggggtcaat                                                    >  1:854643/1‑135 (MQ=25)
                                                                                                  ccgtgtgaaaaaaGCCTTCCGGGGCGGCTGTCGGCC‑AGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAg                                              <  1:1673284/144‑1 (MQ=255)
                                                                                                     gtgAAAAAATCCCTT‑CGCGGCGGCTGTCTGCC‑AGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGaa                                          <  2:1717033/151‑1 (MQ=255)
                                                                                                         aaaatTCCCTC‑CGGGGCGCTTGTCCCCC‑AGTTTTTTTCCCGCCACAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTg                                       <  1:1456434/150‑1 (MQ=255)
                                                                                                         aaaaaTTCCTT‑GGCGGGGGCGGGCTGCCGACCTATTTTCCCCCGGAAAA‑AAACCCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGaaa                                               <  2:815717/143‑1 (MQ=255)
                                                                                                         aaaaaTCCTTC‑CGGGGGGGCTTTCTCCC‑GGTTTTTTTCCCGCCAAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgc                                      <  1:722030/151‑1 (MQ=255)
                                                                                                         aaaaaTCCCTG‑CGCGGCGCCTTTCTGCC‑GGTTATTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGt                                        <  1:507422/149‑1 (MQ=255)
                                                                                                          aaaaTGCCTT‑CGCGGGGGCTGTCTCCC‑ATCTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcg                                     <  1:1443590/151‑1 (MQ=255)
                                                                                                                ccTT‑CGCGGCGGTTGTCTGCC‑AGCTTTTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggcg                                  <  2:394529/148‑1 (MQ=255)
                                                                                                                ccTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAGAGAAAAAAA‑TTTTAAACCCCCCTT‑TTTTGTGGGGATTGCCGGTGGGAAAAAAGTGggaagggaaggggggg                                   >  1:1902162/1‑141 (MQ=255)
                                                                                                                   tccggggcgccTGTCTGCC‑AGTTATTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATc                              <  2:1087478/145‑1 (MQ=255)
                                                                                                                   t‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGGGAGAGAAAAAAA‑TTTTAACCCCCCCTG‑TTTTGTGGGGAATGCCGCTGGGTACAAAGTGGGAAGGGGAgggggggggaac                              >  2:1811705/1‑138 (MQ=255)
                                                                                                                    ccgcgGCGCTTTTCTGCA‑GGTTTTTTTCCCGCCAAAAA‑AAATCCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAAt                               <  1:226325/147‑1 (MQ=255)
                                                                                                                        ggcggcTTCCTGCC‑ATTTTTTTTCCCGCCGAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCt                          <  1:1389652/149‑1 (MQ=255)
                                                                                                                           gcgtttCTGCC‑GGTTTTTTTCCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGc                           <  1:470864/140‑1 (MQ=255)
                                                                                                                                 cggCC‑ACCTTTTTCCCCGCCGCACA‑AACACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                  <  2:1002404/146‑1 (MQ=255)
                                                                                                                                  tccc‑AGTTTTTTTCCCGCCACAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                  <  2:572098/145‑1 (MQ=255)
                                                                                                                                  ctcccggtttttttcCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAc                   <  1:1233884/137‑1 (MQ=255)
                                                                                                                                    cc‑AGTTTTTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg                 <  2:553765/146‑1 (MQ=255)
                                                                                                                                    cc‑AGTTTTTTTCCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt             <  2:772334/150‑1 (MQ=255)
                                                                                                                                     ccgtcttttttcCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCt             <  1:82493/143‑1 (MQ=255)
                                                                                                                                     ccggtttttttcCCCCCGAAAA‑AAACCCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc                <  1:18721/140‑1 (MQ=255)
                                                                                                                                       agtttttttcCCGCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAAt          <  1:331980/146‑1 (MQ=255)
                                                                                                                                            ttttcCCCCCGCAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa       <  1:356676/149‑1 (MQ=255)
                                                                                                                                                 ccgccgAAAA‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGTGAATCAGCTCCAGAACCGCGCTAATTGAGTTCg  <  1:907791/149‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                                
CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCC‑AACAACCGCGCATTTATCCAA‑CGCCATGTGCAAAAATGCCTT‑CGCGGCGGCTGTCTGCC‑AGCTATTTTTCCGCCGCAAC‑AAAACCGCCGTTCTCTCCAGTAGTGGC‑‑GGGGCAAGAGAAATAGC‑TTTAAGCCCGTCATG‑TTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCG  >  NC_000913/357740‑358029

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