Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,562,499 T→A 15.9% I49F (ATC→TTC)  Y75_RS07840 ← ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077791,562,4990TA15.9% 133.0 / 2.3 58I49F (ATC→TTC) Y75_RS07840ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (41/7);  new base A (7/3);  total (48/10)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAACACAAAGCCTTTCACCTCTTTGTTCATCGCCAGTTGGTAG  >  NC_007779/1562268‑1562735
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
gCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCg                                                                                                                                                                                                                                         <  1:871627/239‑1 (MQ=255)
                tATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTTCACGATAGCATCGCTGGCGCGGTCAccc                                                                                                                                                                                                                          >  1:594476/1‑238 (MQ=255)
                tATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTTGCGCGGGCAccc                                                                                                                                                                                                                          >  1:1650275/1‑238 (MQ=255)
                   gCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGGCCAGGCCCAGTTGCTGGCGAATATTTTCCACGATAACATCCCTGGCGCGGTCAccccccc                                                                                                                                                                                                                      >  1:1091565/1‑238 (MQ=255)
                       aGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGGCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGGCACCCGCCCGTa                                                                                                                                                                                                                   >  1:626514/1‑238 (MQ=255)
                        gCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAACTATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGGCCCCGGTGCTGGGCAATATTTTTCACCATAACATCGCCCGCGCGGTCATCCGCCCGTaa                                                                                                                                                                                                                  >  1:502395/1‑238 (MQ=255)
                        gCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTTCACGATAGCATCGCTGGCGCGGTCCCCCGCCCGTaa                                                                                                                                                                                                                  >  1:1454609/1‑238 (MQ=255)
                              tGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTGCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGcc                                                                                                                                                                                                           >  1:347383/1‑239 (MQ=255)
                                   aTCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGGCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTTGCCCGGGCACCCCCCCGGAACCGCGCCGGATc                                                                                                                                                                                                      >  1:1651666/1‑239 (MQ=255)
                                      aggggggCCAGCGGCAGAGCACCAAAAGCCAGCCCCAGCGCCGCCAGAAAAAATACCCCCAACTCTCCCAGCCCGGTCCCCCCGGTACGAGTGGAGGTCCCCAGGACACCATGAAGCAGGTCGCTGACGTGGAGTAAAAACTGTACACACGGGGGCTGGCCCCCCCCCGGTGGCGGGAGAATTTTTCCCCGGAAAGCCGCGCTGGCGCGGCCACCCGCCAGTAACCGCGCCGGATCAc                                                                                                                                                                                                    <  2:801729/232‑1 (MQ=255)
                                             ccAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATc                                                                                                                                                                                            >  1:148269/1‑239 (MQ=255)
                                               agcagcagGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTGCCGTCCCCGTTTCTGGCGAAAATTTTTCACGAAAGCAACGCTGGCGCGGTCCCCCCCCAGAACCCCCGCCCGCACACCAGGGGTCac                                                                                                                                                                                          >  2:640662/1‑238 (MQ=255)
                                                                    aGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTaa                                                                                                                                                                      >  1:368772/1‑238 (MQ=255)
                                                                                gCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGAcgc                                                                                                                                                          >  1:713919/1‑238 (MQ=255)
                                                                                                            gCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAAAATTTTTCACCATAGCATCGCTGGCGCGGGCCCCCCCCCGTAAACGCGCCGGGTCACCAGGGATCAGATGCGGGATAATAAAAGAAAACACGCAGACCCCCGCCACCCCCAGCACCCCCCCCcaaca                                                                                                                             >  2:722336/1‑238 (MQ=255)
                                                                                                               ccGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGt                                                                                                                         >  1:824308/1‑239 (MQ=255)
                                                                                                                               aaTGGATGTTCCCAGGTCCCCATGAAACAGGCCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTCCCACGATAGCATCGCTGGCGCGGt                                                                                                                                                                                                                               <  1:310171/122‑1 (MQ=255)
                                                                                                                                 tGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAg                                                                                                      >  1:1481383/1‑239 (MQ=255)
                                                                                                                                   gATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGGCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGaa                                                                                                    >  1:380378/1‑239 (MQ=255)
                                                                                                                                     tGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCCGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCCGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAgg                                                                                                  >  1:1471573/1‑239 (MQ=255)
                                                                                                                                     tGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCCGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAgg                                                                                                  >  1:666451/1‑239 (MQ=255)
                                                                                                                                     tGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAgg                                                                                                  >  1:1857003/1‑239 (MQ=255)
                                                                                                                                     tGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGGCCCGGCCCCGTTTGTGGCGGATTCTTTCCACCAAAACATCGCTCGCGCGGGCCCCCGCCCGCAACCGCGCCCGGTCCCCCGGGATCCCAAGCGGGATAAGAAACCAAAACACCCAGGAGCCCCCCCCCCCCACCCCCCACCCCCAGCAACGTTCGCGGTAAAAACCTcccaaaa                                                                                                   >  2:1615895/1‑234 (MQ=255)
                                                                                                                                     tGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGGCCAGTCCCAGTTGCTGGCGGATATTTTCCACCAAAGCCTCGCCGGGGCGGTCACCCGCCCGTAAACGCGCCCGCGCACCCGGGGACCGATGGGAGATAATAAACGAAAACACGCCGACGCACCCCACCCCCAGCACCCACCCCCCAAGCCGTTGGCGGAAAAAACACCAAAAgg                                                                                                  >  2:1807063/1‑239 (MQ=255)
                                                                                                                                       tcccaaGGTCACCAGGAAACGGGTCCCTGACGTAGCGGAAAACCGTTCCATACAGTGGCGGGCCCAGCCCCGGTTGCGGGAAATTTTTTCCCAGAATAGCACCGCGGGCGCGGCCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTc                                                                                                <  1:1806756/234‑1 (MQ=255)
                                                                                                                                         cccAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCCGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCAc                                                                                              >  1:76806/1‑239 (MQ=255)
                                                                                                                                         cccAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCCGTTGCTGGCGAATATTTTCCCCGAAAGCATCGCTGGCGCGGTCCCCCGcccg                                                                                                                                                                                                                     >  2:310171/1‑120 (MQ=255)
                                                                                                                                              gTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCGCGTTGGCGTAAAATACTCCAGAAGGGCACGTTAtt                                                                                        >  1:1626656/1‑239 (MQ=255)
                                                                                                                                                  ccATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATAACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGGTGGCGTAAAATACTCCAGAAGGTCACGTTAttt                                                                                       >  1:336309/1‑237 (MQ=255)
                                                                                                                                                   cATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGGCACGTTATTTaa                                                                                     >  1:1184417/1‑237 (MQ=255)
                                                                                                                                                           gggTCGCTGCGGTGGGGTTAAACCTTTACATACGGGGGCGGGTCCAGCCCCGGTTGCTGGAAAATTTTTCCCACGATAGCACCGCGGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtat                                                                            <  1:722336/238‑1 (MQ=255)
                                                                                                                                                             gTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGGGAAAATTTTCCCCGATAGGATCGCTGGCGCGGGCACCCGCCAGGAAACCCGCCGGAGCACCCGGGGTCAGATGGGGGATAATAAACGGAAAACCCCCGACGCCCGCCACCACCAGCACCCAACCCCCAGAGCGTTGGCGTAAAATACCCCCGAAGGGCACGGTATTTTACCATGGttttg                                                                          >  2:412987/1‑239 (MQ=255)
                                                                                                                                                             gTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGGCCAGCCCCAGGTGCTGGGGAATATTTTCCACGATAGACTCGCTCGCGCGGTGACCCCCCAGGAACCGCGCCCGATCACCCCGGATCAGACTGGAGATAATAAACGTAAACACGCAGACGCCCGCCACCACCAACCACCGCCCCCCACCGCGGTGGCGGAAAAAACCCCCCAAAGTCACGTTATTTACTCATGGtaatg                                                                          >  2:983363/1‑239 (MQ=255)
                                                                                                                                                                       tAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGAAAGCATCGCTGGCGCGGTCACCCGCCAGGAACCGCGCCGGATCACCAGGGATCAGATTCGGGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGCAGGTCACGTTATTTACCCTTGTGTTTAATATTGaaa                                                                 >  1:588619/1‑237 (MQ=255)
                                                                                                                                                                         gCGGTAAAACTGTACATACAGTGGCTGGGCCAGTCCCAGTTGCCGGCGAATATTTTCCACGATAACATCGCTGGCGCGGTCACCCCCCAGCAAACGCGCCCGGACACCCAGGGTCCGCTGCGAGATAATAAACGGAAACACCCCGACGCCCCCCCCCACGAGGCCCAACCCCCCGCAGCGGTGGGGCAAAAAACCCCCAAAAGCCACGTTATTTTCTCATCGGATTGATATTTAAAAc                                                               >  2:1222304/1‑238 (MQ=255)
                                                                                                                                                                         gCGGAAAACGTGAACAAACTGTGGCCGCCCCCGCCCCAGTGGCGGGACAATTTTCCCCGCGAAGGCAGCGGTGGGGCGCCCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGAGGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGAc                                                               <  2:76806/238‑1 (MQ=255)
                                                                                                                                                                                   tGTAAATACAGTGGCGGGCCCTCCCCCGGTTGCGGGAAAAATTTCCCCGCGAAAGCCGCGGGGGCGCGCACACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCa                                                    <  1:1235538/239‑1 (MQ=255)
                                                                                                                                                                                     tacatacaGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGGCACCCGCCAGTAACCGCGCCGGATCACCAGGGAt                                                                                                                                                                                             >  1:415156/1‑102 (MQ=255)
                                                                                                                                                                                      acatacaGTGGCTGGTCCAGCCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCGGGCGCGGTCACCCGCCAGTAACCGCGCCgg                                                                                                                                                                                                         <  1:1370800/89‑1 (MQ=255)
                                                                                                                                                                                            aGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACCAAAGCATCGCTGGCGCGGTCACCCCCCCGTAACCGCGCCCGCACACCAGGGATCAGATGCGGGATAATAAACGTAATCACGCAGACGCCCGCCACCACCAGCACCAGCCCCCCCCCGCGTTGGCGTAAAAAACCCCAGAAGGGCACGGTATTTTCTCATGGTATTGATATTTAAGACCTGTTTCAGCATGGGATTGa                                           >  2:750527/1‑239 (MQ=255)
                                                                                                                                                                                               ggCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAac                                         >  1:941595/1‑238 (MQ=255)
                                                                                                                                                                                               ggCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGGCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGGCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAac                                         >  1:287038/1‑238 (MQ=255)
                                                                                                                                                                                               ggCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACAAAAGCATCGCTGGCGCGGGCACCCGCCAGTAAACGCGCCGGGAAACCAGGGAACAGAAGCGGGAGAATAAAAGTAATCCAGCAGACCCCCGCCACCACAAGCACCAGCCCACAGCAGCGTTGGCGGAAAATACTCCCGAAGGGCAAGTTATTTACTCCTGGTATTGATATTGAAGAACCGTTCCAGCATGGGATTGAac                                         >  2:1610753/1‑238 (MQ=255)
                                                                                                                                                                                               ggCTGGGCCCGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCCCCCCCCCGTAACCGCGCCGGAACACCAGGGATCAGATGCGAGATAATAAACGGAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAAAACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCACCCTGGGTGTGAac                                         >  2:1737390/1‑238 (MQ=255)
                                                                                                                                                                                                  tGGTCCAGCCCCGTTTGCTGGAAAATATTTCCCACGATAGATTGGCGGGCGCGGCCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAACACaa                                     <  1:479724/239‑1 (MQ=255)
                                                                                                                                                                                                   ggTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCCGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAACACcaa                                    >  1:853993/1‑239 (MQ=255)
                                                                                                                                                                                                    gTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCCGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGtt                                                                                           >  1:1771717/1‑183 (MQ=255)
                                                                                                                                                                                                          tCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAACACAAAGCCTTTc                             >  2:18473/1‑239 (MQ=255)
                                                                                                                                                                                                          tCCCACTTGCTGGCGAATATTTTCCACGCTAGCATCGATTGCCCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATAAGATGCGCGATAATAATCGTAATCACGCAGACGCCCGCCACCACGAGCACCCGCCCCCAGCAGCGTTGGCGCAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTCTTGATATTAACGACCTGCTCCAGCATGGGATTGATCaaaacatctttg                             >  2:1000190/1‑231 (MQ=255)
                                                                                                                                                                                                              aGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattgat                                                                        >  1:826996/1‑192 (MQ=255)
                                                                                                                                                                                                               gTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGcc                                                                                                                                    >  1:72113/1‑132 (MQ=255)
                                                                                                                                                                                                               gTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAACACAAAGCCTTTCAcc                          >  2:85328/1‑237 (MQ=255)
                                                                                                                                                                                                                    tGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAACACAAAGCCTTTCACCTCTTTGt                   >  1:1414407/1‑239 (MQ=255)
                                                                                                                                                                                                                     ggCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCAc                                                                                                                                                                                                    >  2:453003/1‑63 (MQ=255)
                                                                                                                                                                                                                     ggCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCAc                                                                                                                                                                                                    <  1:453003/63‑1 (MQ=255)
                                                                                                                                                                                                                     ggCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAACACAAAGCCTTTCACCTCTTTGt                   >  1:175243/1‑238 (MQ=255)
                                                                                                                                                                                                                                tCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAACACAAAGCCTTTCACCTCTTTGTTCATCGCCAGtt       >  1:694360/1‑239 (MQ=255)
                                                                                                                                                                                                                                     gATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAACACAAAGCCTTTCACCTCTTTGTTCATCGCCAGTTGGTAg  >  1:227636/1‑239 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
GCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAACAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGTCCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGTAATCACGCAGACGCCCGCCACCACGAGCACCAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAACACAAAGCCTTTCACCTCTTTGTTCATCGCCAGTTGGTAG  >  NC_007779/1562268‑1562735

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: