Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,641,139 G→C 6.1% G251R (GGG→CGG)  PP_4951 → EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,641,1390GC 6.1% 206.5 / 3.5 65G251R (GGG→CGG) PP_4951EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base G (35/26);  new base C (2/2);  total (37/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 = 9.92e-01

ACGGCCCTGGTTTCCTGCTGTTCGGCGCAGCCATGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGG  >  NC_002947/5640990‑5641281
                                                                                                                                                     |                                                                                                                                              
aCGGCCCTGGTTTCCTGCTGTTCGGCGCAGCCATGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCg                                                                                                                                                >  2:173619/1‑150 (MQ=255)
        ggTTTCCTGCTGTTCGTCGCAGCCATGGTACTGATCACCATTGGCATGCGAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGAc                                                                                                                                        >  2:644805/1‑150 (MQ=255)
        ggTTTCCTGCTGTTCGGCGCAGCCATGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCCGGcgcgca                                                                                                                                        <  1:686270/150‑3 (MQ=255)
         gTTTCCTGCTGTTCGGCGCAGCCATGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGAcc                                                                                                                                       >  1:881246/1‑150 (MQ=255)
           ttCCTGCTGTTCGGCGCAGCCATGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCg                                                                                                                                     <  2:319644/150‑1 (MQ=255)
               tgctgTTCGGCGCAGCCATGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGTGCGACCCGGAGc                                                                                                                                 >  1:798567/1‑150 (MQ=255)
               tgctgTTCGGCGCAGCCATGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGc                                                                                                                                 >  1:1090872/1‑150 (MQ=255)
               tgctgTTCGGCGCAGCCATGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGc                                                                                                                                 >  1:1087101/1‑150 (MQ=255)
               tgctgTTCGGCGCAGCCATGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGc                                                                                                                                 >  2:1124295/1‑150 (MQ=255)
               tgctgTTCGGCGCAGCCATGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGc                                                                                                                                 >  2:327348/1‑150 (MQ=255)
                    ttCGGCGCAGCCATGGTACTGATCGCCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTg                                                                                                                            >  2:753619/1‑150 (MQ=255)
                         cgcAGCCATGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTtctc                                                                                                                       <  2:295280/150‑1 (MQ=255)
                          gcAGCCATGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGgc                                                                                                                                             <  1:349476/127‑1 (MQ=255)
                          gcAGCCATGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGgc                                                                                                                                             >  2:349476/1‑127 (MQ=255)
                              ccATGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCAACTCGCTCTCTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCaa                                                                                                                  <  1:346021/150‑1 (MQ=255)
                                aTGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGc                                                                                                                >  2:803286/1‑150 (MQ=255)
                                     aCTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGtt                                                                                                           <  1:173619/150‑1 (MQ=255)
                                      cTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTc                                                                                                          >  2:573659/1‑150 (MQ=255)
                                      cTGATCAACATTGCCATGGACGACCTGGGGGTCCTGCACCGGCCGCTCTTTCGGGTAATGGTGCTGCTGTTCGCCGGCATGGACTGCATGAACGCGGAGTGGGGGAGGCGCGTGTGCTACCCTGTGGCTCTGTTCTCGCCAAGCCTGTTc                                                                                                          >  2:514442/1‑150 (MQ=255)
                                        gATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCcg                                                                                                        <  2:153945/150‑1 (MQ=255)
                                         aTCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCg                                                                                                                               >  2:403548/1‑126 (MQ=255)
                                         aTCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCg                                                                                                                               <  1:403548/126‑1 (MQ=255)
                                         aTCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCcgc                                                                                                       >  1:595831/1‑150 (MQ=255)
                                           cACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCgcgt                                                                                                     <  1:327348/150‑1 (MQ=255)
                                              cATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGgc                                                                                                                                             <  1:821292/107‑1 (MQ=255)
                                              cATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGgc                                                                                                                                             >  2:821292/1‑107 (MQ=255)
                                              cATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCgcgtgcg                                                                                                  >  2:984159/1‑150 (MQ=255)
                                                 tGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTAc                                                                                               <  2:859114/150‑1 (MQ=255)
                                                   ccATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACtt                                                                                             >  2:662863/1‑150 (MQ=255)
                                                     aTGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGa                                                                                                                                   <  2:1077118/110‑1 (MQ=255)
                                                     aTGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGa                                                                                                                                   >  1:1077118/1‑110 (MQ=255)
                                                     aTGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTc                                                                                           <  2:107937/150‑1 (MQ=255)
                                                           ggCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGATGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATc                                                                                     <  1:314155/150‑1 (MQ=255)
                                                            gCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCg                                                                                    <  1:101173/150‑1 (MQ=255)
                                                                 ggCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGtt                                                                                                           <  2:418511/122‑1 (MQ=255)
                                                                 ggCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGtt                                                                                                           >  1:418511/1‑122 (MQ=255)
                                                                 ggCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCAtt                                                                               >  2:741896/1‑149 (MQ=255)
                                                                     aGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCg                                                                           <  2:881246/150‑1 (MQ=255)
                                                                      gCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATTCTCgg                                                                          <  1:741896/150‑1 (MQ=255)
                                                                                      tttcGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTAc                                                                                               <  1:784547/110‑1 (MQ=255)
                                                                                      tttcGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTAc                                                                                               >  2:784547/4‑113 (MQ=255)
                                                                                           ggTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTggg                                                     >  1:698000/1‑150 (MQ=255)
                                                                                              aaTGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAg                                                  <  2:926926/150‑1 (MQ=255)
                                                                                                 gTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGACCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAgcgg                                               >  2:93016/1‑150 (MQ=255)
                                                                                                    gctgctGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGgcgc                                            >  1:546673/1‑150 (MQ=255)
                                                                                                      tgctgTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGgcgcgc                                          >  1:687348/1‑150 (MQ=255)
                                                                                                          gTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGcc                                      >  1:179064/1‑150 (MQ=255)
                                                                                                              gCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTc                                  >  2:1020713/1‑150 (MQ=255)
                                                                                                              gCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTc                                  <  2:546673/150‑1 (MQ=255)
                                                                                                                  gCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTgg                              >  1:253593/1‑150 (MQ=255)
                                                                                                                       gcctgcctGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTTGGCAGCGGCGCGCTGCCGTTCCTGGTGCCa                         >  1:393518/1‑150 (MQ=255)
                                                                                                                         ctgcctgGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGcc                                      >  2:778894/1‑135 (MQ=255)
                                                                                                                         ctgcctgGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGcc                                      <  1:778894/135‑1 (MQ=255)
                                                                                                                              tgGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATtgctgc                  <  1:644805/150‑1 (MQ=255)
                                                                                                                                     gTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGt           <  1:573659/150‑1 (MQ=255)
                                                                                                                                       aCTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTg         >  2:94064/1‑150 (MQ=255)
                                                                                                                                        cTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCt                                <  2:1006419/126‑1 (MQ=255)
                                                                                                                                        cTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCt                                >  1:1006419/1‑126 (MQ=255)
                                                                                                                                           gggtcgcgcCCGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATtgctgct                 >  1:1008746/5‑138 (MQ=255)
                                                                                                                                           gggtcgcgcCCGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATtgctgct                 <  2:1008746/134‑1 (MQ=255)
                                                                                                                                           gggtcgagcCCGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCt     >  1:539906/8‑150 (MQ=255)
                                                                                                                                           gCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGGACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCt     >  1:415763/1‑150 (MQ=255)
                                                                                                                                             tGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTgg   >  1:842431/1‑150 (MQ=255)
                                                                                                                                             tGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACCTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGt                             <  2:529943/124‑1 (MQ=255)
                                                                                                                                             tGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACCTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGt                             >  1:529943/1‑124 (MQ=255)
                                                                                                                                              gcgcgcCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTggg  <  1:753619/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
ACGGCCCTGGTTTCCTGCTGTTCGGCGCAGCCATGGTACTGATCACCATTGCCATGGAAGGCCTGGGCGAGCTGCACCTGCCGCACTTGCGGGTAATGTTGCTGCTGTTCGCCGGCATGGCCTGCCTGGCCGCGTACTGGCTGCGCGCCGGGCGCGACCCGGAGCCGCTGTTCTCGCCAAGCCTGTTCCGCGTGCGTACTTTCGCCATCGGCATCCTCGGCAACCTGTTCGCCCGCCTGGGCAGCGGCGCGCTGCCGTTCCTGGTGCCATTGCTGCTGCAGGTTGCGCTGGG  >  NC_002947/5640990‑5641281

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