Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913433,6340TG31.2% 111.4 / 17.3 80G60G (GGT→GGGribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (14/41);  new base G (19/6);  total (33/47)
Fisher's exact test for biased strand distribution p-value = 5.17e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTC  >  NC_000913/433401‑433871
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
ttCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaa                                                                                                                                                                                                                                                <  2:719653/237‑1 (MQ=255)
      gATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCaa                                                                                                                                                                                                                                        >  1:224014/1‑239 (MQ=255)
        tatCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAg                                                                                                                                                                                                                                     <  1:1322219/239‑1 (MQ=255)
                  tttGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGGGCGACCg                                                                                                                                                                                                                            >  1:205126/1‑238 (MQ=255)
                                 gcgggcCGGGAGGACTAAGCGGGGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGAATCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTAAACCACATGCCGAAGTCCCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCg                                                                                                                                                                                                             <  2:1024995/238‑1 (MQ=255)
                                 gcgcgcCTGGAGGACTAAGCCGTGCAGGACGAGTATTCCATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGa                                                                                                                                                                                                            <  1:146423/239‑1 (MQ=255)
                                   gcgcCTGGAGGACAAGGCCGGGCAGGACGATTATTACAGGGCGCGGGCGCAAAAGCGGGCGAAACGAGGACTTTTTACCACGCACCCAAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCAAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTAGGGGGGGTGAAAAAGCAAAGGGTGCGACCGCCTATGTCACACTCGaa                                                                                                                                                                                                           <  2:205126/238‑1 (MQ=255)
                                         ggaggtATAAGCCGGGCAGGACGATTATTACATGGCGCGGGGGCTAAAGGTGGGGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGGGCGGGTGAACCACATGCCGAAGTCCACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                   <  2:224014/239‑1 (MQ=255)
                                         ggaggaCTAAGCCGTGCAGGACGAGTTTTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGa                                                                                                                                                                                                            <  1:601194/231‑1 (MQ=255)
                                           agcccTAACCCTTGCAGGACCATAATTACATGCCCCGGGCCGTAAAGCTGGCGCACAGGGAGCGTTTTACCCCCCATCCCACCCCTAATGTCGGGTGCGTCATTGTGAAAGAGGGCGAAATGGTGGTTGAAGTTTACCACCACCGGGCGGGTGAGCCTCATGCGAAAGTCCCCGCGTTGCGTATGGGGGGTGAAAAAGCCAAAGGTGCGACCGcc                                                                                                                                                                                                                          <  2:940404/211‑1 (MQ=255)
                                             gaCTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGt                                                                                                                                                                                                                     <  1:1249018/218‑1 (MQ=255)
                                               cTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCCAAAGGGGGAGCCCCCCTTTCACACTCGAACCCtcgaggcccc                                                                                                                                                                                              >  2:373337/1‑230 (MQ=255)
                                                   gCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGGGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTATACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAAGGCGACCGCCCCTGTCGAAACCGAAcccccatgc                                                                                                                                                                                                  >  2:601194/1‑225 (MQ=255)
                                                    ccGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAACGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaa                                                                                                                                                                                                                                                <  1:620763/185‑1 (MQ=255)
                                                     cggagaggaggAGTGATTACGGGGCCCGGGCGCTAAAGCGGGCGCAAGGAGGAGTTTTTACCCCGCTACCCAGCCGGAAGGGGGGTGGAGTCAGTGTCAAAGTTGGGGAACTTGACGGTGGAGGTTCCCACCCCCGTGAGGTTAAGCCACTTGCCAAAGCACACGGTTTGCGTAGGGGGGGTGAAAAACCCAAAGGTCCAACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                         <  2:1080514/233‑1 (MQ=255)
                                                     cGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGCGCGTCATTGTCAACGATGGCGAAATTGACGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGTCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCCAAGGGGCGACCCCCCAAGGcacaa                                                                                                                                                                                                                >  1:940404/1‑214 (MQ=255)
                                                       tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCCAAAGGGGGGACCCCCATGTGCCACTCGCAACCTGTAGCcatccagtgcg                                                                                                                                                                                      >  2:181592/1‑231 (MQ=255)
                                                       tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAGAGGCGAGCGCCCATGTGACACACGaa                                                                                                                                                                                                           >  2:1249018/1‑218 (MQ=255)
                                                              cGAGTATTACGTGGCGCGGGCGATAAAGCGGGCGAAACGAGGACGTTTTACCACGCACCCCAACCCGAATGTCGGGTGCGTCATTGTAAAAGATGGAGAAATTGTCGGTGAAGGTTCCCACCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGGGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                               <  2:158923/238‑1 (MQ=255)
                                                              cGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                               <  1:373337/238‑1 (MQ=255)
                                                                tatATTACAGGGCGCGGGCGATAAAGCTGGCAAAACGAGGACTTTTTACCACGCATCCCAACCCGAATGTGGGGTGCGTCATTGTAAAAGATGGCGAAATTGTCGGTGAGGTTTACCACCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTAGGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCAcc                                                                                                                                                                             <  2:960892/236‑1 (MQ=255)
                                                                 gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCCTCATGGTCGTACGCCACCGt                                                                                                                                                                           >  1:21653/1‑239 (MQ=255)
                                                                           ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaa                                                                                                                                                                                                                                                <  1:40859/163‑1 (MQ=255)
                                                                           ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGaa                                                                                                                                                                                                                                                >  2:40859/1‑163 (MQ=255)
                                                                           ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGAc                                                                                                                                                                   >  1:1293321/1‑238 (MQ=255)
                                                                            gggcgGGCGCTAAAGTGGGCCCAACGAGGACGTTTTCCCACGCATCCCAACCCGAATGTGGGGGGCGTAATTGAAAAAGATGGCGAAATTGTCGGTGAAGGTTCCCCCCACGGGGGGGGTGAACCAGAGGCAGAAGTACACGCGTGGGGTGTGGCGGGGGAAAAACCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCtgtg                                                                                                                                                                     <  2:1285911/233‑1 (MQ=255)
                                                                            gcgcGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGGAAAAGCCCAAGGGGCGGACGCCCCTGTCCCACTCGCACCCCGTAGCCCTCCTGGGCGGACGCCACCGCGCCGTGACGc                                                                                                                                                                 >  2:54714/1‑239 (MQ=255)
                                                                             cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaa                                                                                                                                                                                                                                                <  2:1287238/161‑1 (MQ=255)
                                                                             cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaa                                                                                                                                                                                                                                                >  1:1287238/1‑161 (MQ=255)
                                                                             cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCCCACTCGAACCCTGTAGCCATCATGGTCGGACGCCACCGGGCTGTGACGc                                                                                                                                                                 >  1:190331/1‑238 (MQ=255)
                                                                             cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTCTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCAAAAGGGGGAACGCCTATGTCCACCCCCAACCCCGTAGCCATCATGGGCGGACGCCACCGTGCCGTGAAGc                                                                                                                                                                 >  1:1285911/1‑238 (MQ=255)
                                                                                    gcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAcc                                                                                                                                                                                                                             <  1:993293/172‑1 (MQ=255)
                                                                                              gcgcAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCCAAAGGGGGACCGCCCATGTCACACCCGCACCCCGTTGGCATCATGGGGGTACGCCACCGTGCTGTGGCGCACTCCTccccc                                                                                                                                                      >  2:1002946/1‑231 (MQ=255)
                                                                                                    cGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGGGCGACCCCCTATGTCACACTCGAACCCTGTAGCCCTCATGGTCGTACGCCACCGTGCTGTGACGCACCCCTCGCCGCTGGCCTAgcgc                                                                                                                                          >  1:944312/1‑238 (MQ=255)
                                                                                                            ttttACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCcatca                                                                                                                                                                                            <  1:559087/180‑1 (MQ=255)
                                                                                                            ttttACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTAAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGc                                                                                                                                                                                                                           <  1:726631/150‑1 (MQ=255)
                                                                                                             tttACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCTTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCg                                                                                                                                                                            <  1:631889/195‑1 (MQ=255)
                                                                                                               tACCAGGCATCCCAACCCGAATGTCGGGGGCGTCATTGTAAAAGATGGCGAAATGGTCGGGGAAGGTTACCACCAACGGGCGGGTGACCCACATGCCGAAGTACCCGCGTTGCTAATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTa                                                                                                                                                                                    <  2:473404/185‑1 (MQ=255)
                                                                                                                  cACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                         <  1:535956/177‑1 (MQ=255)
                                                                                                                      cATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGGAAAAACCAAAAGAGCGACCCGCCATGTcaca                                                                                                                                                                                                                 >  2:726631/1‑150 (MQ=255)
                                                                                                                      aaTCCCAACCCGAAGGTCGGGTGGTTAATTGTCAAAGATGGCGAAATGGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACAGGCCGAAGTACACGGGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                     <  2:21653/238‑1 (MQ=255)
                                                                                                                       cctccAACCCAGAATGCCGGGCGCATATTTGACAAAAAGGGCGATAGTGTGGGTGAAGGTACCCCCAAGCGTGCGGGTGACCCAAATGCCGAAGTACCCTCTTTGCGTATGGCGGGTAAAAAAGCAAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                     <  2:1186559/236‑1 (MQ=255)
                                                                                                                         cccAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCAcccgg                                                                                                                                                                          >  1:473404/1‑182 (MQ=255)
                                                                                                                            aaCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGGGACCCCCTCTGTCACCCCCGAACCCCCTGAGCCTCATGGTCGTAGGccc                                                                                                                                                                               >  2:535956/1‑176 (MQ=255)
                                                                                                                                    aTGGCGGGGGCGTCATTGTAAAAGAGGGCAAAATTGTCGGTGAAGGTTCCCACCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCcatcat                                                                                                                                                                                           <  2:329917/158‑1 (MQ=255)
                                                                                                                                     tGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGc                                                                                                       <  1:1132118/239‑1 (MQ=255)
                                                                                                                                              cGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCGATGTCACACTCCAACCCTGTAGCCATCATGGTCGCACGc                                                                                                                                                                                 >  1:329917/1‑158 (MQ=255)
                                                                                                                                               gtcatttaaaaaGATGGCAAAATTGTCGGGAAAGGTTCCCACAAGCGTGGGGGTGAACCACATGCCGAAGTACACGGGTTGCGTATGGCGGGTAAAAAACACAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTgg                                                                                             <  2:190331/230‑1 (MQ=255)
                                                                                                                                                   ttGTCAAAGAGGGCAAAATTGTCGGTGAAGGTTCCCGCCAACGTGGGGGTGACCCACATGCCGAATTACAGGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGGGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                                         <  1:351326/239‑1 (MQ=255)
                                                                                                                                                             gggCGAAAGTGTCGGTGAAGGTTGCCCCAAAGGGGCGGGTGGACCACAGGCGGAGGTGCACGCGTGGCGTAGGGGGGGTAAAAAAGCCAAAGGTGCGGCCGCCTATGTCGCACTCGAACCCTGTAGCCATCATGGTCGTACGCCGCCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTAc                                                                                <  2:875036/237‑1 (MQ=255)
                                                                                                                                                              ggCGAAATTGTCGGTGAAGTTTCCCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGGGGGTGAAAAAGCAAAGGGTGCGACCGCCTATGTCACACTCGAACCCTg                                                                                                                                                                                                      <  2:1209252/122‑1 (MQ=255)
                                                                                                                                                              ggCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGCGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCCGTAGCCATCATTGTCGTACCCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGTTGCAAGACCCTAAc                                                                                                           >  1:61100/1‑210 (MQ=255)
                                                                                                                                                                   aaTCGTCGGTAAGGGTACCCACCAGCGAGGGGGACACCCACATGCCGAAGTACACGCGTTGGGTAGGGCGGTTGAAAAAGCCAAAGGGCCGACCGCCTATGTCACACTCGAACCCTGTAGCCAGCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGc                                                                         <  2:1385990/239‑1 (MQ=255)
                                                                                                                                                                        tCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGCACCCTGTAGCcatcat                                                                                                                                                                                           >  1:1209252/1‑122 (MQ=255)
                                                                                                                                                                                 gTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAAGGGCGACCCCCCATGTCACACTCGCACCCCGTAGCCATCATTGTCGTACGCCCCCCGGCTGTTACGCCCTCCTCCCCCCCGGCGTATCGCCCGGTGGTTCCTCGATGCCAAATCCTCACCCCCCGGTCCCCGGGGGGGGGCCTTTCCGTCTGCACCCGGcgcggcct                                                           >  2:728362/1‑232 (MQ=255)
                                                                                                                                                                                        ccaACGTGGGGGTAACCCACATGCCGAAGTCCCCGCTTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                     <  2:185748/174‑1 (MQ=255)
                                                                                                                                                                                           aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaa                                                                                                                                                                                                                                                <  2:123788/52‑1 (MQ=255)
                                                                                                                                                                                           aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGaa                                                                                                                                                                                                                                                >  1:123788/1‑52 (MQ=255)
                                                                                                                                                                                                  gggTGAACCCAGTGCAGAGGTACGCGCGTGGGGTATGGGGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGc                                            <  1:624068/238‑1 (MQ=255)
                                                                                                                                                                                                   ggTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTTTGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATTGCCGCTGGCGTAGCGCGCGCGGGTGCCTCGATGCAAGATCCTAAc                                                                                                           >  1:185748/1‑174 (MQ=255)
                                                                                                                                                                                                    ggaacccACTTGCAAAAGTGCCCGCGTTGCGTATGGCGGGGAAAAAACCCAAGGGTGGGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg                                              <  1:85687/229‑1 (MQ=255)
                                                                                                                                                                                                            caTGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCCTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGGGGTTGCCTCGATGCAAGATCCTAACCCGCCGGTCGCTGGGCGTGGGCTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCCACGCCtggtgatgat                                 >  2:1007476/1‑238 (MQ=255)
                                                                                                                                                                                                              tGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAAGGGCGACCGCCCCTGTCCCACCCCAACCCCGTAGCCATCATAGTCGTACGCCCCCGGGCCGTGACGCACTCAACGCCGCTGGCGTAGCGGGGGTGGGTGCCTCCGTTCAAAGACCCAAACCCCCGGGCCCCGGGGGGGGGGCTTTCCGTCCGCAACCAGCCGGGATTGACGTCAGCCCCCGcccgccgct                                    >  2:85687/1‑227 (MQ=255)
                                                                                                                                                                                                               gCCGAAGTACACGCGTTGCGTATGGGGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTa                                                                                 <  1:258109/188‑1 (MQ=255)
                                                                                                                                                                                                               gCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCCATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGa                              >  1:572370/1‑239 (MQ=255)
                                                                                                                                                                                                                    aGTCCACGCGTTGCGTATGGCGGGTGAAAAGGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCg                         <  1:181592/239‑1 (MQ=255)
                                                                                                                                                                                                                    aGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCg                         <  1:1002946/239‑1 (MQ=255)
                                                                                                                                                                                                                    aGTACACGCGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCg                         <  1:1007476/239‑1 (MQ=255)
                                                                                                                                                                                                                      tACACGCTTTGCTTATGGGGGGTGAAAAGCCAAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                     <  2:529503/145‑1 (MQ=255)
                                                                                                                                                                                                                         gcgcgTTGCGTAGGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa                    <  1:201838/238‑1 (MQ=255)
                                                                                                                                                                                                                         acGCGTTGCGTAGGGGGGGGAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTGCGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa                    <  1:584393/239‑1 (MQ=255)
                                                                                                                                                                                                                           gcgTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                   >  1:774941/1‑64 (MQ=255)
                                                                                                                                                                                                                           gcgTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                   <  2:774941/64‑1 (MQ=255)
                                                                                                                                                                                                                                gCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCCTCATGGTCGTACGCCACCCGGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAAc                                                                                                           >  1:529503/1‑145 (MQ=255)
                                                                                                                                                                                                                                gCGTATGGCGGGGGAAAAAACCAAAGGTGCGACCGCCTATGTCACACCCGAACCCCGTAGCCATCCTGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCCCTGGCGTAGCGCGCGGGGTTGCCTCGCTGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCACCAGGCTGGCATTGACGTCAGCCCCGGCCTGATGTTGCGTGAAGCCGAGCAATTTAATa             >  2:1362529/1‑239 (MQ=255)
                                                                                                                                                                                                                                 caaaTGGGGGGTAAAAAAGCAAAGGGTGGGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                             <  2:1015239/124‑1 (MQ=255)
                                                                                                                                                                                                                                 cGTATGGGGGGGGAAAAAGCCAAAGGTCCCCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTa                                                                                                                                                                                    <  2:992075/73‑1 (MQ=255)
                                                                                                                                                                                                                                 cGTAGGGCGGGGGAAAAAGCCAAAGGGGCGGCCGCCTATGTCCCACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGGGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGc                                            <  1:601281/207‑1 (MQ=255)
                                                                                                                                                                                                                                            gTGAAAAAGCCCAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGGCGTACGCCACCGTGCTGTGACGCACCCATCCCCCCCGGCGTAGCGCGCGTGGTTGCCCCGATGCAAGATCCTAACCCCCGTGTCGCTGGGGGGGGGCTTTCCCGTCTGCCACAGGCCGGCACTGTCGTCAGCCACGCCCCGATGAGGAGTGAAGCCCAGCAATTGAATCAAGGCATtctc  >  2:737738/1‑239 (MQ=255)
                                                                                                                                                                                                                                            gTGAAAAAGCCAAAGGGGCGACCGCCTATGTTCCACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGGACCCATCGCCGCTGGCGGAGCGCGCGGGGGTGCCTCGATGCAAGa                                                                                                                  >  1:1015239/1‑127 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 
TTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTC  >  NC_000913/433401‑433871

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