Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,355,550 G→A 9.0% A51V (GCG→GTG)  truB ← tRNA pseudouridine(55) synthase TruB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,355,5500GA 9.0% 182.5 / 5.7 66A51V (GCG→GTG) truBtRNA pseudouridine(55) synthase TruB
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base A (3/3);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

GCGGGTCTGCAGCACCTCACCCTCGGCATCGCCGGTGTTGGTAGTCTGCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCTGACATTGCGGCGGATACGTTTGACCTG  >  NC_002947/5355405‑5355695
                                                                                                                                                 |                                                                                                                                                 
gcgGGTCTGCAGCACCTCACCCTCGGCATCGCCGGTGTTGGTAGTCTGCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCa                                                                                                                                                <  7:36939/149‑1 (MQ=255)
         cagcaCCTCACCCTCGGCATCGCCGGTGTTGGTAGTCTGCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGa                                                                                                                                       <  2:129422/149‑1 (MQ=255)
               cTCACCCTCGGCATCGCCGGTGTTGGTAGTCTGCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTAc                                                                                                                                 >  1:127390/1‑149 (MQ=255)
                  aCCCTCGGCATCGCCGGTGTTGGTAGTCTGCCCCCTCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTACGAAAACTTGGTCGCTTCACCGAAGCACAGAGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCgg                                                                                                                              >  7:160968/1‑149 (MQ=255)
                   cccTCGGCATCGCCGGTGTTGGTAGTCTGCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCgg                                                                                                                              >  4:95303/1‑148 (MQ=255)
                      tCGGCATCGCCGGTGTTGGTAGTCTGCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATg                                                                                                                          >  3:14459/1‑149 (MQ=255)
                          cATCGCCGGTGTTGGTAGTCTGCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGACGGTCGCCAGCGGGGCGAGGCTACCGGTATgcccg                                                                                                                      >  6:223592/1‑149 (MQ=255)
                              gCCGGTGTTGGTAGTCTGCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAATCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACAGGTATTGCCGGCCt                                                                                                                  >  8:164768/1‑149 (MQ=255)
                              gCCGGTGTTGGTAGTCTGCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATggccggcc                                                                                                                   <  3:95303/148‑1 (MQ=255)
                                     ttGGTAGTCTGCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCg                                                                                                            <  1:152779/148‑1 (MQ=255)
                                        gTAGTCTGCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCAtt                                                                                                        >  6:182731/1‑149 (MQ=255)
                                           gTCTGCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTcag                                                                                                     >  1:197274/1‑149 (MQ=255)
                                           gTCTGCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTcag                                                                                                     >  8:70410/1‑149 (MQ=255)
                                              tGCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTcagcag                                                                                                  >  7:200702/1‑149 (MQ=255)
                                               gCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGACGGCCTTTTCCGCATTCagcagc                                                                                                 >  7:83362/1‑149 (MQ=255)
                                                cctgaTCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCGATCAAGTACTACGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCagcagc                                                                                                 <  8:160968/144‑1 (MQ=255)
                                                ccccATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGTCTTTTCCGCATTcagca                                                                                                   <  7:182884/146‑1 (MQ=255)
                                                ccccATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTcagcag                                                                                                  <  3:129349/147‑1 (MQ=255)
                                                ccccATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTcagca                                                                                                   <  4:73497/146‑1 (MQ=255)
                                                   cATCTGCATGACGGTTTCATAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCa                                                                                               >  1:27035/1‑147 (MQ=255)
                                                             aCGGTTTCATAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCt                                                                                    <  2:27035/148‑1 (MQ=255)
                                                               ggTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGc                                                                                  >  4:133502/1‑148 (MQ=255)
                                                                 tttCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAg                                                                               >  8:82280/1‑149 (MQ=255)
                                                                   tCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAgg                                                                              <  5:182731/148‑1 (MQ=255)
                                                                     gTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAggcgg                                                                           >  7:219295/1‑149 (MQ=255)
                                                                     gTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCACCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAggcgg                                                                           >  8:256607/1‑149 (MQ=255)
                                                                      tAGCCCTTGACGGAATCAAGCAAGTACTGCGAAAACTTGGTCGCTTCACCTAAGAACAGCGGAAAAACTCCGGTCGCCAGCGCGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCAAACGCACTTTATGCAAggcgt                                                                           >  4:230164/1‑147 (MQ=255)
                                                                           cTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTgg                                                                     <  7:82280/149‑1 (MQ=255)
                                                                           cTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTg                                                                      >  5:72331/1‑148 (MQ=255)
                                                                               tCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGt                                                                 <  3:133502/149‑1 (MQ=255)
                                                                                cGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCACCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTg                                                                                   >  4:165661/1‑130 (MQ=255)
                                                                                cGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCACCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTg                                                                                   >  4:160389/1‑130 (MQ=255)
                                                                                cGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCACCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTg                                                                                   <  3:165661/130‑1 (MQ=255)
                                                                                cGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCACCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTg                                                                                   <  3:160389/130‑1 (MQ=255)
                                                                                  gAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTg                                                                >  6:42869/1‑147 (MQ=255)
                                                                                   aaTCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAg                                                             <  7:70410/149‑1 (MQ=255)
                                                                                   aaTCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAg                                                             <  2:197274/149‑1 (MQ=255)
                                                                                    ctCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGc                                                            <  5:223592/148‑1 (MQ=255)
                                                                                      ccAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGACGcc                                                           >  8:305923/1‑148 (MQ=255)
                                                                                               agtgagAAAACTTGGTCGCTTCACCGAAGCACAGCGTCAACACGCCGGTCGCCAGCGCGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCtt                                                  <  3:230164/143‑1 (MQ=255)
                                                                                                cTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCtttt                                                >  5:303344/1‑147 (MQ=255)
                                                                                                cTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGt                                                <  1:337469/149‑1 (MQ=255)
                                                                                                   cGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTc                                               >  3:81656/1‑147 (MQ=255)
                                                                                                   cGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGa                                             <  4:14459/149‑1 (MQ=255)
                                                                                                      aaaCTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCACCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCa                                          <  7:256607/149‑1 (MQ=255)
                                                                                                                  ttCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGGATTCAGCAGCAAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCGAGAGGCTTGTCGAGCAGGATGATGgcgc                              >  4:81308/1‑149 (MQ=255)
                                                                                                                  ttCACCGAAGCACAGCGGCAACACGCCGCTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTTGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGc                              >  1:176591/1‑149 (MQ=255)
                                                                                                                    cACCTAAGCACAGCGGCAACACGCCGGTAGCCAGGGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCATCCAACGCACTTTCTGCAAGGCGGGATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCt                             <  8:83362/148‑1 (MQ=255)
                                                                                                                       cGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCTGACa                         >  4:164389/1‑149 (MQ=255)
                                                                                                                       cGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCAGGAAGCAGCCAACGCCCTTTATTCAAGGAGGCATATGTGGTGATGACCCCCCGCTGGTTGTGGAGGATGATGAAGCTGAcg                         >  1:107928/1‑148 (MQ=255)
                                                                                                                           gCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCtt                                                  >  3:269533/1‑120 (MQ=255)
                                                                                                                           gCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCtt                                                  <  4:269533/120‑1 (MQ=255)
                                                                                                                             acaGCGGCAACACGCCGGTCGCCAGTGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGgatga                                    <  7:118923/132‑1 (MQ=255)
                                                                                                                             acaGCGGCAACACGCCGGTCGCCAGTGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGgatga                                    >  8:118923/1‑132 (MQ=255)
                                                                                                                              caGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCTGACATTgcggc                  <  7:164768/149‑1 (MQ=255)
                                                                                                                               aGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCTGACATTgcggcg                 <  4:81656/149‑1 (MQ=255)
                                                                                                                               aGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCTGACATTgcggcg                 <  7:305923/149‑1 (MQ=255)
                                                                                                                                 cGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCTGACATTGCGGCGGa               <  8:219295/149‑1 (MQ=255)
                                                                                                                                   gCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCTGACATTGCGGCGGATa             >  5:158862/1‑149 (MQ=255)
                                                                                                                                        acGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCTGACATTGCGGCGGATACGtt         <  5:246355/148‑1 (MQ=255)
                                                                                                                                        acGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCTGACATTGCGGCGGATACGtt         >  6:246355/1‑148 (MQ=255)
                                                                                                                                         cGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCTGACATTGCGGCGGATACGTTTg       <  3:164389/149‑1 (MQ=255)
                                                                                                                                          gCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCAAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCTGACATTGCGGCGGATACGttt        <  8:200702/147‑1 (MQ=255)
                                                                                                                                            cGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCTGACATTGCGGCGGa               <  4:93124/138‑1 (MQ=255)
                                                                                                                                            cGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCTGACATTGCGGCGGa               >  3:93124/1‑138 (MQ=255)
                                                                                                                                              gTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCTGACATTGCGGCGGATACGTTTGACCTg  <  6:303344/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GCGGGTCTGCAGCACCTCACCCTCGGCATCGCCGGTGTTGGTAGTCTGCCCCATCTGCATGACGGTTTCGTAGCCCTTGTCGGAATCCAGCAAGTACTGCGAAAACTTGGTCGCTTCACCGAAGCACAGCGGCAACACGCCGGTCGCCAGCGGGTCGAGGCTACCGGTATGGCCGGCCTTTTCCGCATTCAGCAGCCAACGCACTTTCTGCAAGGCGGCATTGGACGTGAAGCCCAGCGGCTTGTCGAGCAGGATGATGCCGCTGACATTGCGGCGGATACGTTTGACCTG  >  NC_002947/5355405‑5355695

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