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,4090TG97.9% 172.3 / ‑4.2 47intergenic (‑758/+28)sgcX/yjhPputative endoglucanase with Zn‑dependent exopeptidase domain/putative methyltransferase YjhP
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (31/15);  total (31/16)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

AAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCAT  >  NC_000913/4532168‑4532648
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
aaaCCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGccagcc                                                                                                                                                                                                                                           >  1:302714/1‑248 (MQ=255)
     tGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATAc                                                                                                                                                                                                                                   >  1:312920/1‑251 (MQ=255)
                       gggTTTGAGAGTGCCGTTGCTGTCGCTGCATTACGCACAATATGGATCACTGGCACTGTGTGCATCTGGTGCCGACATCTGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCGCCTTATGGTGACATCTTGATTAAGAAACATACTGAATATGTATTGTGCAGTGGATACATCAGATCAAAGCAAAGGGAAAGCCTCGATATGGCCAGCCAGAAATGCCGGGAACCTGTTGTTAcc                                                                                                                                                                                                                 <  2:168269/251‑1 (MQ=255)
                             gagaTTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCg                                                                                                                                                                                                            >  1:252880/1‑250 (MQ=255)
                                   gCCGTTGCTGTCGCTGCATTACGCACAATAGTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGt                                                                                                                                                                                                      <  2:66907/250‑1 (MQ=255)
                                     cGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAACTACCGGGAACCTGTTGTTACCGGGCGATGAGTGCg                                                                                                                                                                                                   >  2:86787/1‑251 (MQ=255)
                                           tGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCg                                                                                                                                                                                                   >  1:220206/1‑245 (MQ=255)
                                           tGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAAGAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGGAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCg                                                                                                                                                                                                   <  2:220206/245‑1 (MQ=255)
                                            gTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAacac                                                                                                                                                                                              <  2:283247/249‑1 (MQ=255)
                                            gTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGc                                                                                                                                                                                            >  2:163129/1‑251 (MQ=255)
                                               gCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGcccc                                                                                                                                                                                         >  2:252885/1‑251 (MQ=255)
                                                  gCATTACGCACAATATTGATCACTGGCACTGTGGGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTTTGCAGGTGATACATCAGATCAAAGCAAAGGTATAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAg                                                                                                                                                                                       <  1:30106/250‑1 (MQ=255)
                                                        cGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGc                                                                                                                                                                                >  1:330553/1‑251 (MQ=255)
                                                        cGCACAATATTGATCACTGGCACTGCGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAACCATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGc                                                                                                                                                                                >  2:169996/1‑251 (MQ=255)
                                                                atTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCcgc                                                                                                                                                                        >  1:196886/1‑251 (MQ=255)
                                                                   gATCACTGGCCCTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGTCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCcgcgcg                                                                                                                                                                     <  1:246040/251‑1 (MQ=255)
                                                                      cACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGt                                                                                                                                                                    >  1:319565/1‑249 (MQ=255)
                                                                      cACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAg                                                                                                                                                                  <  2:274974/251‑1 (MQ=255)
                                                                           gCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACAGCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCAc                                                                                                                                                             <  2:330553/251‑1 (MQ=255)
                                                                                gtgtgCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGc                                                                                                                                                        <  1:210559/251‑1 (MQ=255)
                                                                                        cTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATATTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCCAAGCATTCCCGCGCGTAGGTCCCGTAGCGTTTCGgc                                                                                                                                                >  1:192626/1‑251 (MQ=255)
                                                                                          ggTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGgcgc                                                                                                                                              >  1:319687/1‑251 (MQ=255)
                                                                                          ggTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGgcgc                                                                                                                                              >  2:251416/1‑251 (MQ=255)
                                                                                           gTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATATTGATTAATAACCATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTCCCGGGCGATAAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGcgcg                                                                                                                                             >  2:131076/1‑251 (MQ=255)
                                                                                                 aCATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGACAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCACCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGttt                                                                                                                                       >  2:327330/1‑251 (MQ=255)
                                                                                                              gaaaGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATAATGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGACAGAAATATCGGGAACCTGTTGTTACCGGGAGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGGTATGTTTAGTTCGGCCCTGa                                                                                                                          <  2:236543/249‑1 (MQ=255)
                                                                                                               cAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGAc                                                                                                                         <  1:13976/251‑1 (MQ=255)
                                                                                                               cAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTACTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGAc                                                                                                                         >  2:19542/1‑251 (MQ=255)
                                                                                                                   ccTGACAGAACTGTTCATGTTGAGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGCTGAGTCCGACCACGCCCCAGCCACAGCATTCCCGCGCGGAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTcg                                                                                                                     >  2:240141/1‑251 (MQ=255)
                                                                                                                            aCTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGa                                                                                                              >  1:282042/1‑249 (MQ=255)
                                                                                                                              tGTTCAAGTTGTGACAAATGATGCTCATTCTGGTGACATCTTGATTAATAAACATAATGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCAAGCCCGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAACGCATTACCGAGCGTAGGTCACGTAGCGTTAAGGAGCGATGTTAAGATCGGACCTGACTTCGGCGGCGAAGTc                                                                                                          >  2:98584/1‑251 (MQ=255)
                                                                                                                                               atgatgCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCtt                                                                                         <  1:41533/251‑1 (MQ=255)
                                                                                                                                                                 aCATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATgg                                                                       >  1:291530/1‑251 (MQ=255)
                                                                                                                                                                  cATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGt                                                                      >  1:332801/1‑251 (MQ=255)
                                                                                                                                                                      ttGATTAATAAACATACTGTATATGTCTTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGCACACGCACCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCCTGGTCAGc                                                                  >  2:302122/1‑251 (MQ=255)
                                                                                                                                                                           taataaACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCAttt                                                             >  1:23358/1‑251 (MQ=255)
                                                                                                                                                                                        gAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATAc                                                >  1:223731/1‑251 (MQ=255)
                                                                                                                                                                                         aaTATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATAcc                                               <  2:312920/251‑1 (MQ=255)
                                                                                                                                                                                            atGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCTGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGt                                            <  2:332801/251‑1 (MQ=255)
                                                                                                                                                                                            atGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGACGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCCTTTTGCGGCTTCCTACCGGt                                            >  2:298189/1‑251 (MQ=255)
                                                                                                                                                                                            atGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGt                                            >  1:269330/1‑251 (MQ=255)
                                                                                                                                                                                              gTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTcc                                          >  1:313557/1‑251 (MQ=255)
                                                                                                                                                                                                   gtgCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGcc                                     >  2:64449/1‑251 (MQ=255)
                                                                                                                                                                                                              aCATAAGATCAAAGCAAAGGTAAAGCCTCGATATTGCCAGCCAGAAAGACCGGGAACCTTTTGTTACCGGGCGATGAGTGCGAACCCGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCt                           <  2:113946/250‑1 (MQ=255)
                                                                                                                                                                                                                        aaagcaaagGTAAAGCCTCGATATTGCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCAt                <  2:282042/251‑1 (MQ=255)
                                                                                                                                                                                                                              aagGTAACGCCTCGATATTGCCAGCCAGACATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGACCACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCCTGGTCAGCCATTTTGCGGCTTCCTACCGGTCCCAGCCTTCCTGGTCTGCCCGTACCATTTccac          >  2:19443/1‑251 (MQ=255)
                                                                                                                                                                                                                                      gCCTCGCTATTGCCCGCCAGAAATACCGGGAACCTGTTGTTACCGGACGATGAGTGCGCACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGCCAGCCATTTTGCGGCTTCATCCCGGTCCCCGCCTTCCTGGTCTGCCAGTACCCTTTCCACCACGTCAt  >  2:204283/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
AAACCTGATACCGCTTCTTCCCAGGGTTTGAGATTGCCGTTGCTGTCGCTGCATTACGCACAATATTGATCACTGGCACTGTGTGCATCTGGTGCCGACATCGGCATTCCGCAAGCCTGACAGAACTGTTCATGTTGTGACATATGATGCTCCTTATGGTGACATCTTGATTAATAAACATACTGAATATGTATTGTGCAGGTGATACATCAGATCAAAGCAAAGGTAAAGCCTCGATATTTCCAGCCAGAAATACCGGGAACCTGTTGTTACCGGGCGATGAGTGCGAACACGCCCCAGCCAAAGCATTCCCGCGCGTAGGTCACGTAGCGTTTCGGCGCGATGTTTAGTTCGGCCCTGACTTCGGCGGCGAAGTCATCGTCAGGATTCGCTTCCAGCCAGCGGCGCATGGTCAGCCATTTTGCGGCTTCATACCGGTCCCAGCCTTCCTGGTCTGCCAGTACCATTTCCACCACGTCAT  >  NC_000913/4532168‑4532648

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