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

GATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGCCG  >  NC_000913/213192‑213491
                                                                                                                                                     |                                                                                                                                                       
gATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGaa                                                                                                                                                        <  1:550325/151‑1 (MQ=255)
    ccTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGAGCCGGGGAggg                                                                                                                                                     >  2:256/1‑148 (MQ=255)
     cTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAggg                                                                                                                                                     >  1:1014226/1‑147 (MQ=255)
        aCCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATGGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACa                                                                                                                                                 <  2:875087/150‑1 (MQ=255)
          ccTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTAc                                                                                                                                                  >  2:243798/1‑147 (MQ=255)
             gTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTaa                                                                                                                                           <  1:1215639/151‑1 (MQ=255)
             gTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTaa                                                                                                                                           <  1:6786/151‑1 (MQ=255)
              tttACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTaa                                                                                                                                           >  1:973400/1‑150 (MQ=255)
                  cGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTGCAGCTTaa                                                                                                                                           >  1:851225/1‑146 (MQ=255)
                     ttGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGa                                                                                                                                    <  2:369607/150‑1 (MQ=255)
                      tGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGCgggggg                                                                                                                                  >  1:290903/1‑151 (MQ=255)
                       gggCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAgg                                                                                                                                  <  2:414074/150‑1 (MQ=255)
                          cgcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGCg                                                                                                                                >  1:681848/1‑149 (MQ=255)
                          cgcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGat                                                                                                                              <  2:806218/151‑1 (MQ=255)
                          cgcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGat                                                                                                                              <  1:255354/151‑1 (MQ=255)
                           gcgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTTATTCgggggggg                                                                                                                                >  2:555147/1‑146 (MQ=255)
                            cgTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGatat                                                                                                                            <  2:833090/151‑1 (MQ=255)
                             gTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCCTAATGCTGGAGGCGgtaat                                                                                                                           >  1:1150529/1‑149 (MQ=255)
                               ttGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTc                                                                                                                          <  2:519586/150‑1 (MQ=255)
                                     tCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGAACGCTTAATGCGGGGGGCGGTATTCGcccccc                                                                                                                   >  1:864700/1‑151 (MQ=255)
                                      ccGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTcc                                                                                                                   <  1:595013/150‑1 (MQ=255)
                                      ccGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTcc                                                                                                                   <  2:1135825/150‑1 (MQ=255)
                                         aCGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGGGGGGCGAGCTTAAAGCGGGAGGCGGGATTTCCCCTcc                                                                                                                   >  2:447855/1‑147 (MQ=255)
                                         aCGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCGGGGGAGTGGAGAACTTAATGCGGGGGGGGGTATTCGCCCTCcgc                                                                                                                 >  2:526854/1‑149 (MQ=255)
                                          cGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTg                                                                                                              <  1:1010977/151‑1 (MQ=255)
                                           gCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGc                                                                                                             <  1:1160897/151‑1 (MQ=255)
                                                 aGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAACTTAATGCGGGGGGCGATTTTCGCCCTCCGCGGGCAg                                                                                                           >  2:206692/1‑147 (MQ=255)
                                                   tCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGGACAACTTTAATCGGGGGGGGGTAttttcccccc                                                                                                                   >  2:797622/1‑129 (MQ=255)
                                                    cGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGGGGGGGTATTCGCCCTCCGCGTGCAgaagaa                                                                                                      >  2:392036/1‑149 (MQ=255)
                                                      cAACTGTGGTGGATTAAATCCGTCCCCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCg                                                                                                   <  2:587424/150‑1 (MQ=255)
                                                       aaCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGTACAGCTTAATGCGGGGGGGGATATTTGCCCCCCCCGTGCAGACGAAACgg                                                                                                  >  2:248446/1‑150 (MQ=255)
                                                        aCTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTc                                                                                                <  1:1091276/151‑1 (MQ=255)
                                                            tggtggATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGAAAGCCTAAAGCGGGAGGGGGTATTTGCCCCCCCCGTGGAGACGAAACGGGCAgc                                                                                             >  2:847926/1‑150 (MQ=255)
                                                                         ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAgc                                                                              <  1:332532/151‑1 (MQ=255)
                                                                         ccGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAgc                                                                              <  2:124868/151‑1 (MQ=255)
                                                                             cACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAATTCGGGGGGCGGTATTCGCCCCCCCCGTGCCGGCGAAGGGGGCAGCGGGG‑GTTTCAAAACGcc                                                                           >  2:926180/1‑150 (MQ=255)
                                                                             cACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCa                                                                          <  1:1162929/151‑1 (MQ=255)
                                                                              aCCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCTTAATAAGGGGGGGGGTATTTGCCCTCCCCCGGCAGACGAAGCGGTTAGCGGGC‑GGTTCAAAACGcc                                                                           >  2:737560/1‑149 (MQ=255)
                                                                                         cGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAGCTTAATGCGGGAGGCGGTATTCGCCCCCCGCGCGCAGACGAAGCGGGCAGCGGGC‑GGTTCAAAGCGCCCGGATtgct                                                                  >  2:245724/1‑147 (MQ=255)
                                                                                            aaaCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATAtt                                                           <  1:148836/151‑1 (MQ=255)
                                                                                              aCATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGAAAGCTTTATGGGGGGGGGGATATTTGCCCCCCGCGTGGCGGCCAAACGGGCAGCGGGCGGTTTCAAA‑CGCCAGGGTTGGTGCcta                                                             >  2:833082/1‑147 (MQ=255)
                                                                                                aTTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGTGTACAGCTTAATTCGGGAGGGGATATTTGCCCTCCCCGTGCAGACGAAACGGGCAGCGTGC‑GGTTCAAAGCGCCAGGATTGCTGCATATTTTCg                                                       >  1:257023/1‑151 (MQ=255)
                                                                                                 ttGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGGTGGGGAGGGGACAGCCTAATGCGGGGGGCGATATTTGCCCCCCCCGGGGAGACGAAGCGGGCAGGGGGCGGTTTCAAA‑CGCCAGGAGTGCTGCATATTGTCg                                                       >  2:498774/1‑150 (MQ=255)
                                                                                                  tgtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGTGTACAACTTAATGCGGGAGGCGGTATTCGCCCTCCCCCGGCAGACGAAGCGGGCAGCGTGG‑GGTTCAAAACGCCAGGATTGCTGCATATTTTCggg                                                     >  1:838193/1‑151 (MQ=255)
                                                                                                  tgtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCgg                                                      <  1:56350/150‑1 (MQ=255)
                                                                                                  tgtgCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCgg                                                      <  2:290854/150‑1 (MQ=255)
                                                                                                     gCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCg                                                   <  1:825803/150‑1 (MQ=255)
                                                                                                     gCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGt                                                  <  1:800664/151‑1 (MQ=255)
                                                                                                       cGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGGACAGCTTAATGCGGGGGGCGATATTCGCCCTCCGCGTGCCGACGAAGAGGTCAGCCGGCGGTTTCAAA‑CGCCAGGGTTGCTGCATATTGTCggggg                                                   >  1:293881/1‑146 (MQ=255)
                                                                                                        gTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAc                                               <  2:1016806/151‑1 (MQ=255)
                                                                                                          ggCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAcgg                                             <  1:181483/151‑1 (MQ=255)
                                                                                                           gCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTACAGCTTAATGCGGGAGGCGATATTTCCCCCCCGCGTGCAGACGGAGCGGGCAGCGTGC‑GGTTTAAAGCGCCAGGAGTGCTGCCTATTTTCCGGGGTAAagg                                             >  1:156129/1‑150 (MQ=255)
                                                                                                             aGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAAcggcgg                                          <  2:159942/151‑1 (MQ=255)
                                                                                                              gACGTGGCTTCAACCGCTGGAATTACCGTCGAGGCTGAGAAGTGTACATATTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGTCAGGATTGCAGCATATTGTCGGGCGGAACGGCGGa                                         <  2:342808/151‑1 (MQ=255)
                                                                                                               aCGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGAc                                        <  2:765573/151‑1 (MQ=255)
                                                                                                                    gCTTCAACCGCTGGAATTACCGGCGGGGCTGGGGAGGGTAAAGCCTAATGCGGGGGGGGATATTCCCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCCGGATTGCTGCATATTGTCGGGCGTAACGGCGGGCGTa                                     >  1:326637/1‑149 (MQ=255)
                                                                                                                     cTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGGGGACAGCTTAAAGCGGGGGGGGGTATTTGCCCTCCGCGCGCAGACGAAGCGGGCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTTTCGGGCGTAACCGGGGACGGaa                                    >  1:661939/1‑149 (MQ=255)
                                                                                                                        cAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTaa                               <  2:918132/151‑1 (MQ=255)
                                                                                                                         aaCCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTaaa                              <  2:1215148/151‑1 (MQ=255)
                                                                                                                                      aCCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCa                  <  2:1150343/150‑1 (MQ=255)
                                                                                                                                         ggcggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAag                >  2:1075276/1‑149 (MQ=255)
                                                                                                                                           cggGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCt            <  1:673987/151‑1 (MQ=255)
                                                                                                                                             gggCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTgg          <  2:1075143/151‑1 (MQ=255)
                                                                                                                                               gCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTggg         >  2:1182355/1‑150 (MQ=255)
                                                                                                                                                cTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCg       >  1:75928/1‑151 (MQ=255)
                                                                                                                                                 tGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCg       >  2:1180255/1‑150 (MQ=255)
                                                                                                                                                  gggAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGt      <  1:1018420/150‑1 (MQ=255)
                                                                                                                                                     aaGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTgccg  <  1:788572/151‑1 (MQ=255)
                                                                                                                                                     aaGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTgccg  >  1:113486/1‑151 (MQ=255)
                                                                                                                                                     aaGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTgcc   >  1:656312/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
GATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCAACCGCTGGAATTACCGGCGGGGCTGGGAAGTGTACAGCTTAATGCGGGAGGCGATATTCGCCCTCCGCGTGCAGACGAAGCGGTCAGCGTGC‑GTTTCAAAGCGCCAGGATTGCTGCATATTGTCGGGCGTAACGGCGGACGTAAGCTAAAGAAAATCTGGCAAGAGCTGGGCGTGCCG  >  NC_000913/213192‑213491

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