Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,222,2330GA45.0% 64.1 / 37.5 60pseudogene (60/75 nt)ycgHpseudogene; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (26/7);  new base A (0/27);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 8.34e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAA  >  NC_000913/1222083‑1222377
                                                                                                                                                        |                                                                                                                                                 
gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTg                                                                                                                                                   <  2:235432/151‑1 (MQ=255)
  gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGaa                                                                                                                                                 <  1:832469/151‑1 (MQ=255)
    gAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat                                                                                                                                               >  2:828703/1‑151 (MQ=255)
     aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat                                                                                                                                               >  2:263296/1‑150 (MQ=255)
      aaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAatat                                                                                                                                             <  2:84296/151‑1 (MQ=255)
        gTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATaa                                                                                                                                           <  1:516073/151‑1 (MQ=255)
             aCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCt                                                                                                                                      >  2:310581/1‑151 (MQ=255)
               aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGa                                                                                                                                    <  2:697156/151‑1 (MQ=255)
               aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGa                                                                                                                                    <  1:894090/151‑1 (MQ=255)
                 gATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑c                                                                                                                                 >  2:589174/1‑151 (MQ=255)
                    ggggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑ccc                                                                                                                               <  2:426383/150‑1 (MQ=255)
                       gCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑Cccgccg                                                                                                                           >  1:188678/1‑151 (MQ=255)
                          aGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGATCTGAT‑CCCGCCGtt                                                                                                                         >  2:897371/1‑150 (MQ=255)
                               tttGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGGGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTa                                                                                                                    >  2:831388/1‑150 (MQ=255)
                               ttgggggaAAATTAAACTGGGCGAAAATGGAAC‑CCTGAATCCAGCGGCCAGCGTGGGGGAAATGGGGGCGTGCAGGGGGGTG‑ATAATGGCTACAATGCCCCCGCGGTGATGGGGGGCCTAAAATATAAGTTCTGAT‑CCCGCCGTTAGCt                                                                                                                     <  1:668176/146‑1 (MQ=255)
                                     ggAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAAcc                                                                                                             >  2:424499/1‑151 (MQ=255)
                                       aaaTTAATCTGGGGGAAAAGGGCAC‑CCTGACCCCAGCGCCCAGCGTCGGGGGAATGGGGGCGTGCAGCTGGGTA‑ATAATGGCAACAATCCCCCCGAGGTGATGGTGCGCTTAAAATATAAGTTTTGTC‑CCCGCCGTTAGCTAAAAAACcgc                                                                                                           <  1:581107/151‑1 (MQ=255)
                                       aaaTTAAACTGGGCGTAAAGGGAAC‑CCTGAACCCAGCGCCCAGCGTCGGGGGATTGTGGGCGTGCAGCGGGGTA‑ATAATGGCTACAAGCACCCCGCAGTGAGGGGGGGCCTAAAATATAATTTCTGAC‑CCCGCCGTTAGCTAAAAAACcg                                                                                                            <  1:828687/150‑1 (MQ=255)
                                        aaTAAAACTGGGCATAAATGGAAC‑CCGGAACCCAGCGGCCGGCGTCGGGGGAGTGGGGGCGTGCGGCTGGGTA‑ATAATGCCTACAATCCCCCCGAAGTGATGGTGGGCCTAAAATAAAAGTTCTGCT‑CCCGCCGTTAGCTAAAAAACcgcg                                                                                                          <  1:273240/151‑1 (MQ=255)
                                            aacctgggCGAAAAGGGCAC‑CCTGAACCCGGCGGCCAGCGTCGGGGGAAGGTGGGCGTGCAGGTGGGTG‑ATAATGGCTACAAGGACCCCGCAGTGATGGGGGGCCTAAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGcgtc                                                                                                        <  2:701158/143‑1 (MQ=255)
                                                  gggcgaaaaGGCAACCCTAAACCCAGCGCCCGGCGTCGGGGGAAGGGGGGGGGGCGGCGGGGTG‑ATAATGCCTAAAATGCCCCCGAAGTGAGGGGGGGCGTAAAATATAAGTTCTGCC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCa                                                                                                 <  1:589162/146‑1 (MQ=255)
                                                      tAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc                                                                                            >  2:271166/1‑151 (MQ=255)
                                                        aaTGGCAC‑CCTAAACCCGGCGGCCGGCGTGGGGGGAAGGGGGGCGGGAAGCGGGGTA‑ATAAGGGCTAAAATGCCCCCGCAGTGATGGGGGGCTAAAAATATAAGTTCTGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc                                                                                          <  1:910729/151‑1 (MQ=255)
                                                           ttccgatcTGAACCCAGCGGCCAGCGTGGGGGGAAGGGGGGGGTGCAGCGGGGTA‑ATAATGGCTAAAATGACCCCGGAGTGAGGGGGGGCGAAAAATATATGTTCTGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                               <  1:805726/137‑1 (MQ=255)
                                                           ggCAA‑CCTAAACCCAGCGGCCGGCGTCGGGGGAATGTGGGCGTGCAGCGGGGTA‑ATAATGGCTACAATGCCCCCGCAGTGAGGGTGGGCTTAAAATATAATTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgcg                                                                                         <  1:439146/149‑1 (MQ=255)
                                                            gCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGt                                                                                       >  2:384235/1‑150 (MQ=255)
                                                              aa‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                     >  1:626705/1‑150 (MQ=255)
                                                               acccTAAACCCGGCGGCCGGGGTCGGGGGAATGGGGGCGGGCGGCTGGGTA‑ATAAGGGCTAAAATACCACCGCAGTGGTGGGGGGCTTAAAATATAAGTTCGGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACa                                                                                    <  2:585517/149‑1 (MQ=255)
                                                               a‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                                   >  2:452256/1‑151 (MQ=255)
                                                                 cgtaaagccggcgcCCGGCGTGGGGGGATTGTGGGGGTGCAGCGGGGTA‑ATATTGGCTACAATGCCCCCGAAGGGAGGGGGGGCATAAAATAAATGTCCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata                                                                                  <  1:76734/141‑1 (MQ=255)
                                                                 ccTGACCCCAGCGCCCGGCGTCGGGGGAGGGTGGGCGTGCAGTGGGGTG‑ATAATGGCTAAAATGCCACCGCAGTGAGGGTGGCCCTAAAATATAATTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                                   <  2:27259/150‑1 (MQ=255)
                                                                    gAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                               >  2:805742/1‑135 (MQ=255)
                                                                     aacccggCGGCCAGCGTCGGGGGAATGGGGGCGTGCAGCGGGGTG‑AAAAGGGCTACAAGGACCCCGAAGTGAGGGGGGGCTAAAAATATAATTTCTGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat                                                                                 <  2:835174/142‑1 (MQ=255)
                                                                     aacccAGCGCCCAGCGTGGGGGGATTGGGGGCGTGCAGTTGGGTG‑ATAATGCCTACAATGCCCCCGCAGTGATGGTGGGCATAAAAAATAAGTTCTGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa                                                                              <  1:297949/148‑1 (MQ=255)
                                                                     aaTCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa                                                                               >  2:466838/1‑150 (MQ=255)
                                                                             gccgaGCGTCGGGGGAATGGGGGCGTGCAGTGGGGTA‑ATAAGGGCTAAAATGACACCGCAGTGATGGGGGGCCTAAAATATAAGTTATGAA‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                       <  2:462949/146‑1 (MQ=255)
                                                                                    gTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGcc                                                               >  1:367620/1‑151 (MQ=255)
                                                                                        ggggAATGTGGGGGGGCGGC‑GGGTGAATATTGTTTAAAAGTACACGGGAGTGAGGGTGGGCCGGAAAAATAATTTCTCAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGcc                                                               <  2:160744/147‑1 (MQ=255)
                                                                                         gggAATGGGGGGGTGCAGTGGGGTA‑ATATTGTCTACAATGACCCCGAAGTGTTGGGGGCCCAAAAATAAAGTTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAg                                                          <  1:815986/151‑1 (MQ=255)
                                                                                           aaaggtgGGGGTGCGGGGGGGTG‑AAAATGGCAAAAATGACCCCGAGGTGATGGGGGGCTTAAAATAAAGTTTCGGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGc                                                         <  1:745624/146‑1 (MQ=255)
                                                                                               gtgGGCGTGCAGCTGGGTA‑ATAATGGCTACAATGACCCCGCAGTGATGGTGGGCCTAAAATATAATTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACg                                                      <  2:755904/149‑1 (MQ=255)
                                                                                                     gTGCAGTGGGGTG‑ATAATGCCTAAAAGAACCCGGGAGGGTTGGGGGGCATAAAATATAA‑TTTTGATCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtgt                                              <  1:1301/151‑1 (MQ=255)
                                                                                                      tGCAGCTGGGTG‑ATAATGGCTACAATGACCCCGCAGTGATGGGGGGCCTAAAATATAATTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtgt                                              <  2:708398/150‑1 (MQ=255)
                                                                                                         aGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatc                                          >  2:711832/1‑151 (MQ=255)
                                                                                                                 g‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCg                                  >  2:795682/1‑151 (MQ=255)
                                                                                                                   aaaagggccaaaaagACCCCGCGGTGAGGGGGGGCCTAAAATAAAAGTTCGGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCaa                                        <  1:266501/130‑1 (MQ=255)
                                                                                                                     aaTGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCg                                  >  1:169453/1‑148 (MQ=255)
                                                                                                                       tGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATc                             >  2:386750/1‑151 (MQ=255)
                                                                                                                            aCAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAc                          >  1:490589/1‑149 (MQ=255)
                                                                                                                              aaTGACACCGCAGTGATGGAGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                       >  2:772667/1‑150 (MQ=255)
                                                                                                                                  acacCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTg                  >  2:478529/1‑151 (MQ=255)
                                                                                                                                   ccccGGAGTGGTGGGGGGCCTAAAATATATTTTCTGTC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAttt                   <  1:952424/147‑1 (MQ=255)
                                                                                                                                   ccccGCAGTGAGGGGGGGCAAAAAATATAATTTCTCAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAttt                   <  1:184581/147‑1 (MQ=255)
                                                                                                                                     ccGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg               >  2:493342/1‑151 (MQ=255)
                                                                                                                                      cgcggtgggggtggGCAAAAAATATAGTTTCGGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg               <  1:424488/141‑1 (MQ=255)
                                                                                                                                       gCAGTGTTGGTGGCCCTAAAATAAAATTCCTGAC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg               <  1:831372/149‑1 (MQ=255)
                                                                                                                                       gCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt             >  2:284/1‑151 (MQ=255)
                                                                                                                                           tGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCttt         >  2:891850/1‑151 (MQ=255)
                                                                                                                                               ggtggGCTTAAAAAAAAAGTTCCTAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAg       <  1:82445/149‑1 (MQ=255)
                                                                                                                                               gggggcccTAAAATATAAGTTTTGCC‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt           <  2:199660/139‑1 (MQ=255)
                                                                                                                                                   ggCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCaa  >  1:830253/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                 
GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTG‑ATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGAT‑CCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAA  >  NC_000913/1222083‑1222377

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: ‑