Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913213,3410AG46.7% 68.5 / 22.2 61G337G (GGA→GGG) tilStRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (4/28);  new base G (28/0);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 3.55e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGC  >  NC_000913/213193‑213489
                                                                                                                                                     |                                                                                                                                                       
aTGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAg                                                                                                                                                       <  1:597739/151‑1 (MQ=255)
      cACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTAc                                                                                                                                                  <  1:480128/150‑1 (MQ=255)
         ccTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCt                                                                                                                                              <  2:169846/151‑1 (MQ=255)
           tGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGAAggggga                                                                                                                                                  >  1:764375/1‑141 (MQ=255)
            gTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCtt                                                                                                                                             <  2:93066/149‑1 (MQ=255)
               tACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCGTGGCAGTGGAAACCTTAAtt                                                                                                                                         >  2:755470/1‑149 (MQ=255)
                aCGTTTGGGCGCGTTTGAAATCCGCCGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCCGGGGAgggggcaagtcta                                                                                                                                           >  1:39058/1‑135 (MQ=255)
                 cGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCGGGGGAGGGGTAAGCTTAATGCgg                                                                                                                                      >  2:740633/1‑151 (MQ=255)
                  gTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGGAGGGGAAAGCTTAAAGCggg                                                                                                                                     >  2:877563/1‑151 (MQ=255)
                  gTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGGGGGGACAGCTTaaatcgcg                                                                                                                                     >  1:229823/1‑145 (MQ=255)
                  gTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCggg                                                                                                                                     <  1:185736/151‑1 (MQ=255)
                   tttGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGTGTACAGCTTAATGCggg                                                                                                                                     >  2:950323/1‑150 (MQ=255)
                    ttGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGa                                                                                                                                    <  2:973261/150‑1 (MQ=255)
                     tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGGAGTGGACAGCTTTAACCGGGAgg                                                                                                                                  >  1:110986/1‑151 (MQ=255)
                      gggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAgg                                                                                                                                  <  1:588815/150‑1 (MQ=255)
                         cgcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGat                                                                                                                              <  2:432080/151‑1 (MQ=255)
                          gcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGGGGACAGCTTAATGCgggggggg                                                                                                                                >  1:234374/1‑146 (MQ=255)
                            gTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGCGGGGGGGGTACAGCCTAATTGGGGGGGGGATAtt                                                                                                                           >  1:167208/1‑151 (MQ=255)
                             tttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGCGGTAtt                                                                                                                           >  1:366285/1‑150 (MQ=255)
                                   aTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGGAGGGTACAGATTAATACGGGGGGCGAGATTTGccccc                                                                                                                    >  1:249992/1‑149 (MQ=255)
                                           cTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGGAGGGGTAAGCCTAAATCGGGGGGCGATAtttcgcctcccgctg                                                                                                              >  1:34271/1‑136 (MQ=255)
                                            tATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCGGGGGAGGGGACAGCTTAATTCGGGAGGCGATATTTGCCCCCCGC‑Gt                                                                                                              >  2:300885/1‑147 (MQ=255)
                                               cAGTCGCTACTGAGGTGGATTCAAACCGTCGCCGGGCAAAGAGAAATCATTGTGGCGTGGCAGTCGTGGCGTCCAGCGCAGTAATTACCTGG‑GGGGCTGGGAAGTGTACAACTTAATGCGGGGGGCGATATTCGCCCTCCGC‑Gtacggcg                                                                                                        <  1:516300/150‑7 (MQ=255)
                                               cAGTCGCAACTGTGGTGGATTAAATCCGTCACCGAACAAAGCGAAAACAATGTGCCGTGGCAGACATGGCTTCAACCGCTGGAATTGCCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGAc                                                                                                        <  1:226423/150‑1 (MQ=255)
                                                aGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGGAGGGGAAAGCTCAATTCGGGGGGCGACATTCGccccc                                                                                                                    >  1:528825/1‑136 (MQ=255)
                                                 gTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑aa                                                                                                    <  2:568387/151‑1 (MQ=255)
                                                     cAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGGAGGGGACAGCTTAATGCGGGGGGGGATATTCGCCCTCCGC‑GGGCAGACG‑AAGCgg                                                                                                >  2:944402/1‑151 (MQ=255)
                                                       aCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGAAGTGTACAGCTTAATTCGGGAGGGGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCgg                                                                                                >  2:828652/1‑149 (MQ=255)
                                                          gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAgc                                                                                           <  1:656877/151‑1 (MQ=255)
                                                                   tAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAAGGTACAGCTTAATTCGGGGGGGGATATTCGCCCCCCCC‑GTGGAGACG‑AAACGGGGAGC‑GGGCGgttt                                                                                 >  1:318427/1‑150 (MQ=255)
                                                                         cGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAgcg                                                                           <  2:511130/151‑1 (MQ=255)
                                                                          gTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGGGGTATTTGCCCCCCCC‑GTGCAGACG‑AAGGGGGCCGC‑GGGCGGTTTAAAgcgc                                                                          >  2:647152/1‑151 (MQ=255)
                                                                               cGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGGGGACAGCTTAATGCGGGGGGGGGTATTTGCCCTCCGC‑GGGCAGACG‑GAGCGGTCAGC‑GGGCGGTTCAAAACGCCAgg                                                                      >  2:281599/1‑150 (MQ=255)
                                                                                gggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGGAGTGGAAAGCTTAATGCGGGGGGCGGTATTTTCCCTTCGC‑GTGCAGACG‑AAGCGGGCCGC‑GTGCGGTTCAAAGCGCCAgg                                                                      >  1:820843/1‑149 (MQ=255)
                                                                                gggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGAt                                                                    <  2:556676/151‑1 (MQ=255)
                                                                                     aaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATtgctg                                                               <  2:928749/151‑1 (MQ=255)
                                                                                        cGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTTGGGAGAGTACAGCGTAATGAGGGGGGGGGTATTCGCCCCCCGC‑CTGCCGACG‑AAGCGGGCAGC‑GGGGGGTTCTAAACGCCCgggat                                                                   >  1:323768/1‑141 (MQ=255)
                                                                                          aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGAGGGGAGGGGGCAGCTTAATATGGGGGGCGGTATTTCTCCCCCCC‑GGGTAGGCG‑CAGGGGGACGC‑GTGCGGTTTAAAACCCCCGGActgctg                                                               >  1:66743/1‑146 (MQ=255)
                                                                                          aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCatat                                                          <  2:688621/151‑1 (MQ=255)
                                                                                           aaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATAtt                                                         <  2:633046/151‑1 (MQ=255)
                                                                                                  gtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGCCGGGAGCGGAAAACGTAAAGTGGGGGGGGGCATTTTCCCCCCCC‑GCGGA‑ACGCAAGGGGTCAGC‑GGGCCCTTCAAAACGCCAGGGCTGCTTCACATTGtttggg                                                  >  2:922086/1‑146 (MQ=255)
                                                                                                  gtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCggg                                                   <  2:181939/150‑1 (MQ=255)
                                                                                                    gCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCg                                                 <  2:593700/150‑1 (MQ=255)
                                                                                                     ccGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGTGGACAGCTTAATGCGGGGGGGGATATTTGCCCTCCCC‑GGGCAGACG‑AAACGGGCAGC‑GGGCGGTTTAAAACGCCAGGGTTGCTGCATATTTTTggggg                                                 >  2:65403/1‑147 (MQ=255)
                                                                                                        tGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGGTGGGGAGGGGACAGCTTAATGCGGGAGGGGATATTTCCCCCCCCC‑GGGCAGACG‑AAGCGGGCAGC‑GGGCGTTTTCAAACGCCCGGATTGCTGCCAATTTTCGGGCGGAAcg                                            >  1:462270/1‑151 (MQ=255)
                                                                                                             gACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGCGGTATTTGCCCCCCCC‑GTGGAGACG‑AAGCGGTCCGC‑GGGCGGTTCAAAAAGCCCGGGATGCTGCCAATTTTCGGGGGGAAcggcg                                         >  1:541470/1‑149 (MQ=255)
                                                                                                             gACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGa                                       <  2:449794/151‑1 (MQ=255)
                                                                                                                gTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGt                                    <  1:195670/151‑1 (MQ=255)
                                                                                                                 tGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGCCG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGt                                    <  2:637978/150‑1 (MQ=255)
                                                                                                                    cTTCAACCGCTGGAATTACCGGC‑GGGGGCGGGGAGTGGAAAGCTTAATTCGGGGGGCGCTATTTGCCCCCCCC‑CGGCAGACG‑AAAGGGGCAGCGGGGGGTTTCAAAAGCCCAGG‑TTGCTGCAATTTGTTGGGGGTTAaggcg                                         >  1:600367/1‑142 (MQ=255)
                                                                                                                      tCAACCGCTGGAATTACCGGC‑GGGGCTGGGGAGTGGTCAGCCTTATGAGGGGGGGGTTATTCGCCCTCCCCGGTGC‑GAGG‑AAACGGGTAGC‑GGGGGGTTCAAAGCGCCCGGAATGCTGTATTTTTTCGGGGGGAAAggcgggcg                                     >  1:921379/1‑144 (MQ=255)
                                                                                                                       cAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTaa                             <  1:496497/151‑1 (MQ=255)
                                                                                                                        aaCCGCTGGAATTACCGGC‑GGGGGTGGGGAGGGTACAGATTAATCGGGGAGGGGATAATCGCCCTCCGC‑GGGGGGACG‑AAGCGGGCAGG‑GGGGGGTTCAAAACGCCCGGATTGTTGCATATTTTCGGGGGGAACGGGGGGCgcg                                   >  2:204807/1‑142 (MQ=255)
                                                                                                                            gCTGGAATTACCGGC‑GGGGGTGGGGAGGGTTAAGCCTAATTCGGGGGGGGATATTTGGCCCCCCC‑GGGGAGACG‑AAGCGGGCAGG‑GTGCGGTTCTAAAAGCCAGGATTGGTGGATATTTTTGGGGGGGACGGGGGACGCTAGCTaaaaaa                         >  1:498455/1‑150 (MQ=255)
                                                                                                                            gCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGaaa                        <  2:658593/151‑1 (MQ=255)
                                                                                                                                         gc‑ggGGCTGGGAAGTGCACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGc           <  2:777163/151‑1 (MQ=255)
                                                                                                                                            ggggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTg         <  1:243241/151‑1 (MQ=255)
                                                                                                                                                  gggAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTg   >  2:564363/1‑151 (MQ=255)
                                                                                                                                                  gggAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCAAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGca     <  1:820536/149‑2 (MQ=255)
                                                                                                                                                   ggAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGt    >  1:675984/1‑149 (MQ=255)
                                                                                                                                                   ggAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTgc  <  1:279564/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
ATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGC‑GGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGC‑GTGCAGACG‑AAGCGGTCAGC‑GTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGC  >  NC_000913/213193‑213489

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