Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,895,402 A→C 26.2% N70K (AAT→AAG)  yiaN ← 2,3‑diketo‑L‑gulonate transporter large permease YiaN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077793,895,4020AC26.2% 57.6 / 3.7 44N70K (AAT→AAG) yiaN2,3‑diketo‑L‑gulonate transporter large permease YiaN
Reads supporting (aligned to +/- strand):  ref base A (17/12);  new base C (10/2);  total (28/16)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGGCGCGATAATGCCGCCAGAGGCAATCAGCCCCGCCGCCCGGTTGACCGGATAGTTGGCGCTGCGCATCATCGGCACCAGTAGAGCGGCAACGGCGGCGGTATCCGCAACGGCAGAACCAGAAAGGCTCGCCATAATCATAGCTGCCAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGTCAACAATGCGTTTT‑G‑ACAAGCCGCCCGCATTCATGATTTCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGCCATGATTAAGCTCCTTGCGGTGGCGACGTTA  >  NC_007779/3895166‑3895639
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
tttGGCGCGATAATGCCGCAAGAGGCAATCACCCCCGCCGCCCGGTTGACCGGATAGTTGGCGCTGCGCATCATCGGCACCAGTAGAGCGGCAACGGCGGCGGTATCCGCAACGGCAGACCCAGAAAGGCTCGCCATAATCATAGCTGCCAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGTCAACAATGCGTTTT‑G‑ACAAGCCGCCCGCAtt                                                                                                                                                                                                                                             <  1:694926/238‑1 (MQ=255)
   ggCGCGATAATGCCGCCAGAGGCAATCAGCCCCGCCGCCCGGTTGACCGGATAGTTGGCGCTGCGCATCATCGGCACCAGTAGAGCGGCAACGGCGGCGGTATCCGCAACGGCAGAACCAGAAAGGCTCGCCATAATCATAGCTGCCAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGGCAACAAAGCGGTTT‑T‑GCAAAGCGCCCCCCTTCAt                                                                                                                                                                                                                                          >  1:189476/1‑239 (MQ=255)
                                                         ttGGCGCTGCGCATCATCGGCACCAGTAGAGCGGCAACGGCGGCGGTATCCGCAACGGCAGAACCAGAAAGGCTCGCCATAATCATAGCTGCCAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGGCCAAAATGCGCTTT‑G‑GCAAGCCGCCCGCATTCCTGGTTTTACCCGCCAGCCCAAAAAACGGCATCGCCAGCCGGGGGAAGCTATCgg                                                                                                                                                                                     >  1:2080353/1‑238 (MQ=255)
                                                                 gcgcATCATCGGCACCAGTAGAGCGGCAACGGCGGCGGTATCCGCAACGGCAGAACCAGAAAGGCTCGCCATAATCATAGCTGCCAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGTCAACAATGCGGTTT‑G‑ACAAGCCGCCCCCATTCATGATTTCACCCCCCAGCACAAAAAACGGAATCGCCAGCAGGGGGAAGCCATCGGCTCCGCTc                                                                                                                                                                             >  1:1288037/1‑238 (MQ=255)
                                                                                                             cagcggaaccaCCAGAACGCCGCGACATAAACATACCCCCAAGCGCGCGGAGGCGCCCCCGCCCGCCTGGTTTATCCCCCACCATTTTCATCGGCAGGTAAAAAATGCGTTTG‑A‑AAACGCCCCCCGCATTCATGATTCCCCCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTAGTGCCATGATCTGGACATCAAACATGTccagccag                                                                                                                                <  2:649960/229‑1 (MQ=255)
                                                                                                                   cAGAACCAGAAAGGCTCGCCATAATCATAGCTGCCAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGGCAACAATGCGTTTT‑G‑ACAAGCCGCCCCCCTTCATGATTTCACCCCCCACCCCCAAGAACGGAAACCCCAGCCGGGAGAAGGCATCCGCCCCCTTCACCACCGTTTGTGCCATGATTCGGAAAACAAAAAAGTccagccagcaaact                                                                                                                          >  2:1167875/1‑233 (MQ=255)
                                                                                                                   cAGAACCAGAAAGGCTCGCCATAATCATAGCTGCCAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGGCAACAAAGCGGTTT‑G‑ACAAGCCCCCCGCAATCATGATTTTACCCGCCCGCACAAAGAAAGGAAACGCCACCACGGAGAAGACAACGCCCCCCTTCACCCACGGTTTTGCCCTGATCTCGACATCAAACATGTCCCCCCAGAAAAta                                                                                                                          >  2:1078201/1‑238 (MQ=255)
                                                                                                                            aaaGGCTCGCCATAATCATAGCTGCCAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGTCAACAATGCGTTTT‑G‑ACACGCCGCCCGCATTCATGATTTCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAAcgccgc                                                                                                                  >  1:881838/1‑238 (MQ=255)
                                                                                                                                  tCGCCATAATCATAGCTGCCAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGGCAACAAAGCGTTTT‑G‑ACAAACCGCCCCCCTTCATGATTTCACCCGCCAGCACAAAAAACGGAAACCCCAGCCGGGAGAAGCTATACGCTCCCGTCCCCAGCGGTTGTGCCATGATCTGGACATCAAAAATGTCCACCCAGAACATTAACGCCCCCCCCca                                                                                                            >  2:1533062/1‑238 (MQ=255)
                                                                                                                                   cGCCATAATCATAGCTGCCAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGGCAACAAAGCGGTTT‑G‑AAAAACCCCCCGCCTTCATGATTTTTCCCGCCAGCACAAAAAAAGGAAAAGCCAGCAGGGGGAAACAAACTGCCCCCGTCACCCACGGTTTTGCCCTTGTCTGGAAATCCAAAATGGCCAACCCAAACAATAACGCCGccccccccc                                                                                                          >  2:2231349/1‑234 (MQ=255)
                                                                                                                                              tAGCTGCCAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGCCCACCACGCGGTTT‑G‑ACAAGCAGCCCCCGTTCATCATTTCACCCCCCCGCACACCGAACGGAACCGCCCGCAAGGAGAAACACTCGCCCTCCCTCACCCGCCCTTGTGCCCCGAGCAGGATATCAAACCACTTCACGCAGCACATTCAAGCCGCCCCCGCCAACAACAGGGc                                                                                                >  2:2020002/1‑238 (MQ=255)
                                                                                                                                               cgcccccaccccGCCGGCGTACCCCACCCCGCGTGTTTTATGCCCCCCCAGTTTCAGCGGCAGGTCAACAATGCGTTGT‑G‑ACAAGCCCCCCGCATTCATGATTTCCCCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGccc                                                                                              <  1:1215909/228‑1 (MQ=255)
                                                                                                                                               aGCTGCCAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGTCAACAAAGCGTTTT‑G‑ACAAGCCGCCCGCATTCATGATTTCACCCCCCAGCACAAAGAAAGGAATAGCCAGCAGGGAGAAGCTATCGGCTCCGTTCAACAGAGTTTGGGCCATGAGCTGGa                                                                                                                                                    >  1:301597/1‑185 (MQ=255)
                                                                                                                                                gCTGCCAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGTCAACCAAGCGTTTT‑T‑ACAAACCGCCCCCCTTCATGATTTTACCCCCCAGCACAAAAAAAGGAATCGCCAGCAGGGGGAAGCTATCGGCCCCCCTCACCAGCGTTTTTGCCATGATCTGGACATCAAAAATTTCCAGCCAGAACAATAACTCCCCCCCCCACAACAACAATGCCCa                                                                                             >  2:1695044/1‑239 (MQ=255)
                                                                                                                                                gCTGCCAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGGCCACCAAGCGTTTT‑T‑ACCAACCGCCCCCCTTCCTGATTTTACCCCCCAGCACAAAAAAAGGAATCCCCAGCAGGGGGAAGCTATCGGCCCCGCTCACCACCGGTTTTGCCCAGATATGGGCACCAAAACCATCCAGCCCAAACATTTACGCCCCCCCGCacaccaaacatgcccc                                                                                             >  2:2019864/1‑230 (MQ=255)
                                                                                                                                                     cAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGTCAACAAAGCGTTTT‑T‑ACAAACCCCCCGCATTCATGATTTCACCCGCCAGCACAAAAAAAGGAAACGCCAGCCGGGAGAAACTATCGGCCCCCTTCCCCAGCGCTTGTGGCCTGATCTGGACATCAAAAAAGGCCCGCCAGAACATTAAAGACCCCCCCCACAACAACAAAGCCCAAGCaa                                                                                        >  2:1195938/1‑239 (MQ=255)
                                                                                                                                                                    cccccccccGCCTGGTTTACGCCCCCCCAGTTTCACCGGCAGGACAACAATGCGTTTT‑G‑ACAACCCGCCCGCCTTCATGATTCCCCCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAg                                                                         <  1:1681601/233‑1 (MQ=255)
                                                                                                                                                                       cagcccGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGTCCACAATGCGGTTT‑T‑TCAAACCCCCCGCCTTCATTATTTTACCCGCCCCCACAAAAAACGGAAAAGCCCGCACGGAGAAACTATCCGCCCCGTTCACCAACGTTTGTGCCACGATCGTGACAACAAACATGTCCAGCCAAAACATTAACGCCGCCCCCCCCAACAACAAGGCCCCGGGAAAAAGCAAAcccaaaacca                                                                      >  2:202956/1‑231 (MQ=255)
                                                                                                                                                                         gcccGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGTCAACAAAGCGCTTT‑G‑ACAAAACGCCCGCATTCATGATTTCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCCATCGGCTCCGTTCACCAGCGGTTGTGCCATGATCTGGACATCAAAAATGTCCAGCCAGAACATTAAAGCCGCCCCGCACAACAACAGAGCCCACGCAATAGGCAAACCGaccacaat                                                                     >  1:844289/1‑231 (MQ=255)
                                                                                                                                                                          cccGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGGCAACAAAGCGTTTT‑T‑ACAAACCGCCCGCCTTCATGATTTCACCCGCCAGCACAAAAAACGGAATCGCCAGCAGGGAGAAACTATAGGCTCCCGTCACCAGCGGTTGGGCCATGATCTGGGCATCAAACATGTCCAGCCAGAAAAATAACGCCCCCCCCCCCAAAAACAGAGCCCAAGCAAAAGGCAAACCGGTAACaaaa                                                                    >  2:1983964/1‑236 (MQ=255)
                                                                                                                                                                           ccGCCTGGTTTATGCCCCACCGGTTTCATCGGCAGGTCAACAAGGCGTTTG‑G‑ACAGGCCGCCCGCATTCGTGTTTCCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATacc                                                                  <  1:2019864/239‑1 (MQ=255)
                                                                                                                                                                                   tttATGCCCCACCAGTTTCATCGGCAGGGCAACAATGCGTTTT‑G‑ACAAGCCGCCCGCATTCATGATTTCACCCGCCAGCACAAAGAACGGCATCGCCAGCAGGGAGAAGGTATCGGCCCCGCTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCCCCCCGCACAACAACAATGCACAGGCAATAAGCAAACCGAAAGCAAAACCACCCAaca                                                          >  2:2030685/1‑239 (MQ=255)
                                                                                                                                                                                         ccccACCATTTTCATCGCCAGCCAAACAATGCATTTG‑G‑AAAAGCCGCCCGCATTCATGATTTCCCCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGc                                                    <  1:2069592/239‑1 (MQ=255)
                                                                                                                                                                                            caccaGTTTCATCGGCAGGTCAACAATGCGTTTT‑G‑ACCAACCCCCCCCATTCCTTATTTCCCCCGCCAGCACAAAAAAAGGAAACGCCAGCAGGGGGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGGCAACAAACATTTCCAGCCCGAACATTAACGCCCCCCCGCACAACAACAATGCCCCGGAAAAAGGCCAAACGCAAACAAAAACACCCAAAAAACAAccc                                                  >  2:1225796/1‑238 (MQ=255)
                                                                                                                                                                                             accaGTTTCATCGGCAGGTCAACAATGCGTTTT‑G‑ACAAGCCGCCCGCCTTCATGATTTCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCCTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACCCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGccccg                                                >  1:975457/1‑237 (MQ=255)
                                                                                                                                                                                                        cggaagaacaacaaGGCGTTTG‑G‑AAAGGCCGCCCGATTTCATGATTCCCCCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGca                                     <  1:1551106/231‑1 (MQ=255)
                                                                                                                                                                                                        agaaaaaccaacaaTGCTTTTG‑C‑AAAACCCCCCCGCATGCATGCTTCCCCCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGGGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGca                                     <  2:1028504/231‑1 (MQ=255)
                                                                                                                                                                                                         ggCAGGTCAACAATGCGTTTT‑G‑ACAAACCGCCCGCCTTCATGATTTCACCCGCCAGCACAAAGAACGGAAACGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACAAGTCCAGCCAGAAAATTAACGCCGCCCCCCACAACAACAATGCCCCAGCAAAAGGCAAAACGAAAGCCAAAACACCCAACAAACAACCCAGAAAAAACAGcaa                                    >  2:883868/1‑238 (MQ=255)
                                                                                                                                                                                                          aacagcacaaaaTGCATTTTAG‑ACAAGCCGCCCGATTTCATGATTCCCCCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGc                                      <  1:883868/227‑1 (MQ=255)
                                                                                                                                                                                                           cAGGTCAACAATGCGTTTT‑G‑ACAAACCGCCCGCATTCATGATTTCACCCGCCAGCACAAAGAACGGAATCGCCAGCCGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCCCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGc                                                                        >  1:276667/1‑201 (MQ=255)
                                                                                                                                                                                                            gggTCAACAAGGCGTTTT‑A‑CCAAGCCGCCCGCATCCATGTTTTCCCCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGc                                 <  1:1622644/238‑1 (MQ=255)
                                                                                                                                                                                                               tCAAAAATGCGTTTG‑G‑ACCACCCCCCCCAATCCATGATTCCCCCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCTGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGCCa                               <  1:2231349/238‑1 (MQ=255)
                                                                                                                                                                                                               ccaaaacggCGTTTA‑G‑AAAAGCCCCCCGCATTCAGGATTCCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCGATCGGCTCCGTTCACCAGCGTTTGCGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGCCa                               <  1:1078201/230‑1 (MQ=255)
                                                                                                                                                                                                               cacacaagggttttt‑GAACAACCCCCCCGACTTCATATTTCCCCCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGCCa                               <  2:844289/228‑1 (MQ=255)
                                                                                                                                                                                                                   caaTGCGTTTT‑G‑ACAAGCCGCCCGCATTCATGATTTCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGCCATGAtt                          >  1:415249/1‑239 (MQ=255)
                                                                                                                                                                                                                        cGTTTT‑G‑ACAAGCCGCCCGCATTCATGATTTCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGCCATGATTAAGc                      >  2:103117/1‑238 (MQ=255)
                                                                                                                                                                                                                           ttt‑G‑ACAAGCCGCCCGCATTCATGATTTCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGCCACGAttctgactct                  >  1:1028504/1‑231 (MQ=255)
                                                                                                                                                                                                                             t‑G‑ACAAGCCGCCCGCATTCATGATTTCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGCCATGATTAAGCTCCtt                 >  2:2173694/1‑238 (MQ=255)
                                                                                                                                                                                                                                    aGCCGCCCGCATTCATGATTTCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGCCATGATTAAGCTCCTTGCGGTg           >  1:1655744/1‑239 (MQ=255)
                                                                                                                                                                                                                                     gccgccCGCATTCATGATTTCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGCCATGATTGAGCTCCTTGCGGTg           >  2:1630039/1‑238 (MQ=255)
                                                                                                                                                                                                                                     cccgcccgcATTCAATTTTCCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGCCATGATTAAGCTCCTTGCGGTgg          <  1:157438/238‑1 (MQ=255)
                                                                                                                                                                                                                                      ccgcccgcATTAATGATTCACCCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGCCATGATTAAGCTCCTTGCGGTgg          <  1:1983964/238‑1 (MQ=255)
                                                                                                                                                                                                                                             cATTCATGATTTCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGCCATGATTAAGCTCCTTGCGGTGGCGACGTTa  >  2:2092547/1‑239 (MQ=255)
                                                                                                                                                                                                                                             cATTCATGATTTCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGCCATGATTAAGCTCCTTGCGGTGGCGACGTTa  <  1:103117/239‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
ATTGGCGCGATAATGCCGCCAGAGGCAATCAGCCCCGCCGCCCGGTTGACCGGATAGTTGGCGCTGCGCATCATCGGCACCAGTAGAGCGGCAACGGCGGCGGTATCCGCAACGGCAGAACCAGAAAGGCTCGCCATAATCATAGCTGCCAGCACGCCGACGTAGCCCAGCCCGCCTGGTTTATGCCCCACCAGTTTCATCGGCAGGTCAACAATGCGTTTT‑G‑ACAAGCCGCCCGCATTCATGATTTCACCCGCCAGCACAAAGAACGGAATCGCCAGCAGGGAGAAGCTATCGGCTCCGTTCACCAGCGTTTGTGCCATGATCTGGACATCAAACATGTCCAGCCAGAACATTAACGCCGCCCCGCACAACAACAGTGCCCAGGCAATAGGCAAACCGATAGCAATACCACCCAACAGACAGCCCAGAAAAATCAGCACAGCCATGATTAAGCTCCTTGCGGTGGCGACGTTA  >  NC_007779/3895166‑3895639

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: ‑