Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,515,106 T→C 6.6% E307G (GAG→GGG) ‡ cbpA ← curved DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,515,1060TC 6.6% 171.0 / 2.5 60E307G (GAG→GGG) ‡cbpAcurved DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (25/31);  new base C (2/2);  total (27/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 = 8.59e-01

TTCGACAATGCCGTGCTCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGCTCTGGCTG  >  NC_002947/5514959‑5515252
                                                                                                                                                   |                                                                                                                                                  
ttCGACAATGCCGTGCTCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTtctc                                                                                                                                                   <  3:338681/149‑1 (MQ=255)
            gagCTCTACTATTTCGATCACGCATTCATCCGTGACATCTGCTTCCTGACACAGTGTACGGATGTCCAGTTGTACGATCATGGTTCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGt                                                                                                                                       <  4:112488/147‑1 (MQ=255)
               cTCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCCGAGAAAGCtgcgttc                                                                                                                                     <  5:258901/148‑7 (MQ=255)
                  gacgaTTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCaga                                                                                                                                  <  4:222932/148‑1 (MQ=255)
                      aTTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTcg                                                                                                                             >  7:173670/1‑149 (MQ=255)
                          cGATCACGCATTCAGCCGTGACATCGGCATCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCggg                                                                                                                         >  1:200205/1‑149 (MQ=255)
                           gATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCTGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCggg                                                                                                                         <  2:200205/148‑1 (MQ=255)
                            aTCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCa                                                                                                                       <  8:173670/149‑1 (MQ=255)
                            aTCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCa                                                                                                                       <  5:143202/149‑1 (MQ=255)
                               aCGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGa                                                                                                                     <  7:258863/148‑1 (MQ=255)
                                 gCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTc                                                                                                                  <  5:87023/149‑1 (MQ=255)
                                       aGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGa                                                                                                            >  4:204645/1‑149 (MQ=255)
                                          cGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGAcgc                                                                                                         <  1:276540/149‑1 (MQ=255)
                                                         cTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTc                                                                                                                  <  4:280381/125‑1 (MQ=255)
                                                         cTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTc                                                                                                                  >  3:280381/1‑125 (MQ=255)
                                                              acaGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCtt                                                                                       >  1:254944/1‑147 (MQ=255)
                                                                aGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTcc                                                                                                              <  6:127388/122‑1 (MQ=255)
                                                                aGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTcc                                                                                                              >  5:127388/1‑122 (MQ=255)
                                                                  ggTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTg                                                                                 >  4:157985/1‑149 (MQ=255)
                                                                   gTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTgc                                                                                <  7:140044/149‑1 (MQ=255)
                                                                   gTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTAATGATCACTTACTCCCTTGTGTCCTCTGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTgc                                                                                <  5:193859/149‑1 (MQ=255)
                                                                      cGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTgcgc                                                                              <  2:254944/148‑1 (MQ=255)
                                                                             ccAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGt                                                                      >  1:342784/1‑149 (MQ=255)
                                                                             ccAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCAGTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGt                                                                      >  6:283667/1‑149 (MQ=255)
                                                                                   gAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCcgc                                                                >  5:334510/1‑149 (MQ=255)
                                                                                      cGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGct                                                             <  3:204645/149‑1 (MQ=255)
                                                                                            gggTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCt                                                       <  3:157727/149‑1 (MQ=255)
                                                                                              gTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCtgt                                                     <  3:157985/149‑1 (MQ=255)
                                                                                                gttgctCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGaccacc                                                                                         <  1:120061/109‑1 (MQ=255)
                                                                                                gctgctCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGaccacc                                                                                         >  2:120061/1‑111 (MQ=255)
                                                                                                 ctgctCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCtgttt                                                   <  7:210690/148‑1 (MQ=255)
                                                                                                     tCATGATGACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTgcgc                                                                              <  3:134539/117‑1 (MQ=255)
                                                                                                     tCATGATAACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTgcgc                                                                              >  4:134539/1‑117 (MQ=255)
                                                                                                      cATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGc                                             >  1:17144/1‑149 (MQ=255)
                                                                                                       aTGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGc                                             >  8:47395/1‑148 (MQ=255)
                                                                                                        tGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCa                                           <  2:17144/149‑1 (MQ=255)
                                                                                                          aTCTCTTACTCGATTGTGTCCTCGGGTTGAAGGCTGCTTGCTCTGAATGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCTACGCTGGCGGGATGACCACCTTGAGCTGCGCATAATGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAg                                          <  5:264624/148‑1 (MQ=255)
                                                                                                              cTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCt                                     >  3:80167/1‑149 (MQ=255)
                                                                                                              cTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCt                                     >  6:51395/1‑149 (MQ=255)
                                                                                                                   tCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTTGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGGTCTCCCTGTTTGTTCGCCAGGCCCTTGccc                                >  6:75284/1‑149 (MQ=255)
                                                                                                                   tCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGccc                                >  3:157477/1‑149 (MQ=255)
                                                                                                                    ccATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTgc                                                                                <  1:263295/100‑1 (MQ=255)
                                                                                                                    ccATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTgc                                                                                >  2:263295/1‑100 (MQ=255)
                                                                                                                     cATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGccc                                <  5:51395/147‑1 (MQ=255)
                                                                                                                      aTTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGt                             <  3:66800/149‑1 (MQ=255)
                                                                                                                      aTTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGt                             <  2:172185/149‑1 (MQ=255)
                                                                                                                            tCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACg                        >  5:143530/1‑148 (MQ=255)
                                                                                                                            tCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCATCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACg                        >  4:228865/1‑148 (MQ=255)
                                                                                                                              cTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGACAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAg                     <  5:268778/149‑1 (MQ=255)
                                                                                                                                    ttGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTg               >  1:47977/1‑149 (MQ=255)
                                                                                                                                      gAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTgg              <  5:283667/148‑1 (MQ=255)
                                                                                                                                         ccaaGCTTTCCGAGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCt                >  2:349196/4‑143 (MQ=255)
                                                                                                                                         ccaaGCTTTCCGAGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCt                <  1:349196/140‑1 (MQ=255)
                                                                                                                                          caaGCTTTCCGAGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGc          >  5:150701/3‑148 (MQ=255)
                                                                                                                                            aGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAGGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGctct       >  1:234139/1‑149 (MQ=255)
                                                                                                                                            aGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGctct       <  6:334510/149‑1 (MQ=255)
                                                                                                                                             gCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGc                                        <  6:19020/115‑1 (MQ=255)
                                                                                                                                             gCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGc                                        >  5:19020/1‑115 (MQ=255)
                                                                                                                                             gCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGCTCTg      <  4:80167/149‑1 (MQ=255)
                                                                                                                                              cTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGCTCTgg     <  2:342784/149‑1 (MQ=255)
                                                                                                                                               tttCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGCTCTGGc    >  8:304898/1‑149 (MQ=255)
                                                                                                                                                 tctcGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGCTCTGGCt   <  1:144089/148‑1 (MQ=255)
                                                                                                                                                 tctcGCAAAGCTTGGTCCAGAGTTCGCGGCCAGATTCGTCCGACGCGGGCGGCATGACGACCTTGAGCTGCGCAGAAAGGTTGCCGCGCACTCCCTGTGGGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCCCAGGCCGCTCTGGCTg  >  8:267296/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTCGACAATGCCGTGCTCGACGATTTCGATCACGCATTCAGCCGTGACATCGGCTTCCTGACACAGGGTACGCATGTCCAGTTGAACGATCAGGGTGCTGCTCATGATCACTTACTCCATTGTGTCCTCGGGTTGAACGCAGCTTTCTCGGAAAGCTTGGTCCAGAGTTCGCGGGCAGATTCGTCCGACGCGGGCGGCATGACCACCTTGAGCTGCGCATAAAGGTTGCCGCGCTCTCCCTGTTTGTTCGCCAGGCCCTTGCCCGGTACACGCAGGCGCTGGCCGCTCTGGCTG  >  NC_002947/5514959‑5515252

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