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

TGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTG  >  NC_000913/213194‑213482
                                                                                                                                                   |                                                                                                                                              
tGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAgg                                                                                                                                             >  1:959110/1‑150 (MQ=255)
  ccTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAgtgt                                                                                                                                           <  2:680255/151‑1 (MQ=255)
   cTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAggg                                                                                                                                            >  1:105087/1‑147 (MQ=255)
     cACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACa                                                                                                                                        <  1:734664/151‑1 (MQ=255)
        ccTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCGGGGGAggg                                                                                                                                            >  2:905961/1‑142 (MQ=255)
           gTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCCGGGGGGAGGACAGCttt                                                                                                                                   >  2:10030/1‑149 (MQ=255)
             ttACGTTTGGGCGTGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCCGGGGAgggggcagatt                                                                                                                                    >  2:132216/1‑137 (MQ=255)
             ttACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAAt                                                                                                                                 <  1:676451/150‑1 (MQ=255)
             ttACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGGGTaaaagttaa                                                                                                                                  >  1:453553/1‑141 (MQ=255)
               aCGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAATGc                                                                                                                               >  2:855006/1‑150 (MQ=255)
                   ttGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTAAAGCTTAATTCggggg                                                                                                                          >  2:886441/1‑149 (MQ=255)
                    tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAg                                                                                                                          <  2:303436/150‑1 (MQ=255)
                    tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCAGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAg                                                                                                                          <  2:286824/150‑1 (MQ=255)
                     ggccgcgTTTGAACTCCGACGCTATCAGTCGCAACTGTGGTGGACTAAATCCGTCACCGGGAAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCATCCGGTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGc                                                                                                                        <  2:275140/148‑1 (MQ=255)
                       gcgcGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGGGTGTACAGCGTTATACgggggg                                                                                                                         >  2:735232/1‑148 (MQ=255)
                         gcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGata                                                                                                                    <  2:576315/151‑1 (MQ=255)
                          cgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAAGGTACAGATTAATGCgggggggggta                                                                                                                    >  1:626106/1‑145 (MQ=255)
                           gTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATAtt                                                                                                                  <  2:831689/151‑1 (MQ=255)
                             ttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTc                                                                                                                 <  1:71969/150‑1 (MQ=255)
                               gAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACACCTTAATTCGGGAggaggcattcccc                                                                                                              >  1:831788/1‑140 (MQ=255)
                                   tCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGAGTAAAGCTTAAAGCGGGAGGGGATATTTGCCCTcc                                                                                                          >  1:781451/1‑151 (MQ=255)
                                    ccGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCTCTcc                                                                                                          <  2:582859/150‑1 (MQ=255)
                                     cGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTTAAGCGGGGGGGGATATTTGcccccccc                                                                                                        >  1:92575/1‑149 (MQ=255)
                                     cGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGGGTACAGCCTAATGCGGGAGGCGGTATTCGcccccccc                                                                                                        >  2:913136/1‑149 (MQ=255)
                                            aTCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAg                                                                                                  <  2:202527/150‑1 (MQ=255)
                                                gTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGaa                                                                                             <  1:861006/151‑1 (MQ=255)
                                                gTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGaa                                                                                             <  2:857695/151‑1 (MQ=255)
                                                      aCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTc                                                                                       <  2:631/151‑1 (MQ=255)
                                                         gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAACCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCGCCGCGTGCAGACGAAGCGGTCa                                                                                      <  1:29636/149‑1 (MQ=255)
                                                         gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCa                                                                                      <  1:479374/149‑1 (MQ=255)
                                                            gtggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgcgtg                                                                                 <  2:146281/151‑1 (MQ=255)
                                                                  tAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑Gttt                                                                           <  2:67686/150‑1 (MQ=255)
                                                                     aTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCCCGTGCAGACGAAACGGGCAGCGGGC‑GGTTTaa                                                                        >  1:252134/1‑150 (MQ=255)
                                                                         gTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTTAAGCGGGGGGGGAGATTTGCCCCCCCCGGGCCGACGAAGCGGGCAGCGTGC‑GGTTCTAAacgc                                                                   >  2:785415/1‑151 (MQ=255)
                                                                             ccGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCCCGTGCAGACGAAGCGGGGAGCGGGC‑GGTTTAAAGCGCCAgg                                                               >  1:527127/1‑151 (MQ=255)
                                                                               gggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTTATGCGGGGGGCGGTATTCGCCCTCCCCGTGCAGACGAAGCGGGCAGCGGGC‑GTTTCAAAACGCCAgg                                                               >  1:621304/1‑149 (MQ=255)
                                                                               gggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGAt                                                             <  1:248441/151‑1 (MQ=255)
                                                                                ggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGAtt                                                            <  1:885943/151‑1 (MQ=255)
                                                                                    aaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCATCCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATGCGGGAGGGGGTATTCGCCCCCCGCGGGCAGACGAAGCGGGCGGCGGGG‑GTTTTTAAACGCCAGGAtt                                                            >  1:441778/1‑147 (MQ=255)
                                                                                     aGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGCGATATTTGCCCTCCCCGTGCAGACGAAGCGGTCAGCGGGC‑GGTTCAAAACGCCAGGGTtgctg                                                        >  1:720368/1‑150 (MQ=255)
                                                                                       cGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATTCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGGGC‑GTTTCAAAGCGCCAGGATTGCTGCat                                                     >  1:139443/1‑151 (MQ=255)
                                                                                       cGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCCTAATGCGGGGGGGGGTATTTGCCCCCCCCGGGCAGACGGAAGGGGCAGCGGGG‑GGTTTAAAACGCCAGGATTGCTGCa                                                      >  1:325589/1‑150 (MQ=255)
                                                                                       cGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAAATCGGGAGGGGGTATTCCCCCTCCGCGTGCAGACGAAAGGGGCAGCGGGG‑GGTTTAAAGCGCCCGGGTtgct                                                         >  2:387792/1‑147 (MQ=255)
                                                                                        gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGAAGGGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCCCGGGCAGACGAAGCGGGCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCat                                                     >  1:98744/1‑150 (MQ=255)
                                                                                         aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGCGGGGGGCGGTATTCGCCCCCCGCGTGCAGACGAAGCGGGCAGCGGGCGGTTTCAAA‑CGCCCGGGTtgctg                                                        >  1:283814/1‑146 (MQ=255)
                                                                                             cATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATGCGGGGGGGGGTATTTGCCCTCCCCGGGCAGACGAAGCGGGCAGCGGGC‑GGTTTCAAGCGCCCGGGTTGCTGCATATTGt                                                >  1:746043/1‑150 (MQ=255)
                                                                                                            gACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGGTTTATTCGGGAGGGGATATTTCCCCTCCGCGTGCAGACGAAGCGGGCAGCGTGC‑GGTTCCAAAAGCCAGGAATGTCGCATATTGTTGGGCGGAAcggcggg                                >  1:30592/1‑150 (MQ=255)
                                                                                                            gACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGa                                <  1:331785/151‑1 (MQ=255)
                                                                                                             aCGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGGGTGTACAACTTAATGCGGGGGGGGATATTCGCCCTCCGCGGGCAGACGAAAAGGGGAGCGTGC‑GCTTTAAAACGCCCGGATTGCTGCATATTTTCGGGGGGAAcggcggg                                >  2:840639/1‑149 (MQ=255)
                                                                                                             aCGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGAc                               <  2:958589/151‑1 (MQ=255)
                                                                                                                 ggCTTCAACCGCTGGAAGTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGGGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTaa                           <  2:674669/151‑1 (MQ=255)
                                                                                                                  gCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAg                          <  2:728755/151‑1 (MQ=255)
                                                                                                                     tCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCt                        <  1:370819/150‑1 (MQ=255)
                                                                                                                          cGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGaa                  <  1:64153/151‑1 (MQ=255)
                                                                                                                           gCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGaaa                 <  2:105049/151‑1 (MQ=255)
                                                                                                                             tGGAATTACCGGCGGGGGTGGGGAGTGAACAACTTCAAGCGGGGGGCGAAAATCCCCCTCCGCGTGCAGACGAAGCGGGCCGCGTTC‑GGTTCCAAGGACCCGGGCTGCTGCAAATTGTCGGGGGGAACGGCGGGCCGAAGCTCAACaaaa                >  2:152440/1‑150 (MQ=255)
                                                                                                                              ggAATTACCGGCGGGGCTGGGGAGCGTACAGCCTACTGAGGGGGGGGTTATTCGCCCCCCGCGGGCAGACGAAGCGGGCAGGGGGG‑GGTTCCAAGCGGCCGGAATGGTGCATATTGGCGGGGGGAACGGCGGGCGTAAGCTaaaaaaaaac              >  2:341912/1‑149 (MQ=255)
                                                                                                                                   tACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGc          <  1:485028/150‑1 (MQ=255)
                                                                                                                                    aCCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGTGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCa         <  2:233821/150‑1 (MQ=255)
                                                                                                                                    aCCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTATGCTAAAGAAAATCTGGCa         <  1:589026/150‑1 (MQ=255)
                                                                                                                                    aCCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGCGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCaa        <  2:735884/151‑1 (MQ=255)
                                                                                                                                          ggggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTg  >  1:179263/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTG  >  NC_000913/213194‑213482

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