Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913433,6340TG27.6% 104.6 / 21.2 89G60G (GGT→GGGribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (13/50);  new base G (24/0);  total (37/52)
Fisher's exact test for biased strand distribution p-value = 2.08e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAA  >  NC_000913/433398‑433853
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                              
aGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGGGAAGGTTACCACCAACGGGCGGGGGAACAACATGCCGAAGTAAACGCGTTGCGTATGGCGGGGGa                                                                                                                                                                                                                              >  1:427972/1‑239 (MQ=255)
  tttCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTAATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGACGGACACGCGTTACGTATGGCCGGTGgaa                                                                                                                                                                                                                            >  2:841211/1‑239 (MQ=255)
  tttCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGGGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGAACACGCGTTGCGTATGGCGGGGGaaa                                                                                                                                                                                                                            >  1:1119474/1‑239 (MQ=255)
  tttCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGGGCGTCATTGTCAAAGAGGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGGGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaaa                                                                                                                                                                                                                            >  1:1813020/1‑239 (MQ=255)
   ttCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaa                                                                                                                                                                                                                             <  1:1198387/237‑1 (MQ=255)
   ttCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaa                                                                                                                                                                                                                             <  1:841211/237‑1 (MQ=255)
        aGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACGCGCGTTGCGTATGGCGGGTGAAAAAGCCa                                                                                                                                                                                                                      <  1:1780583/239‑1 (MQ=255)
        aGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGcc                                                                                                                                                                                                                       <  1:125372/238‑1 (MQ=255)
            atCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGGAAAAGCCaaaag                                                                                                                                                                                                                  >  1:374937/1‑237 (MQ=255)
            atCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGGGAACCACATGCCGAAGAACACGCGTTGCGTATGGCGGGGGGAAAAACCaaaag                                                                                                                                                                                                                  >  1:893197/1‑237 (MQ=255)
                  gaagtgggagAAGCGAGCGCGCCCCTGGAGGACTAAGCCGTGAAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATGGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGtc                                                                                                                                                                                                                <  2:893197/226‑2 (MQ=255)
                     tttGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGGGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCAAAAGGGCGAccc                                                                                                                                                                                                          >  1:1646990/1‑237 (MQ=255)
                                 aTCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGATTATTACATGGCGCGGGCGATAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAGGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCacac                                                                                                                                                                                            <  1:1045321/238‑1 (MQ=255)
                                   gggccgcCTGGAGGACTAAGCCGTGCAGGAGGATTATTACAGGGCGCGGGCGCTAGAGCTGGCGCAACGAGGGCGTTTTACCACGCACCCCAACCCGAATGGCGGGTCCGTCATTGACAAAGAGGGCGAAATTGTCGGGGAAGGTTACCACCAACGTGCGGGTGAACCACAGGCCGAAGTACCCGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCg                                                                                                                                                                                         <  2:231026/235‑1 (MQ=255)
                                         ccTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGCGACCGCCTATGTCACACTCGAAcccc                                                                                                                                                                                   >  1:1321003/1‑238 (MQ=255)
                                            gtaggaCTAGGCCGTGCAGTACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGGTTTACCACGCACCCCAACCCGAATGTCGGGTGCGTCATTGTAAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGTTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTAAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                <  2:2092323/237‑1 (MQ=255)
                                                gaCTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGGGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGAGGGGTGAAAAACCAAAAGGGCGACCCCCCATGTCACAACCCGAAcctggaggcc                                                                                                                                                                             >  1:2055703/1‑231 (MQ=255)
                                                   cacgcagTGCAGGGCTGGTATTCGAGGGCGAGGGCGTAAAACTGGGCGCACCGGGTCCTTTTTCCCCGCCATCCAACCCCAAATTGCGGGGGGGTAATTGCAAAAGAGGGCGAAATGGTGGGTGAGGGTTCCCCCCACGGGGCGGGTGACCCCCATGCCAAAGTCCCCGGGTTGCGTATGGGGGGTAAAAAAGCCAAGGGTGGGACCGCCTATGTCACACTCGAACCCTGTAGCcatca                                                                                                                                                                         <  2:667163/233‑1 (MQ=255)
                                                    aaGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGAGCGGCCGCCTATGCCCCCCTCCAACCCTCTAGCcctctt                                                                                                                                                                        >  1:223427/1‑237 (MQ=255)
                                                          tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCCATGTCACACTCGAAACCTGTAGCCCACATGGGCGt                                                                                                                                                                  >  1:407712/1‑239 (MQ=255)
                                                          tGCAGGACGAGTATTACATGGCGCGGGCCCTAAAGCGGGCGCACCGAGGACGTTTTACCACGCATCCAACCCCAAATTTCGGGTGCGTCATGGCCAAAGATGGCGAAATGGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACAGGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                      <  2:1321003/235‑1 (MQ=255)
                                                                 tgagttgtACTTGGCCGGGGCGTTAACGCGGGCGCAAGGAGGACGTTTACCCTCGCATCCCAACCCGAGTGGCGGGTGCTTCGTTGCCAAAGATGGCAGAATTTTCGGTGGAGTTGCCCACCAGCGTGGGGTTGACCCACGTGCCGACGTCCCCCTTTTGCGTATGGCGGGTGAAAAGCCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                            <  2:2055703/231‑1 (MQ=255)
                                                                 cGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGTGGGACCGCCTCTGTCACACCCCAACCCTGGAGCCATCATGGGTGTAAGCCa                                                                                                                                                            >  1:1595137/1‑238 (MQ=255)
                                                                    gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGGGACCGCCTCTGTCACACCCCAACCCTGTAGCCAACATgg                                                                                                                                                                      >  1:42220/1‑225 (MQ=255)
                                                                    gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGGGACCCCCCATGTCACAACCCAACCCCGTAGCCAACATGGGTGTACCCCCCCgg                                                                                                                                                        >  1:183183/1‑238 (MQ=255)
                                                                    gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGCAACCATCATGGTCGTACGCCACCGt                                                                                                                                                        >  1:1072880/1‑239 (MQ=255)
                                                                    gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCACAATCGAACCCTGTACCCCACCTGGTCGTACGCCAcccg                                                                                                                                                        >  1:42878/1‑237 (MQ=255)
                                                                     tcttactgggCGCGGGCTTTAAGGCTGGCGCACCGAGGACGTTTTACTGCGCTCCACACCCGAATTGCCGGGGCCTCCTTTGTCAAAGGGGGCAAATGTGTCGTGAAAGTTTCCCCCCACCGTGCGGGCGCCCCCCTGGCCGACGCCCCCGCGTGGTGTTGGGCGGGAAAAAAAGCCAAGGGTGCGCCCGCCTATGTCACCCTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                       <  2:1813020/231‑1 (MQ=255)
                                                                      aTTACATGGCGGGGGCGATAAAGCTGGCGCAACGAGGACGTTTTCCCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                       <  1:1796735/238‑1 (MQ=255)
                                                                         aCATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCCAAAGGGCGACCCGCTATGGCACCACCGAAcc                                                                                                                                                                                     >  1:1988084/1‑205 (MQ=255)
                                                                               cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAAGCCAACGTGGGGCCGCCCCATGCCAACTCGAAACCCGTAGCCATCCCCGGCCGACGCCACCCTGCtgttcagc                                                                                                                                              >  2:699433/1‑233 (MQ=255)
                                                                               cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGGGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGGAAAAGCCAAAGGGGCGACCGCCCCTGTCACACTCGCACCCCGTAACCATCCAGGTAGGTAGCCACCGAGCTTGGACGc                                                                                                                                              >  2:1256657/1‑238 (MQ=255)
                                                                                cgggcgtcaaaaGGTGGCCCAACGATGCCGTTGTCCCACGCATCCCAACCCGATGGTCGGGTGCGTCATGGCCAAAGATGGCAAAATGGACGGGGAGGGTCACCACCAACGTGCGGGTGACCCACATGCCGAATTACCCTCGTTCCGTATGGCGGGTGAAAAAGCCAAAGGTGCTACCGCCTATGTCACACTCGACCCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCAc                                                                                                                                            <  2:407712/230‑1 (MQ=255)
                                                                                    gcgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAGGGGCGACCGCCTATGTCACACCCCAACCCTGTAGCCATCATGGGCGTACGCCACCGTGCTGTGACGCACTcc                                                                                                                                         >  1:1946252/1‑237 (MQ=255)
                                                                                              tGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGGGCGAACGCCTAGTTCACCCTCCAACCCCGTAg                                                                                                                                                                               >  1:1523039/1‑190 (MQ=255)
                                                                                                                          cATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                    <  1:1274966/239‑1 (MQ=255)
                                                                                                                           gtcacAACCCGAATGTCGGGTGCGTCGTTGGAAAAGAGTGAGAAAGGTGCGGGGAAGTTGCCCACCAACGTGCGGGTGAACCACATGCCGAAGAGCACGCGTGGGGTATGGGGGGTAAAAACACAAGGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                    <  2:115077/234‑1 (MQ=255)
                                                                                                                                       aTGGCGGGTGCGTCATTGTCAAAGAGGGCGAAATTGGCGGTGAAGGTTCCCACCAACGTGCGGGTAAACCACATGCCGAAGTACACGCGTTGCGTAGGGGGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGaa                                                                                                                                                                                       <  1:2024093/143‑1 (MQ=255)
                                                                                                                                                    tCATTGGAAAAGAGGCCGAATTTGTGGGTGAAGGTTACCACCAACGGGGGGGTGACCCACATGCCGAAGTACCCGGTTTGCTTATGGGGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTc                                                                                                                                                                                          <  2:807205/128‑1 (MQ=255)
                                                                                                                                                     taTTGTCAAAGATGGCGAAATTGTGGGGGAAGGTTCCCACCAACGTGCGGGTGAACCACAGGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGc                                                                          <  1:1256657/238‑1 (MQ=255)
                                                                                                                                                        tGTCAAAGATGGCGAAATTGTCGGTGAAGTTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGc                                                                                                                                                                              <  1:1194100/136‑1 (MQ=255)
                                                                                                                                                        gcataaatggggCGAAGTTGGCGGGGATGTTGACCCCCACCGGGCGGTGTAACCCCTGGCCGAAGTCCCCGCGTTGCGTGTGGGGGGAGAAAACGCAAAAGGTGCCACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                        <  2:1119474/228‑1 (MQ=255)
                                                                                                                                                          tCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaa                                                                                                                                                                                                                             <  1:1836266/89‑1 (MQ=255)
                                                                                                                                                          tCAAAGAGGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGGGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                       <  1:2030007/157‑1 (MQ=255)
                                                                                                                                                          acAAGGATGGAGAAGTGGCGGGTGAAGTTTCCCCCCACGGGGCGGGGAAACCGCAGGCCGAAGTCCCCGCTTTGCGAAGGGCGGGTGAAAACGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTc                                                                                 <  2:1390484/226‑1 (MQ=255)
                                                                                                                                                              aGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGGGCGACCGCCTCTGTCACACTCGAACCCTGTa                                                                                                                                                                                >  1:807205/1‑128 (MQ=255)
                                                                                                                                                                aTGGCGAATTTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGtt                                                                                                               <  1:661711/191‑1 (MQ=255)
                                                                                                                                                                  ggCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaaaa                                                                                                                                                                                                                           <  2:72705/83‑1 (MQ=255)
                                                                                                                                                                  ggCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaaaa                                                                                                                                                                                                                           >  1:72705/1‑83 (MQ=255)
                                                                                                                                                                    cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCCAAGGGGCGACCCCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCCGGCGTAGCGCGCGTGGGTGCCTCGATGCAAGATCCTAACCCGCAGGCCGCTGGGCggg                                                                       >  1:1390484/1‑225 (MQ=255)
                                                                                                                                                                          tGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGAATGGCGGGTGGAAAAGCCAAAGGCGGCACCGCCGATTTCTCACTCGGCCCCTGGAACCCTCCTGGTCGTAGCTCAACGTGCTGTGACCCCATCACCCCCACCTGTGTACCGCGGCTGGCTGCCTCCCAGCAAGAACCTCACCCCCCGGGCGCCGGGCGTGGGACTTATCGTCTGGAAc                                                    >  2:396180/1‑239 (MQ=255)
                                                                                                                                                                             cGGTGAGGGTTACCACCAACGTGCGGGTGAACCACAGGCCGAAGTACCCGCGTTGCGTAGGGCGGGTAAAAAAGCCAAAGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAgg                                                 <  1:699433/239‑1 (MQ=255)
                                                                                                                                                                             cGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAGGGGCGACCGCCTCTGTCACCCTCGAACCCCCTAGCCACCATGGTCGTACGCCACCCCGCTGTGACGCACCCCTCGCCGCTGGCGTAGCGGGCGGGGGTGCCTCGCTGCAAGAGCCCAACCCGCAGGTGCCTGGGCGCGGGCTTTACCGTCTGCAACac                                                  >  2:287959/1‑237 (MQ=255)
                                                                                                                                                                              ggTGAAGGTTCCCACCAACGTGCGGGTGAACCACAGGCCGAAGTACACGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                <  2:1211033/112‑1 (MQ=255)
                                                                                                                                                                               gggAAGGTTCCCCCCAACGGGGGGGTGAACCACATGCCGAAGTCCACGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                    <  1:318831/185‑1 (MQ=255)
                                                                                                                                                                               gggAAGGTTACCCCCAAGGGGCGGGTGACCCACATGCCGAAGTACACGCTTTGCGTAGGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGc                                                <  1:287959/236‑1 (MQ=255)
                                                                                                                                                                                     gTTCCCACCAACGTGCGGGTAAACCACAGGCCGAAGTCCACGCGTTGCGTAGGGCGGGTAAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCa                                                                                                                                             <  2:592353/140‑1 (MQ=255)
                                                                                                                                                                                     gTTACCACAAAGGGGGGGGGGAAACACATGCCGAAGTCCAGGAGTGGCGTGTGGGGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc                                                                               <  1:983603/202‑1 (MQ=255)
                                                                                                                                                                                      ttACCACCAACGTGCGGGTAACCCCCATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCAt                                          <  1:1051498/237‑1 (MQ=255)
                                                                                                                                                                                        accaccAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGGGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                      >  1:1211033/1‑112 (MQ=255)
                                                                                                                                                                                         ccaccaACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCCAAGGGGGGACCGCCCATGTCCCACTCCAACCCCGTAGGCCTCATGGTCGTACGCCACCGTGCTGTGACGCCCTCCTCCCCCCCGGCGTCGCGCGCGCGGGTGCCTCGATGCAAGATCCTAAc                                                                                          >  2:318831/1‑187 (MQ=255)
                                                                                                                                                                                             caACGGGCGGGTGAACCCCATGCCGAATTCCCCGCGTTGCGTAGGGGGGTTAAAAAACCAAAAGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGc                                                                          <  1:963528/199‑1 (MQ=255)
                                                                                                                                                                                               aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATcgccg                                                                                                                                   >  1:592353/1‑140 (MQ=255)
                                                                                                                                                                                                 tggggggTGAACCACTAGCCGAAGTACTCGGTTGGCGTAGGGGGGGTGAAAAAGCCAAGGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGcc                               <  1:741296/233‑1 (MQ=255)
                                                                                                                                                                                                 gTGCGGGTGACCCACATGCCGAAGTACCCGCGTTGCGTATGGGGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGa                                                                                                                                                                                        <  2:391437/85‑1 (MQ=255)
                                                                                                                                                                                                  tGCGGGTGAACCACATGCCGAAGTACACGCGTGGCGTATGGGGGGTGAAAAAGCAAAAGGTGCGTCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGc                                                                                                                                                      <  1:218537/118‑1 (MQ=255)
                                                                                                                                                                                                   gggggTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg                            <  1:708796/236‑1 (MQ=255)
                                                                                                                                                                                                       gTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTg                                              <  1:1034858/216‑1 (MQ=255)
                                                                                                                                                                                                       gTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCCCCCCCGAACCCCGTGGCCATCATGGGCGGACGCCACCGTGCTGTGACGCACCCATCGCCCCCGGCGTCGCGCGCGTGGTTGCCTTCGTTCAAGATCCTTACCCGCAGGGCGCTGGGCGGGGGCGTTa                                                               >  2:963528/1‑199 (MQ=255)
                                                                                                                                                                                                        tGAACCACATGCCGAAGTCCCCGCGTTGCTTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCtg                       <  1:1453170/238‑1 (MQ=255)
                                                                                                                                                                                                          aaCCACATGCCGAAGTACGCGCGGTGGGTGTGGGGGGTGAAAAAGCAAAAGGGGCGACCGCCTATGTCACTCTCGAACCCGGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGc                                                <  1:1891726/211‑1 (MQ=255)
                                                                                                                                                                                                           aCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGc                                                                                                                                                                              >  1:391437/1‑85 (MQ=255)
                                                                                                                                                                                                            ccACATGCCGAAGTACACGCGTTGCGTTTGGCGGGTAAAAAAGACAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCcatcattg                                                                                                                                                                      <  2:1723240/92‑3 (MQ=255)
                                                                                                                                                                                                             cacaGGCCGAAGTACCCGCGTTGCGTAGGGCGGGTAAAAAACCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACg                                                                                                                                               <  2:1404617/114‑1 (MQ=255)
                                                                                                                                                                                                                 tGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGGCAAAGGGGCGACCGCCTATGTc                                                                                                                                                                                                 >  1:919679/1‑61 (MQ=255)
                                                                                                                                                                                                                 tGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGGGGCCCCCTATGTCACACTCGAACCCTGTAGCCATCAAGGGCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGGCTCGGTGCCAGATCCTAAAACCCCGGGCGCCGGGGGGGGGCTTTTACCGCTGCCACAGGCTCGGCTTGACGt                                    >  2:1034858/1‑216 (MQ=255)
                                                                                                                                                                                                                    cgattaacacGCTTTGCGTATGGCGGGTGAAAAAGCCAAAGGGCCAACCGCCTATGTCACCCTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAg          <  1:1340684/233‑1 (MQ=255)
                                                                                                                                                                                                                    caagtcacgcgcgTTGTGTAGGGCGGGTGAAAAACCCAAGGGTCGTCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAg          <  1:668234/230‑1 (MQ=255)
                                                                                                                                                                                                                      aaGTCCCCGCGTTGCGTAGGGCGGGTAAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGc                          <  1:396953/221‑1 (MQ=255)
                                                                                                                                                                                                                      aaGTACACGCGTTGCGTATGGCGGGTGAAAAAACCCAAGGGGCGACCGCCCATGTCACACCCGAACCCTGTAGCCATCATGGGCGGACGCCa                                                                                                                                                            >  1:1723240/1‑92 (MQ=255)
                                                                                                                                                                                                                       aGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACCCGAAACCTGTAGCCATCATGGTCGTACGCCACCCGGCTGTGACGCACTCATcgc                                                                                                                                     >  1:1404617/1‑114 (MQ=255)
                                                                                                                                                                                                                        gTACAGGCGTTGCGTAGGGCGGGTAAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCcatca                                                                                                                                                                         <  1:470955/77‑1 (MQ=255)
                                                                                                                                                                                                                         tccacGCTTTGGGTTTGGCGGGTGAAATAGCCAAAGGTGCGACCGCCTATGTCTCCCTCGAACCCTGTATCcatcat                                                                                                                                                                        <  1:1450691/75‑1 (MQ=255)
                                                                                                                                                                                                                            cgtcgTTGGTTAGGGCGGTTAAAAAATCCAAAGGTCCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa  <  1:396180/236‑1 (MQ=255)
                                                                                                                                                                                                                                gTTGCGTATGGCGGGTGAAAAAACCAAAGGTGCGACCGCCCATGTCACACCCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATCGCCCCCGGCGGAGCGCGCGTGGGTGCCTCGATGCCAGATCCTAACCCGCCGGGCGCTGGGCGGGGGCCTTTCCGTCCGCAACCGGCTGGCACTGGCCGCAGGCACGGCCtggtgatga                >  2:396953/1‑221 (MQ=255)
                                                                                                                                                                                                                                    cGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                    <  1:696489/134‑1 (MQ=255)
                                                                                                                                                                                                                                      tATGGGGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc                                                                               <  1:842392/153‑1 (MQ=255)
                                                                                                                                                                                                                                                gTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                <  2:1061413/48‑1 (MQ=255)
                                                                                                                                                                                                                                                gTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                >  1:1061413/1‑48 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                              
AGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAA  >  NC_000913/433398‑433853

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