Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 432,873 T→G 25.5% G65G (GGT→GGG ribD → bifunctional diaminohydroxyphosphoribosylaminopyrimidine deaminase/5‑amino‑6‑(5‑phosphoribosylamino)uracil reductase RibD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007779432,8730TG25.5% 95.7 / 8.9 77G65G (GGT→GGGribDbifunctional diaminohydroxyphosphoribosylaminopyrimidine deaminase/5‑amino‑6‑(5‑phosphoribosylamino)uracil reductase RibD
Reads supporting (aligned to +/- strand):  ref base T (14/41);  new base G (17/4);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 1.59e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-01

ATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACC  >  NC_007779/432642‑433110
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                
gtttGGCGAAGAGATCGCGCGCCGGGAGGACTAAGCCGTGCAGGACTAGTATTACAGGGCGCGGGGTCTAAAGCTGGCGCACCGAGGACGTTTTACCACGAATCCCACCCGGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCTTAGGGCGGGTGAAAAAGCCAAAGGTGCGAc                                                                                                                                                                                                                                             <  2:703791/236‑1 (MQ=255)
 tttGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGATGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCACCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAAACAAAGAGGCGAccc                                                                                                                                                                                                                                           >  1:454616/1‑237 (MQ=255)
       gAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGAAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCACCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGt                                                                                                                                                                                                                                    <  1:774049/239‑1 (MQ=255)
                gcgcgcCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGa                                                                                                                                                                                                                          <  1:267310/239‑1 (MQ=255)
                         gaggaTAAAGCGGTGCAGGAGAGGTATTACAGGGGGGGGGGGTAAAGGCTGGCGAAACGAGGCGTTTTTCCCACGCAGCCCAACCCGAATGTCGGGTGGGTCATTGTCAAAGAGGGCGAAAGTGGCGGTGAAGGTTACCCCGAACGTGGGGGTGACCCAAAGGCGGAAGTACCCGCGTTGCGTATGGGGGGTAAAAACGCAAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                                  <  2:1005385/238‑1 (MQ=255)
                         gaggaCTAAGCCGTGCAGGACGAGTATTCCATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGCCTTTTTCCCACGCATCCCAACCCGAAGGTCGGGTGCGTCATTGTAAAAGATGGCAAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTAAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                                  <  1:1862815/238‑1 (MQ=255)
                         gaggaCTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                                  <  1:916215/238‑1 (MQ=255)
                         gaggaCTAAGCCGTGCAGGACGAGTATTACAGGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTGGGTGAAGTTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTGGCGTATGGCGGGTGAAAAACCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                                  <  1:601268/238‑1 (MQ=255)
                         gagaccgaaGGCGGGCAGGACGAGTATACCAGGGCACGGGCGCTAAGGCGGGCGCAACGAGTACGTTTTACCACGAATCCCAACCGGAATGTCGGGTGCGTCATTGACAAAGATGGGGAAATTGTGGGGGAAGGTTACCACCAACGGGCGGGTAACCCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCGAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                                  <  2:1245966/231‑1 (MQ=255)
                              cTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGGAAAAACCAAAAGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCcatc                                                                                                                                                                                                            >  1:836491/1‑239 (MQ=255)
                                 ggCCGGGCAGGACGGTTATTACATGGCGCGGGCGCAAAAGCTGGCGCAACGAGGACGTTTTCCCCGCCACCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAGGGTTACCACCAGCGTGCGGGTGAACCAAATGCCGAAGTACACGGGTTGCTTATGGCGGGTAAAAAAGCCAAAGGGGCGa                                                                                                                                                                                                                                              <  2:865828/202‑1 (MQ=255)
                                 aGCCGGGCGGGAGGAGTAGTAGCGGGCGCGGGCGCAAAGGGGGGCGTGACGGGGAGGTTTTACCCCGCATCCCAACCCGAAGGGCGGGTGCGTCCTGGTCAAAGAGGGCAAAATGGGCGGTGAGGTTTACCCCCACAGGGGGGGGAAACCACATGCGGAAGTACCCGCGTCGCGTAGGGCGGGGGAAAAGCCGAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCcatcat                                                                                                                                                                                                          <  2:551588/238‑1 (MQ=255)
                                   ccGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGGGCGACCCCCCATTTCACACTCCAACCCAGTTGcatcatcagg                                                                                                                                                                                                      >  1:246032/1‑235 (MQ=255)
                                    cGTTCAGGACGATTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACTTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTATCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAATTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTc                                                                                                                                                                                                      <  1:1681650/239‑1 (MQ=255)
                                    cGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                                        <  1:634884/237‑1 (MQ=255)
                                    cGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                                        <  1:1087905/237‑1 (MQ=255)
                                      tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAGCCAAAAGAGCGGCCGCCTATGTCAACATCGCACCCCGTAGCCCTCCTGGTCGt                                                                                                                                                                                                    >  1:504322/1‑239 (MQ=255)
                                      tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGCGACCGCCCATGTCACACCCGAACCCTGTAGCCATCCTGGGCGt                                                                                                                                                                                                    >  1:963136/1‑239 (MQ=255)
                                      tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAAACAAAAGGGCGACCGCCCATGTCACCCTCGCACCCTCTAGCCATCAATGTCg                                                                                                                                                                                                     >  1:2066055/1‑238 (MQ=255)
                                           gACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGGGCGACCGCCCATGt                                                                                                                                                                                                                                    >  1:865828/1‑203 (MQ=255)
                                             cGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                                              <  1:1321503/238‑1 (MQ=255)
                                               aGTATTACATGGCGCGGGCGCTAAAGCGGGCGCAACGAGAACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCg                                                                                                                                                                                                                           <  1:1366443/207‑1 (MQ=255)
                                                 tATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCg                                                                                                                                                                                                                                               <  1:419713/186‑1 (MQ=255)
                                                   ttACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAGGGGCGACCGCCCATGTCACACCCCAACCCTTTAGCCCTCCAGGTCGTACGCCCCCgcggc                                                                                                                                                                                       >  1:1400691/1‑235 (MQ=255)
                                                     aCATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGTTGCGACCGCCTATGTCACACTCGa                                                                                                                                                                                                                          <  1:1745995/202‑1 (MQ=255)
                                                     aCATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCAAAGGGGCGACCGCCCATGTCACACTCCAACCCTGGAGCCATCATGGTCGGACGCCACCGTGCcgt                                                                                                                                                                                     >  1:1040493/1‑239 (MQ=255)
                                                         ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGTGGGACAGCCCACGTCACACTCGAACCCCGTAGCCATCATAGGCGCACGCCACCCTGCTCTGAc                                                                                                                                                                                  >  1:2089800/1‑238 (MQ=255)
                                                          gcgcGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAAGTGCGACCGCCTATGTCACACCCGAACCCTGTAGc                                                                                                                                                                                                                >  1:1447391/1‑207 (MQ=255)
                                                           cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAACCAAAAGTGCGACCCCCCATTTCACACCCCAACCCCGTAACCAACATTGTCGTAAGCCCCCGGGCGTTGACGc                                                                                                                                                                                >  1:1912775/1‑238 (MQ=255)
                                                             cgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCCAAGGGGCGACCGCCTGTGGCACACTCGAACCCTGTAACCATCATGGGCGTACg                                                                                                                                                                                                 >  1:1482074/1‑219 (MQ=255)
                                                                 cgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATGGTCGGTGAGGGTTCCCCCCAGCGTGCGGGTAAACCACATGCCGAAGTACGCGCTTTGCGTATGGCGGGTGAAAAGGCCAAGGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCAt                                                                                                                                                                          <  1:1597777/238‑1 (MQ=255)
                                                                                        cGTTTTCCCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGTTTACCACCAACGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACg                                                                                                                                                                                 <  1:316298/208‑1 (MQ=255)
                                                                                             tACCACGCATCCCAACCCGAATGTCGGGGGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGGGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                                             <  1:652968/239‑1 (MQ=255)
                                                                                                cccccgcccccACCCGAGTGTGGGGGGCTTCTTGGACAAAGAGGGGAAATGTGTCGGGAAGGGTTCCCCCCAGGGTGGGGGTGACCCACAGGACGAAGCACGCGTGTTGCGTATGGCGGGTGAAAAACCCAAAGGTCGGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGa                                                                                                                                          <  2:1400691/232‑1 (MQ=255)
                                                                                                     aTCCCAACCCGAAGGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCACCGTGCGGGTGAACCACATGCCGAAGTACACGCGTGGCGTATGGCGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                                      <  1:856613/238‑1 (MQ=255)
                                                                                                       cccAACCCGAATGTCGGGTGCGTCATTGTCAAAGAGGGCGAAATTGTCGGTGAAGGTTCCCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGaa                                                                                                                                                                                                                         <  1:1588335/153‑1 (MQ=255)
                                                                                                           aCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGTATGGGGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCt                                                                                                                              <  1:1645709/239‑1 (MQ=255)
                                                                                                                           cGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGTCGGGTGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGGACGCCACCGGGCTGTGACGCACTCATCGCCGCTGGCGGAGCGCGCGTGGGTTCCCCGATGCAAGATCCTAACCCGCAGGTCGCTg                                                                                                              >  1:1631617/1‑239 (MQ=255)
                                                                                                                           cGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGGGCGACCGCCTATGTCACACTCCAACCCTGTAGCCCTCATGGGCGTTCGCCCCCGTGCtgtg                                                                                                                                                                                    >  1:2010111/1‑170 (MQ=255)
                                                                                                                              cATTGTAAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGc                                                                                                           <  1:290671/239‑1 (MQ=255)
                                                                                                                                  gTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGTGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTgg                                                                                                       >  1:1668142/1‑239 (MQ=255)
                                                                                                                                   tCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCACCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCtgtg                                                                                                                                                                                    <  1:309157/162‑1 (MQ=255)
                                                                                                                                      aaGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGCATGGCGGGGGGAAAAGAAAAAAGGGGGGGCGACCATGTTACACTAGAACACTGAAGCAAACAAGGGCGAAAGGCACCGGGCACGGACGTACTCAACACCGCCGGCGTAGAGGGAGCGGCAGCCTCGATGCAACAGTGCAAACCCACGCGTGCTGGGCGTgggcgt                                                                                                   >  2:644650/1‑235 (MQ=255)
                                                                                                                                          tGGCAAAATTGTCGGTGAAGGTTACCACCACCGTGCGGGTGAACCCCATGCCGAAGTCCCCGCGTTGCGTATGGCGGGTGAAAAAGCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTa                                                                                                 <  1:734669/237‑1 (MQ=255)
                                                                                                                                                       cGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCTTGTAGCCATCATGGGCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAg                                                                                                                                    >  1:1138179/1‑191 (MQ=255)
                                                                                                                                                         gTGAAGGTTACCACCATCGTGCGGGTGAACCCCATGCCGAAGTCCCCGCGTTTGTTTTGGCGGTTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGc                                                                                  <  1:644650/238‑1 (MQ=255)
                                                                                                                                                         gTGAAGGTTACCACCAGCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGc                                                                                  <  1:921523/238‑1 (MQ=255)
                                                                                                                                                                    caccaACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGACAAAGGGGCGACCGCCCATGTCACACTCGAAACCTGTAGCCATCATGGGCGTGCGCCCCCGGGCTGTGAAGCACTCATCGCCGCTGGCGTAGCGCGCGGGGGTGCCTCGATGCaaaa                                                                                                                                   >  1:1560898/1‑177 (MQ=255)
                                                                                                                                                                         aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGGAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAAACCTGTAGCCATCATGGGCGTACGCCAcc                                                                                                                                                                                            >  1:2023099/1‑117 (MQ=255)
                                                                                                                                                                         aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGcc                                                                 >  1:1522324/1‑239 (MQ=255)
                                                                                                                                                                         aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGGGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGc                                                                  >  1:1822983/1‑238 (MQ=255)
                                                                                                                                                                          cGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATcgc                                                                                                                                                                       <  1:809955/137‑1 (MQ=255)
                                                                                                                                                                             gCGGGTAAACCACATGCCGAAGCCCACGCGTTGCGTAGGGCGGGTAAAAAACCAAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg                                                             <  1:1185880/238‑1 (MQ=255)
                                                                                                                                                                             aagggTGAACCACATGCCGAAGTACACGCGTTGCGTAGGGGGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg                                                             <  1:1428077/236‑1 (MQ=255)
                                                                                                                                                                              cGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAACCAAAAGGGCGACCGCCCATGTCACACACCCACCCCGTAGGCAACATCGTCGGACGCCACCCCGCTGTGTCGCCCCCCACGCCGCCGGCCGAGCCCCCGCGCGTGCCTCCCTGCAAAGTCCCAAACCGCCCGCCGCCGCGCGGTGGCTTTTCCGGCTGCCACCGGCTTGCATTGACGTCCACCACgg                                                            >  2:952840/1‑238 (MQ=255)
                                                                                                                                                                                    accccaaTGCCGAAGTCCGCGGGTTGCGTAGGGCGGGTAAAAAGGCCAAGGGTGCGGCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCtgatg                                                     <  1:952840/233‑1 (MQ=255)
                                                                                                                                                                                    aaCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAACCAAAAGGGGGACCCCCCATGTCACACCCCGAACCCGTAGCACTCATGGGCGTACGCCACCCTGCTGTGAAGCACCCCTCCCCGCCCGCGTa                                                                                                                                                             >  2:809955/1‑137 (MQ=255)
                                                                                                                                                                                         aatgcgaaaGTACTCGCGTGGCGTAGGGCGGGTAAAAAACCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAg                                                <  1:2042995/232‑1 (MQ=255)
                                                                                                                                                                                           tGCCGAAGTACCCGTGTTGCGTAGGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                                                         <  1:1992635/181‑1 (MQ=255)
                                                                                                                                                                                           tGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGGGACCGCCTATGTCACACTCGAACCCTGTAACCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCCCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGGGGACTTTACCGTCTGCAACAGGCTGGCATTGa                                                                         >  1:1472278/1‑213 (MQ=255)
                                                                                                                                                                                           tGCCGAAGTACACGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGAc                                                                        >  1:182395/1‑214 (MQ=255)
                                                                                                                                                                                                    ccacGCGTTGCGTATGGGGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg                                      <  1:768768/237‑1 (MQ=255)
                                                                                                                                                                                                     cacGCGTTGCGTATGGCGGGTGAAAAAACCAAAAGTGCGACCGCCCCTGTTACACACGAACCCCGTAGACATCATGGTCGGACGCCACCGTGCTGTGACGCACTCATCGCCCCCGGCGGAGCGCGCGGGGGTGCCTCGATGCAACAACCTAACCCGCCGGGCGCTGGGGGTGGGCTTTTcc                                                                                               >  2:1992635/1‑181 (MQ=255)
                                                                                                                                                                                                     cacGCGTTGCGTATGGCGGGGGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCCTTGACGTCAGCCACGGCCtgatgatga                                                 >  1:275320/1‑226 (MQ=255)
                                                                                                                                                                                                      ccgcgTTGCGTAGGGGGGGTGAAAAAGCAAAGGGGGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa                                   <  1:1950715/238‑1 (MQ=255)
                                                                                                                                                                                                      acggtttgtttATGGGGGTAGAAAAACCAAAGGGTGCGCCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa                                   <  1:1924405/229‑1 (MQ=255)
                                                                                                                                                                                                              cGTATGGGGGGTAAAAAAGCCAAGGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc                                                                                                                <  1:2112139/155‑1 (MQ=255)
                                                                                                                                                                                                               atATGGGGGGTAAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATaa                           <  1:436928/237‑1 (MQ=255)
                                                                                                                                                                                                                 aTGGCGGGTGAAAAAGCCAAGGGGGCGACCGCCTATGTCACACTCGATCCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAg                         <  1:870471/238‑1 (MQ=255)
                                                                                                                                                                                                                    gCGGGTGAAAAAGCCAAAGGGGGGACCCCCCCTGTCACACTCGAACCCTGTAGGCATCATGGTCGTACGCCACCGTGCTGTGACGCCCTCATTCCCCCTCGCCTAGCGCCCGGGGGTGGCCCGATGCAAAATCCTTACCCCCAGGGCCGGCGGCGTGGGCGTTCCCGTCTGCCACACGCTGGGCTTGACGTCCGCCCCGGGCCTGCTGTGGGGGAAAGCCGAGCAAttaacaaaagtgt                      >  2:1263858/1‑228 (MQ=255)
                                                                                                                                                                                                                        gTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCaa                                                                                       >  1:869127/1‑170 (MQ=255)
                                                                                                                                                                                                                           aaaaaGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAg              <  1:1560960/239‑1 (MQ=255)
                                                                                                                                                                                                                           aaaaaGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAg              <  1:1458470/239‑1 (MQ=255)
                                                                                                                                                                                                                           aaaaaCCAAAAGGTGCAACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAg              <  1:1247795/239‑1 (MQ=255)
                                                                                                                                                                                                                                  cAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATgc       >  1:1099479/1‑239 (MQ=255)
                                                                                                                                                                                                                                   aaaGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATgc       >  1:1064769/1‑238 (MQ=255)
                                                                                                                                                                                                                                       gTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCa    <  1:837831/237‑1 (MQ=255)
                                                                                                                                                                                                                                        tGCGGCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCAcc  <  1:1374012/238‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                
ATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTGAAAAAGCCAAAGGTGCGACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCACC  >  NC_007779/432642‑433110

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