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→GGGtilStRNA(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

TGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTG  >  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)
                                                                  tAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGttt                                                                           <  2:67686/150‑1 (MQ=255)
                                                                     aTCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCCCGTGCAGACGAAACGGGCAGCGGGCGGTTTaa                                                                        >  1:252134/1‑150 (MQ=255)
                                                                         gTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTTAAGCGGGGGGGGAGATTTGCCCCCCCCGGGCCGACGAAGCGGGCAGCGTGCGGTTCTAAacgc                                                                   >  2:785415/1‑151 (MQ=255)
                                                                             ccGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCCCGTGCAGACGAAGCGGGGAGCGGGCGGTTTAAAGCGCCAgg                                                               >  1:527127/1‑151 (MQ=255)
                                                                               gggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTTATGCGGGGGGCGGTATTCGCCCTCCCCGTGCAGACGAAGCGGGCAGCGGGCGTTTCAAAACGCCAgg                                                               >  1:621304/1‑149 (MQ=255)
                                                                               gggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGAt                                                             <  1:248441/151‑1 (MQ=255)
                                                                                ggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGAtt                                                            <  1:885943/151‑1 (MQ=255)
                                                                                    aaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCATCCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATGCGGGAGGGGGTATTCGCCCCCCGCGGGCAGACGAAGCGGGCGGCGGGGGTTTTTAAACGCCAGGAtt                                                            >  1:441778/1‑147 (MQ=255)
                                                                                     aGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGCGATATTTGCCCTCCCCGTGCAGACGAAGCGGTCAGCGGGCGGTTCAAAACGCCAGGGTtgctg                                                        >  1:720368/1‑150 (MQ=255)
                                                                                       cGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATTCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGGGCGTTTCAAAGCGCCAGGATTGCTGCat                                                     >  1:139443/1‑151 (MQ=255)
                                                                                       cGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCCTAATGCGGGGGGGGGTATTTGCCCCCCCCGGGCAGACGGAAGGGGCAGCGGGGGGTTTAAAACGCCAGGATTGCTGCa                                                      >  1:325589/1‑150 (MQ=255)
                                                                                       cGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAAATCGGGAGGGGGTATTCCCCCTCCGCGTGCAGACGAAAGGGGCAGCGGGGGGTTTAAAGCGCCCGGGTtgct                                                         >  2:387792/1‑147 (MQ=255)
                                                                                        gAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGAAGGGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCCCGGGCAGACGAAGCGGGCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCat                                                     >  1:98744/1‑150 (MQ=255)
                                                                                         aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGCGGGGGGCGGTATTCGCCCCCCGCGTGCAGACGAAGCGGGCAGCGGGCGGTTTCAAACGCCCGGGTtgctg                                                        >  1:283814/1‑146 (MQ=255)
                                                                                             cATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATGCGGGGGGGGGTATTTGCCCTCCCCGGGCAGACGAAGCGGGCAGCGGGCGGTTTCAAGCGCCCGGGTTGCTGCATATTGt                                                >  1:746043/1‑150 (MQ=255)
                                                                                                            gACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGGTTTATTCGGGAGGGGATATTTCCCCTCCGCGTGCAGACGAAGCGGGCAGCGTGCGGTTCCAAAAGCCAGGAATGTCGCATATTGTTGGGCGGAAcggcggg                                >  1:30592/1‑150 (MQ=255)
                                                                                                            gACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGa                                <  1:331785/151‑1 (MQ=255)
                                                                                                             aCGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGGGTGTACAACTTAATGCGGGGGGGGATATTCGCCCTCCGCGGGCAGACGAAAAGGGGAGCGTGCGCTTTAAAACGCCCGGATTGCTGCATATTTTCGGGGGGAAcggcggg                                >  2:840639/1‑149 (MQ=255)
                                                                                                             aCGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGAc                               <  2:958589/151‑1 (MQ=255)
                                                                                                                 ggCTTCAACCGCTGGAAGTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGGGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTaa                           <  2:674669/151‑1 (MQ=255)
                                                                                                                  gCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAg                          <  2:728755/151‑1 (MQ=255)
                                                                                                                     tCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCt                        <  1:370819/150‑1 (MQ=255)
                                                                                                                          cGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGaa                  <  1:64153/151‑1 (MQ=255)
                                                                                                                           gCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGaaa                 <  2:105049/151‑1 (MQ=255)
                                                                                                                             tGGAATTACCGGCGGGGGTGGGGAGTGAACAACTTCAAGCGGGGGGCGAAAATCCCCCTCCGCGTGCAGACGAAGCGGGCCGCGTTCGGTTCCAAGGACCCGGGCTGCTGCAAATTGTCGGGGGGAACGGCGGGCCGAAGCTCAACaaaa                >  2:152440/1‑150 (MQ=255)
                                                                                                                              ggAATTACCGGCGGGGCTGGGGAGCGTACAGCCTACTGAGGGGGGGGTTATTCGCCCCCCGCGGGCAGACGAAGCGGGCAGGGGGGGGTTCCAAGCGGCCGGAATGGTGCATATTGGCGGGGGGAACGGCGGGCGTAAGCTaaaaaaaaac              >  2:341912/1‑149 (MQ=255)
                                                                                                                                   tACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGc          <  1:485028/150‑1 (MQ=255)
                                                                                                                                    aCCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGTGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCa         <  2:233821/150‑1 (MQ=255)
                                                                                                                                    aCCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTATGCTAAAGAAAATCTGGCa         <  1:589026/150‑1 (MQ=255)
                                                                                                                                    aCCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGCGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCaa        <  2:735884/151‑1 (MQ=255)
                                                                                                                                          ggggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTg  >  1:179263/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
TGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTG  >  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: