Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,033,1020TC46.3% 83.6 / 24.7 67S322P (TCC→CCC) hyaAhydrogenase 1, small subunit
Reads supporting (aligned to +/- strand):  ref base T (33/3);  new base C (0/31);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 5.46e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAGGATAAA  >  NC_000913/1032961‑1033248
                                                                                                                                                |                                                                                                                                                  
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:579274/1‑151 (MQ=255)
                                                                    tgCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGcc                                                                       >  1:1176216/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                                    >  2:139765/1‑149 (MQ=255)
                                                                                                           aGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTAc                                    >  1:111725/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)
                                                                                                                             cccaaaagttgctccTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCa                     <  1:36887/133‑1 (MQ=255)
                                                                                                                              cGCAAAT‑GGTTACTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAg                <  1:291239/150‑1 (MQ=255)
                                                                                                                                aaaaT‑GGGTCCTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCa                     <  2:885728/142‑1 (MQ=255)
                                                                                                                                  aaT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCaa            >  1:1153501/1‑150 (MQ=255)
                                                                                                                                   aT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATg          >  2:384388/1‑151 (MQ=255)
                                                                                                                                      gggTCCTCATCCCGCCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAg        <  2:626547/151‑1 (MQ=255)
                                                                                                                                      gggTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAAt           <  2:847125/148‑1 (MQ=255)
                                                                                                                                      gggTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATg          >  1:311373/1‑149 (MQ=255)
                                                                                                                                         tACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAGGAt     >  2:847800/1‑151 (MQ=255)
                                                                                                                                            tCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAGGATaa   >  2:977834/1‑150 (MQ=255)
                                                                                                                                             cATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAGGATaaa  >  1:739887/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGG‑TTTCTGGGAT‑CGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAAT‑GGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAGGATAAA  >  NC_000913/1032961‑1033248

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