Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913433,6490TG39.7% 41.4 / 24.2 71G65G (GGT→GGGribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (5/36);  new base G (21/6);  total (29/42)
Fisher's exact test for biased strand distribution p-value = 6.71e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAG  >  NC_000913/433419‑433874
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                    
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)
                                   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)
                                                         ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGAc                                                                                                                                                                      >  1:1293321/1‑238 (MQ=255)
                                                          gggcgGGCGCTAAAGTGGGCCCAACGAGGACGTTTTCCCACGCATCCCAACCCGAATGTGGGGGGCGTAATTGAAAAAGATGGCGAAATTGTCGGTGAAGGTTCCCCCCACGGGGGGGGTGAACCAGAGGCAGAAGTACACGCGTGGGGTGTGGCGGGGGAAAAACCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCtgtg                                                                                                                                                                        <  2:1285911/233‑1 (MQ=255)
                                                          gcgcGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGGAAAAGCCCAAGGGGCGGACGCCCCTGTCCCACTCGCACCCCGTAGCCCTCCTGGGCGGACGCCACCGCGCCGTGACGc                                                                                                                                                                    >  2:54714/1‑239 (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)
                                                                                                                                                                                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)
                                                                                                                                                                                                                             aaaaaGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAg  <  1:54714/239‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                    
TTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAG  >  NC_000913/433419‑433874

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: