Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913231,6940CT45.2% 72.0 / 29.6 62F191F (TTC→TTT) yafDendo/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base C (7/27);  new base T (28/0);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 2.59e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGA  >  NC_000913/231554‑231843
                                                                                                                                            |                                                                                                                                                        
ccTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGa                                                                                                                                               <  2:186512/151‑1 (MQ=255)
      acCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACCTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTCAATGCCT‑gggg                                                                                                                                              >  2:801502/1‑144 (MQ=255)
         cGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCtcggggcggag                                                                                                                                         >  1:346934/1‑139 (MQ=255)
             tGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGAATTTAAAGCCT‑GGGGCCGGAgga                                                                                                                                      >  1:979763/1‑145 (MQ=255)
              gTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGTTTCAATGCCC‑GGAGCCGGAGAA‑gg                                                                                                                                   >  2:935077/1‑148 (MQ=255)
              gTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAATTTAATGCCT‑TGGGCCGTAGaaag                                                                                                                                    >  1:801111/1‑146 (MQ=255)
                   tggtggTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAATTTAATAGCT‑GGGGCCGTTGAA‑AGagt                                                                                                                                >  1:103940/1‑144 (MQ=255)
                    ggtggtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTTATGCCC‑GGAGGCGTAGaa                                                                                                                                      >  2:275985/1‑140 (MQ=255)
                    ggtggtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAAcgcg                                                                                                                          <  1:579158/151‑1 (MQ=255)
                    ggtggtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGCCGATTTTCATGACC‑GGGGCCGCAGAA‑GGacggaag                                                                                                                            >  2:81324/1‑143 (MQ=255)
                       ggtTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGt                                                                                                                         <  2:1130137/149‑1 (MQ=255)
                          tAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAATTCTATCCCC‑GGGGCCGGAGAA‑GGATGaa                                                                                                                              >  1:1041537/1‑141 (MQ=255)
                           aaTATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTTATGCCT‑GGGGGCGCAGAA‑GGAAGAAAGCGtt                                                                                                                        >  2:782708/1‑146 (MQ=255)
                             tataCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGGTTTAATTCCT‑TGAGGCCGAGAAGGGATGAAGGGGtttttg                                                                                                                    >  1:628706/1‑145 (MQ=255)
                              ataCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGATTTAATGCCT‑TGGGCCGTAGAA‑AGATGAAAGCGTTATATCGCt                                                                                                                >  2:807060/1‑151 (MQ=255)
                              ataCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCT‑GGGGGCGTAGGA‑GGAAGAAAGCGGTATATTGCt                                                                                                                >  2:732615/1‑151 (MQ=255)
                               taCACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCT‑GGGGCCGGAGAA‑GGATGAAAGCGGTTTATTGCtt                                                                                                               >  2:114712/1‑151 (MQ=255)
                                 cacGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTg                                                                                                             <  1:541375/151‑1 (MQ=255)
                                  acGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGGGCCGGAGGA‑GGATGAACGCGGTATATTGCttt                                                                                                              >  2:1013079/1‑149 (MQ=255)
                                                     gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCT‑GGGGGCGCAGAA‑GGGGGAACGCGGTTTATTGCTTTTCGCGGGGAAAGTTGCTg                                                                                            >  2:21996/1‑148 (MQ=255)
                                                     gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTgcg‑c                                                                                        <  2:500805/151‑1 (MQ=255)
                                                      ggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCC‑GGGGCCGCAGAA‑GGGTGAAAACGGTTTTTCGCCTTTCGCCGGGAAAGTCg                                                                                               >  1:931690/1‑144 (MQ=255)
                                                       gCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCa                                                                                      <  2:356626/151‑1 (MQ=255)
                                                         gTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAgg                                                                                    <  2:861435/151‑1 (MQ=255)
                                                           ggATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTg                                                                                  <  2:773699/151‑1 (MQ=255)
                                                            gATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGAATTTAATTCCC‑TGAGCCCGAAAA‑GGAAGAACGCGGTATATAGCTTTTCGCCGGGAAAGTCGCCGCG‑CCCgggg                                                                                  >  1:120244/1‑148 (MQ=255)
                                                                 cTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAGTTTAAAGCCT‑GGGGCCGGAGAA‑GGATGGACGCGGTATATCGCTTTTCGCGGGGAATGTCGCCGCG‑CCCGGGGCGgttt                                                                            >  1:694900/1‑151 (MQ=255)
                                                                  tataGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTTATTCCT‑GGAGCCGGAGAA‑GGATGAAAGCGGTATATTGCTTTTCGCGGGGAAAGTTGCCGCG‑CCCgggggggtttt                                                                           >  1:58967/1‑142 (MQ=255)
                                                                          gCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCgatgat                                                                   <  1:873822/151‑1 (MQ=255)
                                                                           cAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTCAATGCCT‑GGAGCCGTAGAA‑AGAAGAACGAGGTATATAGCTTTTCGCGGGGAAAGTCGCCGCG‑CCAgggggggttttcc                                                                         >  2:671828/1‑133 (MQ=255)
                                                                             gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCT‑GGAGCCGTAGAA‑GGAAGAACGCGTTATATCGCCTTTCGCGGGGAATGTCGCCGCG‑CCAGGGGGGTTTTTCCgatga                                                                    >  2:521103/1‑147 (MQ=255)
                                                                              ttACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCT‑GGGGCCGGAGAA‑GGGTGAAAGCGTTATTTCGCTTTTCGCGGGGAATGGCGCTGCG‑CCCGGGGCGGTTTTCCGATGATCAg                                                                >  2:212148/1‑150 (MQ=255)
                                                                                aCTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAgcg                                                              <  2:891679/150‑1 (MQ=255)
                                                                                    ccTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCT‑GGAGCCGTAGAA‑GGAGGAACGCGGTATATCGCTTTTCGCGGGGAAAGGCGCCGCG‑CCCGGGGGGGTTTTCCCATGATCAGCGcc                                                            >  1:1318/1‑148 (MQ=255)
                                                                                      tATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGATTTAATGCCC‑GGGGCCGTAGAA‑GGATGAAAGCGGTATATTGCTTTTCGCGGGGAAAGTTGCTGCG‑CCAGGGGGGGTTTTCCGGTGATCAGCGCCGTCgg                                                       >  1:612458/1‑151 (MQ=255)
                                                                                       aTTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAAAGCCT‑GGAGCCGTAGAA‑GGATGAAAGCGGTATATTGCTTTGCGGGGGGAAAGTTGGTGCG‑CCAGGGGGGGTTTTACGGTGATTAGGGCCGGCggg                                                      >  2:74346/1‑151 (MQ=255)
                                                                                       aTTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCgg                                                       <  1:902182/150‑1 (MQ=255)
                                                                                       aTTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCgg                                                       <  2:650974/150‑1 (MQ=255)
                                                                                       aTTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCg                                                        <  2:965990/149‑1 (MQ=255)
                                                                                         tGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGc                                                     <  1:488173/150‑1 (MQ=255)
                                                                                          ggCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCg                                                    <  2:1120549/150‑1 (MQ=255)
                                                                                            cGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGATTTAATACCTCGGGGCCGCGGAA‑GG‑GGAAAGCGGTTTATTCCTTTTCGCGGGGAAAGTCGCCGCG‑CCCGGGGCGGTTTTCCgatg                                                                     >  1:598794/1‑131 (MQ=255)
                                                                                             gATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTg                                                <  1:316153/151‑1 (MQ=255)
                                                                                             gATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTg                                                <  1:801474/151‑1 (MQ=255)
                                                                                                  gATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAGTTCCATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGGTATATCGCTTTGCGCGGGGAATGTCGCTGCG‑CCCGGGGCGTTTTTCCGCTGATCAGCGCCGTCGGGCGTTTGGTCg                                            >  2:1011915/1‑150 (MQ=255)
                                                                                                      gCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCT‑GGAGGCGTAGAA‑GGATGAACGCGGTATATCGCTTTGCGCGGGGAATGTCGCTGCG‑CCAGGGGGGTTTTATCGATGATCAGCGCCGGCGGGCGGTTGGTCGccc                                         >  2:105852/1‑149 (MQ=255)
                                                                                                       cTCACCACAGCGGCCCGGTCATTATGGCGGGGGAGTTTAAAGCCC‑GGGGCCGGCGAA‑GGGAGGAAGCGTTATATTGCTTTGCGCGGGAAAAGGCGCTGAGCCCAGGGGGGTTTTACGGA‑GATAAGCGCCGGCCGGGGGTTGGTCggc                                          >  2:545562/1‑145 (MQ=255)
                                                                                                         caccacAGCGGCCCGGTCATTATGGCGGGGGATTTTAAATCCT‑CGAGGCCGAGAA‑GGATGAACGCGTTTTATCGCTTTTCGCGGGGAAAGTCGCTGCG‑GCAGGGGCGGTTTTCCCGTTATCAGCGGCGGCGGGCGTTTTGTCGCCCGCTc                                     >  1:24009/1‑150 (MQ=255)
                                                                                                          accacAGCGGCCCGGTCATTATGGCGGGGGGTTTTTATTCCC‑CGAGGCGGGGGA‑GGGAGAAAGCGGTATATTGCTTTGTGCGGGGAAAGTCGCTGCG‑CCCGGGGGGGTTTTCCGATGGTCTGCGCCGTCGGGGGGTTGGTCGccc                                         >  1:46589/1‑145 (MQ=255)
                                                                                                          accacAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCTCGGGAAATGTCGCTGCA‑CCGGGTGCGTTTTACCGATGATCAGCCCCGTCGGGCGTTTGGTCTCCCGCTCg                                    <  1:937449/150‑1 (MQ=255)
                                                                                                            cacaGCGGCCCGGTCATTATGGCGGGAGATTTTAAAGCCC‑CGGGGCGGAGAA‑GGGAGAAAGCGGTTTATCGCTTTTCGCGGGGAAAGTTGCTGCG‑CCCGGGGCGGTTTTCCGAAGATCAGCGCCCCCGGGGGGTTGGGCGgccccccggt                                  >  1:50928/1‑143 (MQ=255)
                                                                                                               aGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTg                              <  2:897043/151‑1 (MQ=255)
                                                                                                                  ggCCCGGTCATTATGGCGGGAGATTTTAATTCCC‑GGAGCCGGAGAA‑GGATGAACGCGTTATATTGCTTTTCGCGGGGAATGTCGCTGCG‑CCAGGGGCGTTTTTTCGATGATCAGCGCCGGCGGGCGGTTTGGCGCCCCCCCGAtttttttt                           >  1:959235/1‑151 (MQ=255)
                                                                                                                    cccGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGtttt                          <  1:419200/150‑1 (MQ=255)
                                                                                                                       ggTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTAc                      <  2:467552/151‑1 (MQ=255)
                                                                                                                          cATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCg                    <  1:111094/150‑1 (MQ=255)
                                                                                                                          cATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGt                   <  2:431952/151‑1 (MQ=255)
                                                                                                                           aTTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTTGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTg                  <  2:1099040/151‑1 (MQ=255)
                                                                                                                              aTGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTc               <  2:971176/151‑1 (MQ=255)
                                                                                                                                ggCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCt              <  1:1091390/150‑1 (MQ=255)
                                                                                                                                 gCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTg             >  2:561209/1‑150 (MQ=255)
                                                                                                                                     gagaTTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGt        <  1:233408/151‑1 (MQ=255)
                                                                                                                                            cAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGa  >  2:193879/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                        
CCTGACACCCGCCTGTTGATGGTGGTTAATATACACGCCGTCAACTTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCT‑GGAGCCGTAGAA‑GGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCG‑CCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGA  >  NC_000913/231554‑231843

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