Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,532,409 T→G intergenic (‑758/+28) sgcX ← / ← yjhP putative endoglucanase with Zn‑dependent exopeptidase domain/putative methyltransferase YjhP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,532,4090TG94.2% 161.9 / ‑0.9 52intergenic (‑758/+28)sgcX/yjhPputative endoglucanase with Zn‑dependent exopeptidase domain/putative methyltransferase YjhP
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (26/23);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.44e-01

TAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCC  >  NC_000913/4532167‑4532652
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
tAAACCTGATACCGCTTCTTCCCAGGGTTTGCGATTGCCGTTGATGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCCGATCAAAGCAAAGGTAAAGCCTCGATATTGccagcc                                                                                                                                                                                                                                               >  2:290376/1‑249 (MQ=255)
   aCCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGaaa                                                                                                                                                                                                                                          <  1:216860/251‑1 (MQ=255)
            cgattcttcCATGGGTTTGATAGTGTCGTTGCTGGCGCTGCATTACGCACAATAGTTAGCACTGGCACAGTGTGCATCTTTTGCCGACATCTGCATTCCGCAAGTCTGACAGAACTGGTCATGTTGTGTCATATGATGCTCCTGATGGTGACATCTTTATTAATAAACATACTGAATATGTATTGTTAAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCTATATTGCAAGCCAGAAATACCGGGaa                                                                                                                                                                                                                                 <  2:47099/248‑1 (MQ=255)
            cGCTTCTTTCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGATATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAAGAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCATCCAGAAATACCGGGaa                                                                                                                                                                                                                                 <  1:262600/251‑1 (MQ=255)
                  ttcCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCtgtt                                                                                                                                                                                                                           >  1:191607/1‑251 (MQ=255)
                                   tGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGGGATGAGt                                                                                                                                                                                                          >  1:120583/1‑251 (MQ=255)
                                    gCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCGCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTg                                                                                                                                                                                                         <  2:159037/251‑1 (MQ=255)
                                     ccGTTGCAGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGc                                                                                                                                                                                                        <  1:7696/251‑1 (MQ=255)
                                          gCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCAGTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAaca                                                                                                                                                                                                   <  1:215493/251‑1 (MQ=255)
                                                             caATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGACCTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGAAAAGCCTCGATATTGCCCGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCTTTc                                                                                                                                                                                >  2:232730/1‑251 (MQ=255)
                                                                      tCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGt                                                                                                                                                                        >  2:291192/1‑250 (MQ=255)
                                                                         cTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGt                                                                                                                                                                    >  1:147838/1‑251 (MQ=255)
                                                                            gCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGACCCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCAc                                                                                                                                                                 >  2:279007/1‑251 (MQ=255)
                                                                             cACTGTGTGCATCTGGTGCCGACATCTGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGCCATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGAGAAGGTCACg                                                                                                                                                                >  2:272145/1‑251 (MQ=255)
                                                                                 gtgtgCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGc                                                                                                                                                            <  1:12676/251‑1 (MQ=255)
                                                                                   gtgCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCAGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGt                                                                                                                                                          >  1:22193/1‑251 (MQ=255)
                                                                                    tgCATCTGGTGCAGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATAAATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACACCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGAAGCgct                                                                                                                                                         >  2:96391/1‑249 (MQ=255)
                                                                                                     acGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGt                                                                                                                                        <  1:180074/250‑1 (MQ=255)
                                                                                                              cGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTg                                                                                                                               <  1:76458/251‑1 (MQ=255)
                                                                                                                cAAGCCTGACAGAACTGTTCATGTTGTGAAATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGAGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCAAGCCAAAGCATTCCCGCGCGTAGGTCCCGTAGCGTTTCGGCGTGTTGTTTAGTTAGGCCCTGAc                                                                                                                             >  2:90916/1‑251 (MQ=255)
                                                                                                                 aagccctaCAGAACTGTTAATGTTGTGACATATGATGTTCCTTATGGTGACATCTTGGTTAATAAATATACTGAATATGTATTGTGCAGGGGATACAGCAGATCAAAGCAAAGGTAAAGCCGAGAAATTGCCAGCCAGAATTACCGGGAACCTGTTGTTCCCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCAGTCCCGCGCGTAGGTCATGTAGCGTTTCGTCGCGATGTTTAGTTCTGCCCTgcct                                                                                                                            <  2:200586/244‑4 (MQ=255)
                                                                                                                            aaCTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAGTATGTATTGTGCAGGTAATACATCAGATAAAAGCAAAGGTAACGCCTCGATATTGCAAGCCAGAAATACCGGGAACCTGATGTTACCGATCGATGAGCGCGAACACGCCCCCGCCACAGCATTCCCCCGCGTAGGACACGTCGAGTTTCGGCGCTATGTTTAGTTCGGCCCtcaatgaggcagcgct                                                                                                                 >  2:165330/1‑236 (MQ=255)
                                                                                                                                  tCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAACGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATc                                                                                                           >  2:126958/1‑251 (MQ=255)
                                                                                                                                   cATGTTGTGACATATGATGCTCCTTATGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGt                                                                                                         >  1:144396/1‑251 (MQ=255)
                                                                                                                                    agcttgtGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTTTGCAGGTGATCCAGCAGATAAAAGCAAAGGGAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGGTCGGCCCTGACGTCGGCGGCGAAGTCATCGt                                                                                                         <  2:154637/248‑1 (MQ=255)
                                                                                                                                            acaaatgattcctcTTATGGTGACATATTGATTAATAAACATACTGAATATGTAGTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGATCCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCTCGTAGGTCACGTAGCGTTTCGGCGCTATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTc                                                                                                 <  2:152768/238‑1 (MQ=255)
                                                                                                                                            aCATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTc                                                                                                 <  1:296080/251‑1 (MQ=255)
                                                                                                                                                     gCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTccag                                                                                         <  1:291192/250‑1 (MQ=255)
                                                                                                                                                                          aTTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAt                                                                   >  1:111907/1‑251 (MQ=255)
                                                                                                                                                                          aTTAATAAACATACTGAATATGTATTGTGCAGGTAATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTCCCAGCCAGCGGCGCATGGTCAGCCAt                                                                   <  1:231370/251‑1 (MQ=255)
                                                                                                                                                                            taataaACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAtt                                                                  <  1:102754/250‑1 (MQ=255)
                                                                                                                                                                                 aaCATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCAAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTCGCGTTTCGGCGCGCTGTTTAGTTCGGCCCTGACTTCGGCGGGGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCCTGGTCAGCCATTTTGCgg                                                            >  2:66922/1‑251 (MQ=255)
                                                                                                                                                                                   cATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGCTGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCt                                                          >  1:291582/1‑251 (MQ=255)
                                                                                                                                                                                   cATACTGAATATGTAGTTTGAAGGTGATACTACAGATCAAAGCAAATTTAAAGTCACGATATTGCCAGCCAGAAATACCGGGAACCTGTTGGTACCGGGCTAGGAGTGCGAACACGCCCCAGCCAAAGCATGCCCGCGCGTAGGTCACGTAGCGGTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTTCGGCt                                                          <  2:287922/251‑1 (MQ=255)
                                                                                                                                                                                    aTACTGAATATGTATTGGGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCtt                                                         <  1:161250/251‑1 (MQ=255)
                                                                                                                                                                                       cTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAATGCGAAAACGCCCCAGCCAAAGCATTCCCGCGAGTAGGTCACGTAGCGTTTCGGCGCGATGTTTCGTTCGGCCCTGACTTCGGCGGCCCAGACATCGACAGGATTCCGTTCCACCCCTCGACGCATGACCAGCCATTTTGCGGcatct                                                       >  2:36961/1‑246 (MQ=255)
                                                                                                                                                                                         gAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATAc                                                    >  1:263183/1‑251 (MQ=255)
                                                                                                                                                                                                    gtgCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGGAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGGCAGCCATTTTGCGGCTTCATACCGGTCCCCGcc                                         >  2:234964/1‑251 (MQ=255)
                                                                                                                                                                                                      gCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACTTAGTGGTTCGGCGCGATGTTTAGTTCGGCCCTGTCTTCGGCGGCGAAGTCAGCGTCAGGATTCGCTTCCAGTCAGCGGCGCATGGTCAGcc                                                                     <  2:21390/221‑1 (MQ=255)
                                                                                                                                                                                                      gCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGcc                                                                     >  1:21390/1‑221 (MQ=255)
                                                                                                                                                                                                      gCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCtt                                       <  1:290376/251‑1 (MQ=255)
                                                                                                                                                                                                       cAGGTGATACATCAGATCAAAGCAAAGGGAAAGCCTCGATATTGCCAGCCAGAACTACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCAACGTCAGGGTTCGCTTCCCGCCAGCGGCGCCTGGGCCGCCGTTTTGCGGCATCCTACCGGGCCCAGCCTTc                                      >  2:55451/1‑251 (MQ=255)
                                                                                                                                                                                                         ggTGAAACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGACCCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCt                                    <  2:144396/251‑1 (MQ=255)
                                                                                                                                                                                                           tGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCCTACCGGTCCCCGCCTCCCTgg                                  >  2:118937/1‑251 (MQ=255)
                                                                                                                                                                                                            gATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCAAGCCAAAGCATTCCCGAGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGt                                 >  1:129088/1‑251 (MQ=255)
                                                                                                                                                                                                                  tCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGcc                            >  1:122892/1‑250 (MQ=255)
                                                                                                                                                                                                                        caaagcaaagGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGCACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTACGTCCCGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGTCGGCGAAGTCATCGTCCGGATTCGCTTCCAGCCAGCGGCGCAAGGTAAGCCTTATTGCTGCTTCATACCGGTCCCAGCATTCCTGATCTGCCAGTAcc                      >  2:226425/1‑250 (MQ=255)
                                                                                                                                                                                                                          aagcaaagGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCAtt                   <  1:126958/251‑1 (MQ=255)
                                                                                                                                                                                                                             caaagGTAAAGCCTCGATATTGCCAGCCAGCAATACCGGGAACATGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGCAGTCATCGTCAGGATTCGCTTCCAGCCAGAGGCGCATGGTCAGCCATTTTGCCGCTTCATACCGGTCCCAGCCTTCCTTGTCTGCCAGTACCATTTcc                >  2:38091/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTTCCATTTCcaccac           <  1:55451/251‑1 (MQ=255)
                                                                                                                                                                                                                                         cTCGATCTTGCCAGCCAGAAATACCGGGACCCTGTTGTTACCGGGCGATGAGGGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGCCTTCGGCGGCGAAGTCATCGTCAGGGTTCGCTTCCAGCCAGCGCCGCATGGTCAGCCCTTTTGCGGCTTCATACCGGTCCCACCCTTCCTGGTCTGCCGGTACCATTTCCACCACGTCAt      >  2:138246/1‑249 (MQ=255)
                                                                                                                                                                                                                                           cGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTTCCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGcc  <  1:53641/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
TAAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCATAGCC  >  NC_000913/4532167‑4532652

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