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

ATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGC  >  NC_000913/213193‑213489
                                                                                                                                                     |                                                                                                                                                       
aTGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAg                                                                                                                                                       <  1:597739/151‑1 (MQ=255)
      cACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTAc                                                                                                                                                  <  1:480128/150‑1 (MQ=255)
         ccTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCt                                                                                                                                              <  2:169846/151‑1 (MQ=255)
           tGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGAAggggga                                                                                                                                                  >  1:764375/1‑141 (MQ=255)
            gTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCtt                                                                                                                                             <  2:93066/149‑1 (MQ=255)
               tACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCGTGGCAGTGGAAACCTTAAtt                                                                                                                                         >  2:755470/1‑149 (MQ=255)
                aCGTTTGGGCGCGTTTGAAATCCGCCGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCCGGGGAgggggcaagtcta                                                                                                                                           >  1:39058/1‑135 (MQ=255)
                 cGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCGGGGGAGGGGTAAGCTTAATGCgg                                                                                                                                      >  2:740633/1‑151 (MQ=255)
                  gTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGAAAGCTTAAAGCggg                                                                                                                                     >  2:877563/1‑151 (MQ=255)
                  gTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGGGGGGACAGCTTaaatcgcg                                                                                                                                     >  1:229823/1‑145 (MQ=255)
                  gTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCggg                                                                                                                                     <  1:185736/151‑1 (MQ=255)
                   tttGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCggg                                                                                                                                     >  2:950323/1‑150 (MQ=255)
                    ttGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGa                                                                                                                                    <  2:973261/150‑1 (MQ=255)
                     tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGGACAGCTTTAACCGGGAgg                                                                                                                                  >  1:110986/1‑151 (MQ=255)
                      gggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAgg                                                                                                                                  <  1:588815/150‑1 (MQ=255)
                         cgcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGat                                                                                                                              <  2:432080/151‑1 (MQ=255)
                          gcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGCgggggggg                                                                                                                                >  1:234374/1‑146 (MQ=255)
                            gTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCGGGGGGGGTACAGCCTAATTGGGGGGGGGATAtt                                                                                                                           >  1:167208/1‑151 (MQ=255)
                             tttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGCGGTAtt                                                                                                                           >  1:366285/1‑150 (MQ=255)
                                   aTCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGATTAATACGGGGGGCGAGATTTGccccc                                                                                                                    >  1:249992/1‑149 (MQ=255)
                                           cTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGTAAGCCTAAATCGGGGGGCGATAtttcgcctcccgctg                                                                                                              >  1:34271/1‑136 (MQ=255)
                                            tATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCGGGGGAGGGGACAGCTTAATTCGGGAGGCGATATTTGCCCCCCGCGt                                                                                                              >  2:300885/1‑147 (MQ=255)
                                               cAGTCGCTACTGAGGTGGATTCAAACCGTCGCCGGGCAAAGAGAAATCATTGTGGCGTGGCAGTCGTGGCGTCCAGCGCAGTAATTACCTGGGGGGCTGGGAAGTGTACAACTTAATGCGGGGGGCGATATTCGCCCTCCGCGtacggcg                                                                                                        <  1:516300/150‑7 (MQ=255)
                                               cAGTCGCAACTGTGGTGGATTAAATCCGTCACCGAACAAAGCGAAAACAATGTGCCGTGGCAGACATGGCTTCAACCGCTGGAATTGCCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGAc                                                                                                        <  1:226423/150‑1 (MQ=255)
                                                aGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGAAAGCTCAATTCGGGGGGCGACATTCGccccc                                                                                                                    >  1:528825/1‑136 (MQ=255)
                                                 gTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGaa                                                                                                    <  2:568387/151‑1 (MQ=255)
                                                     cAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAATGCGGGGGGGGATATTCGCCCTCCGCGGGCAGACGAAGCgg                                                                                                >  2:944402/1‑151 (MQ=255)
                                                       aCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGAAGTGTACAGCTTAATTCGGGAGGGGATATTCGCCCTCCGCGTGCAGACGAAGCgg                                                                                                >  2:828652/1‑149 (MQ=255)
                                                          gtggtgGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAgc                                                                                           <  1:656877/151‑1 (MQ=255)
                                                                   tAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAAGGTACAGCTTAATTCGGGGGGGGATATTCGCCCCCCCCGTGGAGACGAAACGGGGAGCGGGCGgttt                                                                                 >  1:318427/1‑150 (MQ=255)
                                                                         cGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAgcg                                                                           <  2:511130/151‑1 (MQ=255)
                                                                          gTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGGGGTATTTGCCCCCCCCGTGCAGACGAAGGGGGCCGCGGGCGGTTTAAAgcgc                                                                          >  2:647152/1‑151 (MQ=255)
                                                                               cGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGCGGGGGGGGGTATTTGCCCTCCGCGGGCAGACGGAGCGGTCAGCGGGCGGTTCAAAACGCCAgg                                                                      >  2:281599/1‑150 (MQ=255)
                                                                                gggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGGAAAGCTTAATGCGGGGGGCGGTATTTTCCCTTCGCGTGCAGACGAAGCGGGCCGCGTGCGGTTCAAAGCGCCAgg                                                                      >  1:820843/1‑149 (MQ=255)
                                                                                gggCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGAt                                                                    <  2:556676/151‑1 (MQ=255)
                                                                                     aaGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATtgctg                                                               <  2:928749/151‑1 (MQ=255)
                                                                                        cGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTTGGGAGAGTACAGCGTAATGAGGGGGGGGGTATTCGCCCCCCGCCTGCCGACGAAGCGGGCAGCGGGGGGTTCTAAACGCCCgggat                                                                   >  1:323768/1‑141 (MQ=255)
                                                                                          aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGAGGGGAGGGGGCAGCTTAATATGGGGGGCGGTATTTCTCCCCCCCGGGTAGGCGCAGGGGGACGCGTGCGGTTTAAAACCCCCGGActgctg                                                               >  1:66743/1‑146 (MQ=255)
                                                                                          aaaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCatat                                                          <  2:688621/151‑1 (MQ=255)
                                                                                           aaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATAtt                                                         <  2:633046/151‑1 (MQ=255)
                                                                                                  gtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGCCGGGAGCGGAAAACGTAAAGTGGGGGGGGGCATTTTCCCCCCCCGCGGAACGCAAGGGGTCAGCGGGCCCTTCAAAACGCCAGGGCTGCTTCACATTGtttggg                                                  >  2:922086/1‑146 (MQ=255)
                                                                                                  gtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCggg                                                   <  2:181939/150‑1 (MQ=255)
                                                                                                    gCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCg                                                 <  2:593700/150‑1 (MQ=255)
                                                                                                     ccGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAGCTTAATGCGGGGGGGGATATTTGCCCTCCCCGGGCAGACGAAACGGGCAGCGGGCGGTTTAAAACGCCAGGGTTGCTGCATATTTTTggggg                                                 >  2:65403/1‑147 (MQ=255)
                                                                                                        tGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAATGCGGGAGGGGATATTTCCCCCCCCCGGGCAGACGAAGCGGGCAGCGGGCGTTTTCAAACGCCCGGATTGCTGCCAATTTTCGGGCGGAAcg                                            >  1:462270/1‑151 (MQ=255)
                                                                                                             gACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGGGGCGGTATTTGCCCCCCCCGTGGAGACGAAGCGGTCCGCGGGCGGTTCAAAAAGCCCGGGATGCTGCCAATTTTCGGGGGGAAcggcg                                         >  1:541470/1‑149 (MQ=255)
                                                                                                             gACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGa                                       <  2:449794/151‑1 (MQ=255)
                                                                                                                gTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGt                                    <  1:195670/151‑1 (MQ=255)
                                                                                                                 tGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGCCGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGt                                    <  2:637978/150‑1 (MQ=255)
                                                                                                                    cTTCAACCGCTGGAATTACCGGCGGGGGCGGGGAGTGGAAAGCTTAATTCGGGGGGCGCTATTTGCCCCCCCCCGGCAGACGAAAGGGGCAGCGGGGGGTTTCAAAAGCCCAGGTTGCTGCAATTTGTTGGGGGTTAaggcg                                         >  1:600367/1‑142 (MQ=255)
                                                                                                                      tCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGTCAGCCTTATGAGGGGGGGGTTATTCGCCCTCCCCGGTGCGAGGAAACGGGTAGCGGGGGGTTCAAAGCGCCCGGAATGCTGTATTTTTTCGGGGGGAAAggcgggcg                                     >  1:921379/1‑144 (MQ=255)
                                                                                                                       cAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTaa                             <  1:496497/151‑1 (MQ=255)
                                                                                                                        aaCCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGATTAATCGGGGAGGGGATAATCGCCCTCCGCGGGGGGACGAAGCGGGCAGGGGGGGGTTCAAAACGCCCGGATTGTTGCATATTTTCGGGGGGAACGGGGGGCgcg                                   >  2:204807/1‑142 (MQ=255)
                                                                                                                            gCTGGAATTACCGGCGGGGGTGGGGAGGGTTAAGCCTAATTCGGGGGGGGATATTTGGCCCCCCCGGGGAGACGAAGCGGGCAGGGTGCGGTTCTAAAAGCCAGGATTGGTGGATATTTTTGGGGGGGACGGGGGACGCTAGCTaaaaaa                         >  1:498455/1‑150 (MQ=255)
                                                                                                                            gCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGaaa                        <  2:658593/151‑1 (MQ=255)
                                                                                                                                         gcggGGCTGGGAAGTGCACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGc           <  2:777163/151‑1 (MQ=255)
                                                                                                                                            ggggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTg         <  1:243241/151‑1 (MQ=255)
                                                                                                                                                  gggAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTg   >  2:564363/1‑151 (MQ=255)
                                                                                                                                                  gggAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCAAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGca     <  1:820536/149‑2 (MQ=255)
                                                                                                                                                   ggAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGt    >  1:675984/1‑149 (MQ=255)
                                                                                                                                                   ggAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTgc  <  1:279564/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
ATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGCGTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGC  >  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: