Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,728,840 A→G 8.7% intergenic (‑203/+121) opgG ← / ← dtd glucans biosynthesis protein G/D‑tyrosyl‑tRNA(Tyr) deacylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,728,8400AG 8.7% 192.8 / 7.5 69intergenic (‑203/+121)opgG/dtdglucans biosynthesis protein G/D‑tyrosyl‑tRNA(Tyr) deacylase
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base G (4/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

AATCACGAGCAACCCTCGCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAATGTTACCGG  >  NC_002947/5728692‑5728987
                                                                                                                                                    |                                                                                                                                                   
aaTCACGAGCAACCCTCGCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTa                                                                                                                                                     >  2:129048/1‑149 (MQ=255)
    aCGAGCAACCCTCGCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGaa                                                                                                                                                 >  3:77092/1‑149 (MQ=255)
         cAACCCTCGCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAg                                                                                                                                            <  1:145729/149‑1 (MQ=255)
          aaCCCTCGCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAgg                                                                                                                                           >  8:146159/1‑149 (MQ=255)
             ccTCGCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTc                                                                                                                                         <  8:188549/148‑1 (MQ=255)
             ccTCGCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTc                                                                                                                                         >  7:188549/1‑148 (MQ=255)
                 gCCGAAAACAAATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATg                                                                                                                                     <  1:49708/148‑1 (MQ=255)
                  ccGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCa                                                                                                                                   <  4:128238/149‑1 (MQ=255)
                       aaCAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAgcgcg                                                                                                                              >  3:342525/1‑149 (MQ=255)
                          aGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGa                                                                                                                           >  7:263786/1‑149 (MQ=255)
                               aTGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATc                                                                                                                      <  8:261831/149‑1 (MQ=255)
                                   aaCCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGa                                                                                                                           <  5:130325/140‑1 (MQ=255)
                                   aaCCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGa                                                                                                                           >  6:130325/1‑140 (MQ=255)
                                   aaCCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCcag                                                                                                                   <  7:326649/148‑1 (MQ=255)
                                   aaCCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCcag                                                                                                                   >  5:35665/1‑148 (MQ=255)
                                      cAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAgg                                                                                                               >  3:364884/1‑149 (MQ=255)
                                       aaCGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGCACCGGACCTACGAACGGAGGTCGATGCCGCGCGCGATTATCGAGCTGGt                                                                                                              >  4:284589/1‑149 (MQ=255)
                                          ggCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGc                                                                                                           <  3:191407/149‑1 (MQ=255)
                                           gCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCg                                                                                                          <  8:88683/149‑1 (MQ=255)
                                            ccAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGt                                                                                                         <  1:129048/149‑1 (MQ=255)
                                                  aCAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAAc                                                                                                   <  2:69668/149‑1 (MQ=255)
                                                    aaaTGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCGATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGc                                                                                                 <  6:298636/149‑1 (MQ=255)
                                                         ccGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCa                                                                                             <  4:364884/148‑1 (MQ=255)
                                                             gTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTg                                                                                        <  4:221896/149‑1 (MQ=255)
                                                              ttGGCGAGCACTGTTCGGACCCCGAGCGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTgg                                                                                       <  4:100941/149‑1 (MQ=255)
                                                               tGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGa                                                                                      <  7:146159/149‑1 (MQ=255)
                                                                 gCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAg                                                                                     >  3:359013/1‑148 (MQ=255)
                                                                 gCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAg                                                                                     >  2:77730/1‑148 (MQ=255)
                                                                 gCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGc                                                                                    <  6:35665/149‑1 (MQ=255)
                                                                         tGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTAtt                                                                             >  6:91872/1‑148 (MQ=255)
                                                                             cGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATc                                                                         <  4:342525/148‑1 (MQ=255)
                                                                              ggACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCg                                                                        >  4:300734/1‑148 (MQ=255)
                                                                                aCCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTAc                                                                     <  3:276539/149‑1 (MQ=255)
                                                                                             gTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAgg                                                                                                               >  1:126831/1‑94 (MQ=255)
                                                                                             gTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAgg                                                                                                               <  2:126831/94‑1 (MQ=255)
                                                                                             gTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTg                                                          <  1:27728/147‑1 (MQ=255)
                                                                                             gTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTg                                                          >  2:27728/1‑147 (MQ=255)
                                                                                             gTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGtt                                                        >  8:323105/1‑149 (MQ=255)
                                                                                                  gATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCt                                                   >  8:92174/1‑149 (MQ=255)
                                                                                                    ttCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCtgt                                                 >  5:335377/1‑149 (MQ=255)
                                                                                                       cccAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGaa                                              >  8:17742/1‑149 (MQ=255)
                                                                                                       cccAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGaa                                              >  5:76265/1‑149 (MQ=255)
                                                                                                       cccAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGa                                               >  3:82764/1‑148 (MQ=255)
                                                                                                        ccAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGaac                                             <  4:157157/149‑1 (MQ=255)
                                                                                                         cAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCACGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGaac                                             >  2:240971/1‑148 (MQ=255)
                                                                                                          aaaCGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGaacaa                                           >  6:306800/1‑149 (MQ=255)
                                                                                                            aCGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACg                                         <  5:91872/149‑1 (MQ=255)
                                                                                                                   cAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGtttt                                   <  7:17742/148‑1 (MQ=255)
                                                                                                                   cAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGAGATTATCCCGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTATCCTGTGAACAAGGGGttttt                                  >  6:225363/1‑149 (MQ=255)
                                                                                                                    aGACAATGCTCGAGTCAGAGTGAAACGTACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGtttttt                                 <  8:263786/149‑1 (MQ=255)
                                                                                                                    aGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGtttttt                                 <  4:359013/149‑1 (MQ=255)
                                                                                                                    aGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGtttttt                                 <  4:82764/149‑1 (MQ=255)
                                                                                                                    aGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCACGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGtttttt                                 <  1:240971/149‑1 (MQ=255)
                                                                                                                         aTGCTCGAGTCAGAGTGAAACTGACCTACGAACGCAGGTCCATGCATCGCTCGATTATCCAGCAGGATGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCt                            <  5:285237/149‑1 (MQ=255)
                                                                                                                            cTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCa                          <  3:300734/148‑1 (MQ=255)
                                                                                                                               gAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATAtt                      >  3:192725/1‑149 (MQ=255)
                                                                                                                               gAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATAtt                      >  3:187583/1‑149 (MQ=255)
                                                                                                                                    agagTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGtttttt                                 >  7:34425/1‑133 (MQ=255)
                                                                                                                                    agagTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGtttttt                                 <  8:34425/133‑1 (MQ=255)
                                                                                                                                         tgcatggACCTGCGTTCGTAGGTCCATGTAGCGCGCGATTATACAGAAGGTTGCGTGACAACGCATCAGGATGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAAAGGGTTTTTTGCCGTCATATTTGTAACATa             >  6:90204/6‑148 (MQ=255)
                                                                                                                                           aaCGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAAt          <  7:92174/149‑1 (MQ=255)
                                                                                                                                           aaCGGACCAACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGTTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAAt          >  7:138543/1‑149 (MQ=255)
                                                                                                                                            atggACCTGCGTTCGTAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAAt          >  2:203079/3‑148 (MQ=255)
                                                                                                                                            atggACCTGCGTTCGTAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAAt          <  1:203079/146‑1 (MQ=255)
                                                                                                                                            atggACCTGCGTTCGTAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAATg         >  5:225094/3‑149 (MQ=255)
                                                                                                                                             tggACCTGCGTTCGTAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAAt          >  4:79865/2‑147 (MQ=255)
                                                                                                                                             tggACCTGCGTTCGTAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAAt          <  3:79865/146‑1 (MQ=255)
                                                                                                                                               gACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAATGTTa      >  5:284493/1‑149 (MQ=255)
                                                                                                                                                 ccTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAATGTTAcc    <  3:48059/149‑1 (MQ=255)
                                                                                                                                                   tACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAATGTTACCgg  >  1:49247/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AATCACGAGCAACCCTCGCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAATGTTACCGG  >  NC_002947/5728692‑5728987

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: