Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,033,0870CA42.6% 97.9 / 27.2 61Q317K (CAA→AAA) hyaAhydrogenase 1, small subunit
Reads supporting (aligned to +/- strand):  ref base C (29/6);  new base A (0/26);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.37e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGGGCCTACCACCTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAG  >  NC_000913/1032937‑1033234
                                                                                                                                                        |                                                                                                                                                    
aaaGGGCCTACCACCTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGc                                                                                                                                                      >  1:949538/1‑151 (MQ=255)
 aaGGGCCTACCACCTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGc                                                                                                                                                      >  2:584974/1‑150 (MQ=255)
 aaGGGCCTACCACCTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCa                                                                                                                                                     <  1:82360/151‑1 (MQ=255)
     gCCTACCACCTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑g                                                                                                                                                >  2:285011/1‑151 (MQ=255)
        tACCACCTATAACGCCTGTTCCTCAACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGt                                                                                                                                             >  2:797601/1‑151 (MQ=255)
            accTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTc                                                                                                                                         <  2:817330/151‑1 (MQ=255)
                        tGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCAccgccg                                                                                                                             >  1:799530/1‑151 (MQ=255)
                         gTTCCTCCACACGATGGAATGATGGCGTTTCTTTCGCAATCCAGACTGGTCACGGCTGCATGAGCTGTGCGGAAAATGG‑TTTCTCCGAT‑CGCGGTTCGTTCTGCAGCCGCATGATCGAGATTCCGCAAAT‑GGGTACGCATTccgccg                                                                                                                                <  1:281097/147‑1 (MQ=255)
                              tccACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCg                                                                                                                       >  2:539798/1‑151 (MQ=255)
                                cacaCGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTc                                                                                                                     >  1:241682/1‑151 (MQ=255)
                                 acacGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCg                                                                                                                    >  2:549986/1‑151 (MQ=255)
                                   acGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGt                                                                                                                  >  2:1109162/1‑151 (MQ=255)
                                       tGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTa                                                                                                               >  1:981656/1‑150 (MQ=255)
                                       tGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTa                                                                                                               >  2:550222/1‑150 (MQ=255)
                                       tGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTa                                                                                                               >  2:143867/1‑150 (MQ=255)
                                          aattgtgGTGTTTCTTCCCCAACCCAGTGGTGCCACGGCCGCCTGGGCGGGGGGGAAATTGT‑TTTCGGGGAT‑CGCGTTCCTTTCTCCAGCCGTGTGGTCTATTTCCCGAAAAG‑GGCTCCTCTTTCCACCGCCGATACCGTCGGTTTAAcc                                                                                                            <  1:763676/145‑1 (MQ=255)
                                          aaTGTTGGGTTTTTTTTCCCAACCCAGTTGGTTCAGGGCCGCCGGGGTTGGGCGGAAAATGG‑TTTCGGGGAT‑CGGGTTCCGTTCAACACCCGTGTGTTCTATATCCCCAAATT‑GGGTCCTTTTCCCCCCGCCGATACCGTCGGTTTAACcg                                                                                                           <  2:711878/151‑1 (MQ=255)
                                                ggCTTTTTTTCCCCAACCCAGTCTGGTCACGGCTGCCGGGGTTGTGGGAAAAATGT‑TTTCGGGGAT‑CGCGGTTCTTTCTACAGCCGGGTGGTCGATTTCCCGAAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACcgcg                                                                                                         <  2:1186442/147‑1 (MQ=255)
                                                  tttttcttccCCAACCCAGTCGGGTCAGGGCTGCGGGGGCTGTGCGAAAAATGG‑TTTGGGGGAT‑CGCGGTTCTTTCTACCGCCGCGTGGTCGATATTCCGAAAAT‑GGGTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTgg                                                                                                    <  2:52981/148‑1 (MQ=255)
                                                    ttttttccCCAACCCAGGCGGGCCCCGGCTCCCGGGGGTGTGGGAAAATTGT‑TTTCGGGGCT‑CGCGGTTCGTTCTACGGCGGCGGGGTGAATATCCCAAAATG‑GGGTCCTCACCCCACCGCCGATACCGTCGGTTTAACCGCGCTTgg                                                                                                    <  1:768079/148‑1 (MQ=255)
                                                     ttttttcCCAACCCAGTCTGGTCACGGCTGCCTGGGCTGTGGGGAAAATGT‑TTTCGGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCAATATCCCGCAAGT‑GGTTACTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCgt                                                                                                 <  2:857563/147‑1 (MQ=255)
                                                      tctttcCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCgtgg                                                                                               >  1:228082/1‑151 (MQ=255)
                                                         ttcCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGtggtgg                                                                                            >  1:1140170/1‑151 (MQ=255)
                                                                  cAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTg                                                                                   >  1:550835/1‑151 (MQ=255)
                                                                   aGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGt                                                                                  >  2:134710/1‑151 (MQ=255)
                                                                        gttcccGGCTGCCGGGGCTGGGGGAAAAATGT‑TTTTGGGGAT‑CGCGTTTCTTTCTACGGCGGGGTGGTCGTTATCCCGAAAGG‑GGTTCCTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTg                                                                                <  1:1173043/143‑1 (MQ=255)
                                                                        ggTCGCGGCTGCGGGGGCTGGGCGAAAATTGT‑TTTCGGGGAT‑CGCGGTTCGTCCTACCGCCGCGTGGTCGATATCCCGAAAAG‑GGGTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTgg                                                                               <  1:169742/149‑1 (MQ=255)
                                                                        ggTCACGGCTGCCTGGGGTGTGCGGAAAATGT‑TTTCTGGGAT‑CGCGGTTCTTTCTACAGCCGCGTGGTCGATATCCCGAAATT‑GGGTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGtt                                                                                 <  2:268434/147‑1 (MQ=255)
                                                                            aCGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGc                                                                          >  1:510644/1‑150 (MQ=255)
                                                                                     tGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCg                                                                 >  1:776270/1‑150 (MQ=255)
                                                                                     tGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGc                                                                >  1:673159/1‑151 (MQ=255)
                                                                                      gggCTGTGCGAAAAGTGG‑TTTCTGGGCC‑CGCGGTTCTTCCCACGGCGGCGGGGTCGATATCCCGAAAAG‑GGGTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTc                                                                  <  2:402730/148‑1 (MQ=255)
                                                                                      gggCTGTGCGAAAAATGT‑TTTCGGGGTC‑CGGGGTTCGTTCAACACCCGGGGGGTCGATATCCCGAAAAG‑GGGTACTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCg                                                                 <  2:379853/149‑1 (MQ=255)
                                                                                      gggCTGGGGGAAAAATGT‑TTTCGGGGAG‑CGGGTTTCTTTCTACGGCCGGGTGGTCAATATCCCGAAAAG‑GGGTCCTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAg                                                                    <  2:903238/146‑1 (MQ=255)
                                                                                       gggttgtgggAAAATGGTTTTCGGGGAC‑CGCGTTTCTTTCAACGGCGGGGGGGTCGATTTCCCGAAAAT‑GGGTCCTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGcc                                                               <  2:646485/143‑1 (MQ=255)
                                                                                       ggCGGTGGGAAAATTGT‑TTTCGGGGAC‑CGCGTTTCTTCCTACAGCCGCGTGTTCGATATCCCGAAAGT‑GGGTACTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCa                                                              <  2:286188/151‑1 (MQ=255)
                                                                                         cTGTGCAAAAAATTT‑TTTCTGGGGC‑GGCGGTTCGTTCCGCCGCGGCGTGGTCTTTTTTCCGAAATT‑GGGTCCTCATCCCGCCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCa                                                              <  1:271889/149‑1 (MQ=255)
                                                                                           gtgCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGc                                                                >  1:1022111/1‑145 (MQ=255)
                                                                                            tgCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGcc                                                         >  1:1176216/1‑151 (MQ=255)
                                                                                            tgCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGcc                                                         >  1:579274/1‑151 (MQ=255)
                                                                                                gAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTg                                                     >  2:1043691/1‑151 (MQ=255)
                                                                                                  aaaTGG‑TTTTGGGGGTCCGCGGTTCTTTCTACGGCGGGGTGGTCAATTTCCCGAAAGG‑GGGTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGtt                                                      <  2:306015/149‑1 (MQ=255)
                                                                                                    aTGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGAcc                                                  >  1:521831/1‑150 (MQ=255)
                                                                                                             tGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGAc                                         >  2:846268/1‑151 (MQ=255)
                                                                                                                 aT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTc                                      >  2:360421/1‑150 (MQ=255)
                                                                                                                     gcgGTTCTTTCAACGGCCGGGGGTTCAATATCCCGAAAGG‑GGGTCCCCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATaa                                  <  2:456714/151‑1 (MQ=255)
                                                                                                                       ggTTCGTTCTCCAGCCGCGTGGCCGATTTCCCGAAATG‑GGGTACTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAAc                                 <  2:64491/150‑1 (MQ=255)
                                                                                                                         ttctttctACAGCCGGGTGGTCGATTTCCCGAAAAT‑GGTTACTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAAc                                 <  2:58524/148‑1 (MQ=255)
                                                                                                                            gttcTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCa                            >  1:475374/1‑150 (MQ=255)
                                                                                                                             ttctccggcggggTGGTCGATATCCCGAAAGT‑GGGTACTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAAcc                                <  2:649094/133‑1 (MQ=255)
                                                                                                                               ctacgcccGCGTGGTCAATTTCCCGAAAAT‑GGCTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAg                              <  1:797598/139‑1 (MQ=255)
                                                                                                                               cTACAGCGGGGGGGTCAATATCCCGAAAAG‑GGTTCCCCACCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGc                             <  1:1168214/146‑1 (MQ=255)
                                                                                                                                  cggcggggTGGTCGATA‑TCCGAAAATGGGGTACCAATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAAcc                         <  2:731613/140‑1 (MQ=255)
                                                                                                                                   aGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTAc                      >  1:111725/1‑149 (MQ=255)
                                                                                                                                   aGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTAc                      >  2:139765/1‑149 (MQ=255)
                                                                                                                                     ccGCGTGGTCAATATCCCGCAAAT‑GGGTACTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGaa                  <  2:482587/151‑1 (MQ=255)
                                                                                                                                        ggtgctggttatTCCAAAAAT‑GGCTTCCCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCt                        <  2:196848/133‑1 (MQ=255)
                                                                                                                                         gTGGTCAATATCCCGAAAAT‑GGGTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTa                       <  2:1032619/142‑1 (MQ=255)
                                                                                                                                           ggTCGATTTTCCGAAAAT‑GGGTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAAc                <  2:1033553/147‑1 (MQ=255)
                                                                                                                                              cGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACAt         >  2:257684/1‑151 (MQ=255)
                                                                                                                                                 ttccccGAAAGT‑GGTTACTCACCCCGCCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACa          <  1:865962/143‑1 (MQ=255)
                                                                                                                                                   ttCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGcc    >  1:1081205/1‑151 (MQ=255)
                                                                                                                                                      cGCAAAT‑GGTTACTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAg  <  1:291239/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
AAAGGGCCTACCACCTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAG  >  NC_000913/1032937‑1033234

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