Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,758,465 G→A 5.8% H233H (CAC→CAT PP_2412 ← ParA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,758,4650GA 5.8% 209.3 / 2.1 69H233H (CAC→CATPP_2412ParA family protein
Reads supporting (aligned to +/- strand):  ref base G (39/26);  new base A (2/2);  total (41/28)
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.65e-01

TGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGATTGACTACCACACCTTCCACCACCAGCGCGGGG  >  NC_002947/2758323‑2758612
                                                                                                                                              |                                                                                                                                                   
tGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGAt                                                                                                                                               >  8:312359/1‑149 (MQ=255)
   ttCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTcg                                                                                                                                            <  3:303518/149‑1 (MQ=255)
      gACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGATGCGATTcgcg                                                                                                                                          >  8:150705/1‑148 (MQ=255)
      gACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGATGCGATTcgag                                                                                                                                          >  8:145584/1‑146 (MQ=255)
        ccGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCAt                                                                                                                                       <  8:302817/149‑1 (MQ=255)
           gTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTc                                                                                                                                             <  5:260315/140‑1 (MQ=255)
           gTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTc                                                                                                                                             >  6:260315/1‑140 (MQ=255)
               agcagGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATc                                                                                                                                <  1:13755/149‑1 (MQ=255)
                         aGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGTTCGAACTGCt                                                                                                                        >  7:182758/1‑147 (MQ=255)
                           ccGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAgg                                                                                                                    >  5:58320/1‑149 (MQ=255)
                                     cccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGAt                                                                                                                                               <  6:47200/112‑1 (MQ=255)
                                     cccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGAt                                                                                                                                               >  5:47200/1‑112 (MQ=255)
                                     cccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAg                                                                                                                     <  4:295395/138‑1 (MQ=255)
                                     cccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAg                                                                                                                     >  3:295395/1‑138 (MQ=255)
                                      ccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGgcag                                                                                                         >  2:261818/1‑149 (MQ=255)
                                      ccAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGgcag                                                                                                         >  1:25174/1‑149 (MQ=255)
                                       cAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGgcag                                                                                                         >  5:321550/1‑148 (MQ=255)
                                        tgcGGCCTGACAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGgcagca                                                                                                       >  6:100740/3‑149 (MQ=255)
                                        aaCGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGgcagca                                                                                                       >  1:142967/1‑149 (MQ=255)
                                          cGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCAcc                                                                                                     >  6:129835/1‑149 (MQ=255)
                                           ggCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCg                                                                                                    >  6:73293/1‑149 (MQ=255)
                                               tGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCAt                                                                                                >  6:276962/1‑149 (MQ=255)
                                                 tCACGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCAt                                                                                                <  5:276962/147‑1 (MQ=255)
                                                 acAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCAt                                                                                                <  5:100740/146‑1 (MQ=255)
                                                   aGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCt                                                                                            >  1:202021/1‑149 (MQ=255)
                                                     gCCGCACGCTCAAGCACAACCAGCAAATCGATGAATTCGAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCATc                                                                                          >  8:176707/1‑149 (MQ=255)
                                                      ccGCACGCTCAAGCACATCCAGCAAATCGATGAATTCGAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCg                                                                                         <  8:182758/149‑1 (MQ=255)
                                                        gCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGtt                                                                                       <  2:25174/149‑1 (MQ=255)
                                                            gCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTgcg                                                                                    >  6:319931/1‑148 (MQ=255)
                                                                aaGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGgcagca                                                                                                       <  6:207824/125‑1 (MQ=255)
                                                                aaGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGgcagca                                                                                                       >  5:207824/1‑125 (MQ=255)
                                                                aaGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGtt                                                                               >  6:128839/1‑149 (MQ=255)
                                                                  gCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGgca                                                                                                          >  2:315728/1‑120 (MQ=255)
                                                                  gCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGgca                                                                                                          <  1:315728/120‑1 (MQ=255)
                                                                  gCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGa                                                                             >  6:141543/1‑149 (MQ=255)
                                                                         cAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGAcc                                                                       >  1:11103/1‑148 (MQ=255)
                                                                          aGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGAcc                                                                       >  6:306891/1‑147 (MQ=255)
                                                                            cAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAgg                                                                   >  3:242886/1‑149 (MQ=255)
                                                                               aTCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCt                                                                >  3:180299/1‑149 (MQ=255)
                                                                                  gATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAg                                                                                 >  7:46099/1‑129 (MQ=255)
                                                                                   aTGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTg                                                            >  5:88625/1‑149 (MQ=255)
                                                                                       aTTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATCGTGCGATTCGCGCATCTTGATCCAACTGCTCAGGTACACCGGCAGCACCGGCATCCCTTCGTTGCGCAGTTGATCGACCAACGTCTGGTGCAgc                                                        >  2:97168/1‑149 (MQ=255)
                                                                                           cAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCg                                                    >  5:184676/1‑149 (MQ=255)
                                                                                            aGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGt                                                   <  2:142967/149‑1 (MQ=255)
                                                                                              ggCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCTTGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTgc                                                 <  1:97168/149‑1 (MQ=255)
                                                                                               gCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGTCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTgcg                                                >  7:42046/1‑149 (MQ=255)
                                                                                                    tttATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCgg                                           <  1:261818/149‑1 (MQ=255)
                                                                                                      tATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCa                                                                                  <  1:287488/108‑1 (MQ=255)
                                                                                                      tATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCa                                                                                  >  2:287488/1‑108 (MQ=255)
                                                                                                         gCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAAc                                      >  3:82129/1‑149 (MQ=255)
                                                                                                           cgcgGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACt                                     >  1:376460/1‑148 (MQ=255)
                                                                                                            gcggcgCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACTCCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGa                                   <  8:23332/149‑1 (MQ=255)
                                                                                                            gcggcgCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGa                                   >  8:115076/1‑149 (MQ=255)
                                                                                                              ggcgCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGAtt                                 <  4:180299/149‑1 (MQ=255)
                                                                                                              ggcgCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGAtt                                 <  5:319931/149‑1 (MQ=255)
                                                                                                                  ccAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGATTGAc                              >  4:21793/1‑148 (MQ=255)
                                                                                                                       tGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGATTGACTACcac                        <  5:128839/149‑1 (MQ=255)
                                                                                                                         cACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGATTGACTACcaca                       >  3:65360/1‑148 (MQ=255)
                                                                                                                          aCCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGATTGACTACCACAcc                     >  2:364639/1‑149 (MQ=255)
                                                                                                                          aCCAATGGCACCGAGGCATGATGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGATTGACTACCACAcc                     <  7:145584/149‑1 (MQ=255)
                                                                                                                          aCCAATGGCACCGAGGCATGATGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGATTGACTACCACAcc                     <  7:150705/149‑1 (MQ=255)
                                                                                                                           ccAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGATTGACTACCACACCt                    >  3:32303/1‑149 (MQ=255)
                                                                                                                             aaTGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGATTGACTACCACACCTTc                  <  7:312359/149‑1 (MQ=255)
                                                                                                                              aTGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGATTGACTACCACACCTTcc                 <  5:129835/149‑1 (MQ=255)
                                                                                                                                   aCCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGATTGACTACCACACCTTccaccac            <  5:73293/149‑1 (MQ=255)
                                                                                                                                    ccGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGTCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGATTGACTACCACACCTTccaccacc           <  8:42046/149‑1 (MQ=255)
                                                                                                                                         gCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCgg                                           <  3:204673/112‑1 (MQ=255)
                                                                                                                                         gCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCgg                                           >  4:204673/1‑112 (MQ=255)
                                                                                                                                          cATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGATTGACTACCACACCTTCCACCACCAGcgcg     <  6:88625/149‑1 (MQ=255)
                                                                                                                                              gTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTc                                                                                          >  3:301536/1‑60 (MQ=255)
                                                                                                                                              gTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTc                                                                                          <  4:301536/60‑1 (MQ=255)
                                                                                                                                              gTGCGATTCGCGCATCTTGATAGAACTGCTGAGGTTGACCGGCAGCAACGGCATGCCTTCGTTGAGCAGTTGATCGACCAAGGTCTGGTGCAGCTCCGTTCGCCCGGCGAACTGATTGACTACCACACATTCCACCACCAgctagttg  <  4:32303/148‑7 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TGATTCGACCGGTGCAGCAGGTTACAGCCGGAGTGCCCCCAACGGCCTGTCAGGCCGCACGCTCAAGCACATCCAGCAAATCGATGAATTCCAGGGCCAGTTTATGCCGCGGCGCCAAATGCACCAATGGCACCGAGGCATGGTGCGATTCGCGCATCTTGATCGAACTGCTCAGGTACACCGGCAGCACCGGCATGCCTTCGTTGCGCAGTTGATCGACCAAGGTCTGGTGCAGCGCCGTGCGCCCGGCGAACTGATTGACTACCACACCTTCCACCACCAGCGCGGGG  >  NC_002947/2758323‑2758612

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: