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

AAAGGGCCTACCACCTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAG  >  NC_000913/1032937‑1033234
                                                                                                                                                        |                                                                                                                                                    
aaaGGGCCTACCACCTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGc                                                                                                                                                      >  1:949538/1‑151 (MQ=255)
 aaGGGCCTACCACCTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGc                                                                                                                                                      >  2:584974/1‑150 (MQ=255)
 aaGGGCCTACCACCTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCa                                                                                                                                                     <  1:82360/151‑1 (MQ=255)
     gCCTACCACCTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATg                                                                                                                                                >  2:285011/1‑151 (MQ=255)
        tACCACCTATAACGCCTGTTCCTCAACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGt                                                                                                                                             >  2:797601/1‑151 (MQ=255)
            accTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTc                                                                                                                                         <  2:817330/151‑1 (MQ=255)
                        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:1176216/1‑151 (MQ=255)
                                                                                            tgCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGcc                                                         >  1:579274/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                      >  1:111725/1‑149 (MQ=255)
                                                                                                                                   aGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTAc                      >  2:139765/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)
                                                                                                                                                      cGCAAATGGTTACTCTTCCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAg  <  1:291239/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
AAAGGGCCTACCACCTATAACGCCTGTTCCTCCACACGCTGGAATGATGGCGTTTCTTTCCCAATCCAGTCTGGTCACGGCTGCCTGGGCTGTGCGGAAAATGGTTTCTGGGATCGCGGTTCGTTCTACAGCCGCGTGGTCGATATTCCGCAAATGGGTACTCATTCCACCGCCGATACCGTCGGTTTAACCGCGCTTGGCGTGGTGGCAGCGGCTGTTGGTGTGCACGCAGTCGCCAGCGCCGTTGACCAGCGCAGACGTCATAACCAGCAACCTACAGAAACCGAACATCAGCCAG  >  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: