Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913433,6490TG33.9% 56.5 / 14.2 61G65G (GGT→GGGribDfused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (8/31);  new base G (17/3);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 2.91e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGC  >  NC_000913/433415‑433881
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
cgAATTTGGGAAAGAACCCGCGCCCCTGGAGGCCTAAGCCGTGCAGGACGAGTTTTACAGGGCGCGGGCCCTAAAGCTGGCGCACCGAGACCGTTTTCCCCCGCACCCCAACCCAAATGTCGGGTGCGTCATTGTAAAAGAGGGCAAAATTGTGGGTAAAGGTTCCCCCCACCGGGCGGTTGAACCACATGCCGAAGTACACGCGTTGCGTATGGGGGGTGAAAAAGCCAAAGGTGGGa                                                                                                                                                                                                                                       <  2:673891/238‑1 (MQ=255)
   aTTTGGGGAAGAGATCGCGCGCGTGGAGGACTAAGCGGTGAAGGACGAGTATTAAAGGGCGCGGGCGCAAAAGCGGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTAATTGTCAAAGAGGGCAAAATTGTCGGTGAGGGTTCCCACCAGCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGAcc                                                                                                                                                                                                                                     <  2:207713/238‑1 (MQ=255)
   aTTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGAcc                                                                                                                                                                                                                                     <  1:242430/238‑1 (MQ=255)
                  cgcgcgCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACAAGCCGAAGTACACGCGTTGCGTATGGCGGCGGGAAAAAGCCAAACGGGGGCCGCCCTCTGTCacacca                                                                                                                                                                                                                      >  2:115133/1‑237 (MQ=255)
                          tGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGCAAAAGCCCAAGGTGCGACCCCCTCTTTCccaccccaaacccgt                                                                                                                                                                                                             >  1:30787/1‑228 (MQ=255)
                            gaggaCTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                            <  1:771060/238‑1 (MQ=255)
                               gaCTAGGCCGTGCAGGACGAGTATTCCATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAAcc                                                                                                                                                                                                                 <  1:423233/230‑1 (MQ=255)
                                 cTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGGCAAAGGTGGGACCCCCTATGTTACAACCGAACCCTGTAGCcact                                                                                                                                                                                                      >  1:328639/1‑237 (MQ=255)
                                   aaGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGGAAAAACCAAAGGTGCGAGCGCCTATGTCACACTCGAACCCTCTAGCcat                                                                                                                                                                                                       >  1:647077/1‑236 (MQ=255)
                                      ccGTGCAGGACGATTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                            <  1:191359/228‑1 (MQ=255)
                                                cGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTGCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTAGGGGGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                                        <  1:99582/238‑1 (MQ=255)
                                                cGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAGAAAACCGAAAGGGGGGCCGCCTACGTCACAATGCACCCCAGTATCCACCATgg                                                                                                                                                                                                  >  2:191359/1‑228 (MQ=255)
                                                   gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAGGGGGGACCGCCCATGTTACACTCGAACCCCGTAGCCATCATGGGCGGACGcccc                                                                                                                                                                                       >  2:904541/1‑234 (MQ=255)
                                                   gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGGGCGACCGCCTATGTCACACTCCAACCCTGTAGCCATCATGGTCGTACCCCACCGt                                                                                                                                                                                    >  1:746851/1‑239 (MQ=255)
                                                     aTTCCATGGCGGGGGCGCAAAGGCGGGCGCAGGGAGACCGTTTTACCACGCATCCCAACCCGAGTGTCGGGTGCGGAATTGTAAAAGAGGGCGAAAGTGGCGGGGAAGGTTGCCCCCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTAAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                                   <  2:328639/238‑1 (MQ=255)
                                                              cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAAGGGCGACCGCCTATGTCACACTCCAACCCTGTAGCCCTCGTGGTCGTAccccc                                                                                                                                                                                        >  1:246869/1‑223 (MQ=255)
                                                                cgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAGCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAAGGGAGAACGCCTATGTCACACTCCAAACCCGTAGACATCAAGGGCGTAAGCCACCGGGCTGTGACGCAc                                                                                                                                                                        >  1:443496/1‑238 (MQ=255)
                                                                     gcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGt                                                                                                                                                                                                                             <  1:254217/180‑1 (MQ=255)
                                                                      cTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAAGGGCGACCCCCTATGTCCCACTCGAACCCTGTAACCATCATGGTCGTACGCCACCCGGCTGTGACGCCCACATcg                                                                                                                                                                  >  1:328916/1‑238 (MQ=255)
                                                                               gcgcAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAAACAAAAGGGGGACCGCCTCTGTTACACCCGaa                                                                                                                                                                                                                   >  2:254217/1‑180 (MQ=255)
                                                                                     cGAGGACGTTTTACCACGCATCCAACCCCGAATGTCGGGTGCGGCATGGAAAAAGAGGGCAAAATTGTCGGTGAAGGTTACCCCCAACGGGCGGGTAAACCACAGGCCGAAGTACCCGCGTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                                   <  2:427909/206‑1 (MQ=255)
                                                                                          aCGTTTTACCACGCACCCCAACCCGAATGTCGGGTGCGTCATTTTAAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACAGGCCGAAGTACACGCTTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCt                                                                                                                                                                                                                                 <  2:898331/155‑1 (MQ=255)
                                                                                            gTTTTACCACCAATCCCAACCCGAAGGACGGGTGCGTCATTGTAAAAGAGGGCAAAATTGGCGGGAAGGTTTACCCCCAAGGTGCGGGTGAACCACATGCCGAAGTACCCGCGTGGCGTAGGGCGGGTGAAAAAGCAAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTg                                                                                                                                                                                   <  2:158309/199‑1 (MQ=255)
                                                                                                tACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                                       <  1:241038/239‑1 (MQ=255)
                                                                                                    aCGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACt                                                                                                                                                                                                                       >  1:898331/1‑155 (MQ=255)
                                                                                                      gCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCa                                                                                                                                                                         >  1:158309/1‑199 (MQ=255)
                                                                                                       cATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                                <  1:904541/239‑1 (MQ=255)
                                                                                                         gcccAACCCGAATGTCGGGGGCGTCATTGTCAAAGAGGGCGAAATTGTCGGTGAGGGTTACCGCCAACGTGCGGGTGACCCACATGCGGAAGTACACGCGTGGCGTATGGCGGGTGAAAAACCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                                        <  2:211909/180‑1 (MQ=255)
                                                                                                              aCCCGAAGGGCGGGTGGGTAAGTGGAAAAGATGGCGAAATTGTGTGGGAAGGTTACCACGAACGGGCGGGTAAACCCCGAGCCGAAGTACACGCGTTGCGTAGGGCGGGTGAAAAAGCAAAGGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCt                                                                                                                         <  2:918443/239‑1 (MQ=255)
                                                                                                                  gAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCCCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAAc                                                                                                                      >  1:879487/1‑238 (MQ=255)
                                                                                                                   aaTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCCATGTCACACTCGAACCCTGTAGCCATCATGGTTGTACGCCACCCTGccttg                                                                                                                                                                              >  1:211909/1‑177 (MQ=255)
                                                                                                                                 cATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGc                                                                                                      <  1:483489/239‑1 (MQ=255)
                                                                                                                                 atttGTAAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGGGCGGGTGACCCACATGCCGAAGTACACGCTTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTca                                                                                                                                                                                                                           <  1:426580/120‑1 (MQ=255)
                                                                                                                                   ttgtaaagaaTGGCAAATTTGTCGGTGAAGGTTACCGCCAACGTGCGGGTGAACCACATGCCAAAGTACGCGCTTTACGTATGGCGGGTAAAAAAGCAAAAGTTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                                                   <  2:443496/230‑1 (MQ=255)
                                                                                                                                   ttGACAAAGTGGGGAAATTGGGGGGGAAAGTTTCCCCCCAGCGTGCGGGGAAACCACATGCCAAAGTCCCCGCGTTGCGTAGGGGGGGGAAAAAAGCCAAAGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                                                   <  2:950828/239‑1 (MQ=255)
                                                                                                                                   gtgtgAAAGATGGCGAAATGGTCGGTGAAGGTTACCCCCAACGGGGGGGGAACCCACAGGCGGAAGTACCCGGGTTGCGTATGGCGGGTAAAAAAGCCAAAGTTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                                                   <  2:394412/238‑1 (MQ=255)
                                                                                                                                           gATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGCGCGAACGCCCATGTTCCCCTCGAAcc                                                                                                                                                                                                                 >  2:426580/1‑122 (MQ=255)
                                                                                                                                             tGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTa                                                                                           <  1:514804/237‑1 (MQ=255)
                                                                                                                                              ggCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCAAACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGGTGTGACGCACTCCTCCCCCCTGGCGGAGCGCGCGTGGTTGCCTCGATGGAAGATCCTAAACCGCAGGGCGCTGGGCGCTGGCCtttcc                                                                                          >  2:257899/1‑236 (MQ=255)
                                                                                                                                              ggCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGTGCGACCGCCTATGTTCCACTCGAACCCTGTAGCCCTCATGGTGCTAAGCCACCCGGCTGTGACCCACTCATCCCCCCTGGCGTAGGCCGCGTGGTTGCCCCGCGTCAAAATCCTAACCCCCCGGGCGCCGGGCGCgggcctttac                                                                                         >  2:1017029/1‑230 (MQ=255)
                                                                                                                                               gCGAAATGGTCGGTGAAGGTTCCCACCAACGGGCGGGTGAACCACATGCCGAAGTACACGCGTTGCTAATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTATGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTa                                                                                           <  1:868472/235‑1 (MQ=255)
                                                                                                                                                cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAAACAAAGGGGTGATCCCCTATGTTACACCCCAAGCCTGTAACCATCCCGGTCCCATGCCCACGGGCTGTGTCGCACTCATCGCCGCTGGCGCAGCCCGCGGGGTTGCCTCCCTGCTAGATCCCAAACCGCGGGTCGACGGGccggggcgtt                                                                                             >  2:776765/1‑224 (MQ=255)
                                                                                                                                                         cGGTGAGGGTTCCCCCCAAGGGGGGGGTAACCCACACGCCGAAGTACCCGCGTTGCGTGTGGGGGGTGAAAAAGCCAAAGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAgg                                                                             <  2:746851/239‑1 (MQ=255)
                                                                                                                                                          ggTGAAGGTTCCCCCAAACGTGGGGGTAAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAACCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAgg                                                                             <  1:888164/238‑1 (MQ=255)
                                                                                                                                                          ggTAAAGGTTACCACCAACGTGGGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGGGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAgg                                                                             <  1:645425/238‑1 (MQ=255)
                                                                                                                                                            gtaaaGGTTACCCCCAACGTGCGGGTGACCCCCATGCCGAGTTACACGCGTTGCATATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGc                                                                            <  1:776765/235‑1 (MQ=255)
                                                                                                                                                                      accaccAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGGAAAAGCCAAAGGTGCGACCGCCTCTGTCACACTCGAACCCCGTAGCCATCCTGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCCGGTCGCTGGGCGGGGGCTTTACCCTCTGCCACAGGCCGGGAtttg                                                                   >  2:239071/1‑236 (MQ=255)
                                                                                                                                                                      accaccAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGGGCGACCGCCCATTTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCCCCGCCGGCGTAGCGCGCGGGGGTGCCTCGATGCAAGAACCTAACCCCCAGGGCGCCGGGGGGGGGCGTTTCCGTCGGCAAAAGGCTGGc                                                                        >  2:969789/1‑233 (MQ=255)
                                                                                                                                                                              cGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                          <  1:969789/239‑1 (MQ=255)
                                                                                                                                                                               gTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                                                   <  1:560627/197‑1 (MQ=255)
                                                                                                                                                                                 gCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg                                                        <  1:345694/238‑1 (MQ=255)
                                                                                                                                                                                   ggggaaaCCACGTGCCAAAGTACGCGGGGTGCGTTGGGGGATAAAAAAAGCCAAAGGGGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGc                                                      <  2:328916/233‑1 (MQ=255)
                                                                                                                                                                                     gcgAACCACATGCCAAAGGACACGGTTGTGGATGGGGGCGAAAAAAAGCCAAGGGGGCGCCCGCCTAGGTCACACGCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg                                                        <  2:647077/232‑1 (MQ=255)
                                                                                                                                                                                          ccACGTGCCGAAGTACACGCGTGGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                                                   <  2:108322/186‑1 (MQ=255)
                                                                                                                                                                                          ccACATGCCGAAGTACACGAGTGGCGGATGGGGGGTGAAAAGCCAAAAGGTGCCGCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                          <  1:481890/227‑1 (MQ=255)
                                                                                                                                                                                                  cGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGGGCGACCCCCCATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGGGCTGTGACGCACCCCTCCCCGCTGGGGTAGCGCGCGGGGGTGCCCCCGTGCAAGAACCTAAACCCCCGGGCGCCGGGGGGGGGCTTTACCGTCTGCAACAAGCTGGGATTGGCGGCCGCCACGGCCTGTTGTTGAGTGaa                                       >  2:909852/1‑238 (MQ=255)
                                                                                                                                                                                                    aaGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATTTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTAcc                                                                                         >  1:108322/1‑186 (MQ=255)
                                                                                                                                                                                                    aaGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGAACGCCCATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCCCTGGCGGAGCGCGCGGGGGTGCCTCGCTGCAAGGTCCTAACCCCCAGGTCGCTGGGCGGGGGCTTTCCCGTCTGCAACAGGCCGGCCTTGACGGCAGCCACGGcccggtt                                                >  2:481890/1‑222 (MQ=255)
                                                                                                                                                                                                                  cGTATGGGGGGTAAAAAACCAAAGGGCGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATa                       <  2:879487/238‑1 (MQ=255)
                                                                                                                                                                                                                       ggCGGGGGAAAAAACCAAAGGGGCGACCCCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCCTCGCCGCTGGCGGAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCTGGGCGGGGACTTTACCGTCTGCAACAGGATGGGATTGACGTCAGCCCCGGCCTGGTGGTGGGGGAAGCCGGGCAATTTAATAAAgg                   >  2:960254/1‑237 (MQ=255)
                                                                                                                                                                                                                                        cAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATgc  >  1:418446/1‑239 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                         
AGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGC  >  NC_000913/433415‑433881

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: