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

TGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAGGATAAA  >  NC_000913/1032961‑1033248
                                                                                                                                                |                                                                                                                                                  
tGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCAccgccg                                                                                                                                           >  1:799530/1‑151 (MQ=255)
 gTTCCTCCACACGATGGAATGATGGCGTTTCTTTCGCAATCCAGACTGGTCACGGCTGCATGAGCTGTGCGGAAAATGGTTTCTCCGATCGCGGTTCGTTCTGCAGCCGCATGATCGAGATTCCGCAAATGGGTACGCATTccgccg                                                                                                                                              <  1:281097/147‑1 (MQ=255)
      tccACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCg                                                                                                                                     >  2:539798/1‑151 (MQ=255)
        cacaCGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTc                                                                                                                                   >  1:241682/1‑151 (MQ=255)
         acacGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCg                                                                                                                                  >  2:549986/1‑151 (MQ=255)
           acGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGt                                                                                                                                >  2:1109162/1‑151 (MQ=255)
               tGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTa                                                                                                                             >  1:981656/1‑150 (MQ=255)
               tGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTa                                                                                                                             >  2:550222/1‑150 (MQ=255)
               tGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTa                                                                                                                             >  2:143867/1‑150 (MQ=255)
                  aattgtgGTGTTTCTTCCCCAACCCAGTGGTGCCACGGCCGCCTGGGCGGGGGGGAAATTGTTTTCGGGGATCGCGTTCCTTTCTCCAGCCGTGTGGTCTATTTCCCGAAAAGGGCTCCTCTTTCCACCGCCGATACCGTCGGTTTAAcc                                                                                                                          <  1:763676/145‑1 (MQ=255)
                  aaTGTTGGGTTTTTTTTCCCAACCCAGTTGGTTCAGGGCCGCCGGGGTTGGGCGGAAAATGGTTTCGGGGATCGGGTTCCGTTCAACACCCGTGTGTTCTATATCCCCAAATTGGGTCCTTTTCCCCCCGCCGATACCGTCGGTTTAACcg                                                                                                                         <  2:711878/151‑1 (MQ=255)
                        ggCTTTTTTTCCCCAACCCAGTCTGGTCACGGCTGCCGGGGTTGTGGGAAAAATGTTTTCGGGGATCGCGGTTCTTTCTACAGCCGGGTGGTCGATTTCCCGAAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACcgcg                                                                                                                       <  2:1186442/147‑1 (MQ=255)
                          tttttcttccCCAACCCAGTCGGGTCAGGGCTGCGGGGGCTGTGCGAAAAATGGTTTGGGGGATCGCGGTTCTTTCTACCGCCGCGTGGTCGATATTCCGAAAATGGGTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTgg                                                                                                                  <  2:52981/148‑1 (MQ=255)
                            ttttttccCCAACCCAGGCGGGCCCCGGCTCCCGGGGGTGTGGGAAAATTGTTTTCGGGGCTCGCGGTTCGTTCTACGGCGGCGGGGTGAATATCCCAAAATGGGGTCCTCACCCCACCGCCGATACCGTCGGTTTAACCGCGCTTgg                                                                                                                  <  1:768079/148‑1 (MQ=255)
                             ttttttcCCAACCCAGTCTGGTCACGGCTGCCTGGGCTGTGGGGAAAATGTTTTCGGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCAATATCCCGCAAGTGGTTACTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCgt                                                                                                               <  2:857563/147‑1 (MQ=255)
                              tctttcCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCgtgg                                                                                                             >  1:228082/1‑151 (MQ=255)
                                 ttcCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGtggtgg                                                                                                          >  1:1140170/1‑151 (MQ=255)
                                          cAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTg                                                                                                 >  1:550835/1‑151 (MQ=255)
                                           aGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGt                                                                                                >  2:134710/1‑151 (MQ=255)
                                                gttcccGGCTGCCGGGGCTGGGGGAAAAATGTTTTTGGGGATCGCGTTTCTTTCTACGGCGGGGTGGTCGTTATCCCGAAAGGGGTTCCTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTg                                                                                              <  1:1173043/143‑1 (MQ=255)
                                                ggTCGCGGCTGCGGGGGCTGGGCGAAAATTGTTTTCGGGGATCGCGGTTCGTCCTACCGCCGCGTGGTCGATATCCCGAAAAGGGGTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTgg                                                                                             <  1:169742/149‑1 (MQ=255)
                                                ggTCACGGCTGCCTGGGGTGTGCGGAAAATGTTTTCTGGGATCGCGGTTCTTTCTACAGCCGCGTGGTCGATATCCCGAAATTGGGTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGtt                                                                                               <  2:268434/147‑1 (MQ=255)
                                                    aCGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGc                                                                                        >  1:510644/1‑150 (MQ=255)
                                                             tGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCg                                                                               >  1:776270/1‑150 (MQ=255)
                                                             tGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGc                                                                              >  1:673159/1‑151 (MQ=255)
                                                              gggCTGTGCGAAAAGTGGTTTCTGGGCCCGCGGTTCTTCCCACGGCGGCGGGGTCGATATCCCGAAAAGGGGTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTc                                                                                <  2:402730/148‑1 (MQ=255)
                                                              gggCTGTGCGAAAAATGTTTTCGGGGTCCGGGGTTCGTTCAACACCCGGGGGGTCGATATCCCGAAAAGGGGTACTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCg                                                                               <  2:379853/149‑1 (MQ=255)
                                                              gggCTGGGGGAAAAATGTTTTCGGGGAGCGGGTTTCTTTCTACGGCCGGGTGGTCAATATCCCGAAAAGGGGTCCTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAg                                                                                  <  2:903238/146‑1 (MQ=255)
                                                               gggttgtgggAAAATGGTTTTCGGGGACCGCGTTTCTTTCAACGGCGGGGGGGTCGATTTCCCGAAAATGGGTCCTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGcc                                                                             <  2:646485/143‑1 (MQ=255)
                                                               ggCGGTGGGAAAATTGTTTTCGGGGACCGCGTTTCTTCCTACAGCCGCGTGTTCGATATCCCGAAAGTGGGTACTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCa                                                                            <  2:286188/151‑1 (MQ=255)
                                                                 cTGTGCAAAAAATTTTTTCTGGGGCGGCGGTTCGTTCCGCCGCGGCGTGGTCTTTTTTCCGAAATTGGGTCCTCATCCCGCCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCa                                                                            <  1:271889/149‑1 (MQ=255)
                                                                   gtgCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGc                                                                              >  1:1022111/1‑145 (MQ=255)
                                                                    tgCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGcc                                                                       >  1:579274/1‑151 (MQ=255)
                                                                    tgCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGcc                                                                       >  1:1176216/1‑151 (MQ=255)
                                                                        gAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTg                                                                   >  2:1043691/1‑151 (MQ=255)
                                                                          aaaTGGTTTTGGGGGTCCGCGGTTCTTTCTACGGCGGGGTGGTCAATTTCCCGAAAGGGGGTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGtt                                                                    <  2:306015/149‑1 (MQ=255)
                                                                            aTGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGAcc                                                                >  1:521831/1‑150 (MQ=255)
                                                                                     tGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGAc                                                       >  2:846268/1‑151 (MQ=255)
                                                                                         aTCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTc                                                    >  2:360421/1‑150 (MQ=255)
                                                                                             gcgGTTCTTTCAACGGCCGGGGGTTCAATATCCCGAAAGGGGGTCCCCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATaa                                                <  2:456714/151‑1 (MQ=255)
                                                                                               ggTTCGTTCTCCAGCCGCGTGGCCGATTTCCCGAAATGGGGTACTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAAc                                               <  2:64491/150‑1 (MQ=255)
                                                                                                 ttctttctACAGCCGGGTGGTCGATTTCCCGAAAATGGTTACTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAAc                                               <  2:58524/148‑1 (MQ=255)
                                                                                                    gttcTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCa                                          >  1:475374/1‑150 (MQ=255)
                                                                                                     ttctccggcggggTGGTCGATATCCCGAAAGTGGGTACTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAAcc                                              <  2:649094/133‑1 (MQ=255)
                                                                                                       ctacgcccGCGTGGTCAATTTCCCGAAAATGGCTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAg                                            <  1:797598/139‑1 (MQ=255)
                                                                                                       cTACAGCGGGGGGGTCAATATCCCGAAAAGGGTTCCCCACCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGc                                           <  1:1168214/146‑1 (MQ=255)
                                                                                                          cggcggggTGGTCGATATCCGAAAATGGGGTACCAATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAAcc                                       <  2:731613/140‑1 (MQ=255)
                                                                                                           aGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTAc                                    >  2:139765/1‑149 (MQ=255)
                                                                                                           aGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTAc                                    >  1:111725/1‑149 (MQ=255)
                                                                                                             ccGCGTGGTCAATATCCCGCAAATGGGTACTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGaa                                <  2:482587/151‑1 (MQ=255)
                                                                                                                ggtgctggttatTCCAAAAATGGCTTCCCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCt                                      <  2:196848/133‑1 (MQ=255)
                                                                                                                 gTGGTCAATATCCCGAAAATGGGTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTa                                     <  2:1032619/142‑1 (MQ=255)
                                                                                                                   ggTCGATTTTCCGAAAATGGGTCCTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAAc                              <  2:1033553/147‑1 (MQ=255)
                                                                                                                      cGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACAt                       >  2:257684/1‑151 (MQ=255)
                                                                                                                         ttccccGAAAGTGGTTACTCACCCCGCCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACa                        <  1:865962/143‑1 (MQ=255)
                                                                                                                           ttCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGcc                  >  1:1081205/1‑151 (MQ=255)
                                                                                                                             cccaaaagttgctccTCATCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCa                     <  1:36887/133‑1 (MQ=255)
                                                                                                                              cGCAAATGGTTACTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAg                <  1:291239/150‑1 (MQ=255)
                                                                                                                                aaaaTGGGTCCTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCa                     <  2:885728/142‑1 (MQ=255)
                                                                                                                                  aaTGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCaa            >  1:1153501/1‑150 (MQ=255)
                                                                                                                                   aTGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATg          >  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)
                                                                                                                                                |                                                                                                                                                  
TGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAGGCAATGAGGATAAA  >  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: