Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,531,058 A→C 7.5% R421R (CGT→CGG braE ← high‑affinity branched‑chain amino acid ABC transporter permease BraE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,531,0580AC 7.5% 219.3 / 7.9 80R421R (CGT→CGGbraEhigh‑affinity branched‑chain amino acid ABC transporter permease BraE
Reads supporting (aligned to +/- strand):  ref base A (37/37);  new base C (3/3);  total (40/40)
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.05e-01

GGCAAATATCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGAGCCCATGC  >  NC_002947/5530909‑5531206
                                                                                                                                                     |                                                                                                                                                    
ggCAAATATCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATa                                                                                                                                                      >  2:2243991/1‑150 (MQ=255)
 gCAAATATCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCGGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATAc                                                                                                                                                     >  1:1185905/1‑150 (MQ=255)
       atCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCa                                                                                                                                               <  2:670948/150‑1 (MQ=255)
        tCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAg                                                                                                                                              <  1:2367006/150‑1 (MQ=255)
         cttgTTGCGCCCGACCTTCAGTCTGACTTCGCTTAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCTACAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGTGTCACACCGGTACGGCTGATACGGATCAGc                                                                                                                                             <  2:930718/147‑1 (MQ=255)
          tGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGcc                                                                                                                                            >  1:363865/1‑150 (MQ=255)
                       ccTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCa                                                                                                                               <  2:198102/150‑1 (MQ=255)
                           cAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGatc                                                                                                                           <  2:1474308/150‑1 (MQ=255)
                                 gACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCGGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGatcatcatc                                                                                                                     <  2:1185905/150‑1 (MQ=255)
                                 gACGTCGCTGAGCGCCGTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCTGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGatcatcatc                                                                                                                     <  1:368658/150‑1 (MQ=255)
                                  aCGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGAtcatcatca                                                                                                                    >  1:400231/1‑150 (MQ=255)
                                   cGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATcatcag                                                                                                                   <  1:145379/150‑1 (MQ=255)
                                         tGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATcatcagcaccat                                                                                                             >  2:2254287/1‑150 (MQ=255)
                                           ggcgcCGTGAGGCCACCGAACTGGATCATCAGGTTGTCGACCGAGAGAAGGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATcatcagcaccatca                                                                                                           <  2:402645/149‑1 (MQ=255)
                                                cTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATGAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACg                                                                                                      >  1:691919/1‑150 (MQ=255)
                                                cTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACg                                                                                                      >  2:2373704/1‑150 (MQ=255)
                                                 ttGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGc                                                                                                     <  2:363865/150‑1 (MQ=255)
                                                     tGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTAGTACGAATCATCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGaa                                                                                                 >  2:1270894/1‑150 (MQ=255)
                                                      gCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAAc                                                                                                <  1:1612534/150‑1 (MQ=255)
                                                      gCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAAc                                                                                                <  1:2373704/150‑1 (MQ=255)
                                                        cACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGCTACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAAcag                                                                                              >  1:1836132/1‑150 (MQ=255)
                                                          ccGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCAcgc                                                                                                                                       >  2:19124/1‑107 (MQ=255)
                                                          ccGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCAcgc                                                                                                                                       <  1:19124/107‑1 (MQ=255)
                                                          ccGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGAACAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAAcagca                                                                                            <  1:2243991/150‑1 (MQ=255)
                                                             aaCTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAgcagca                                                                                         <  2:1836132/150‑1 (MQ=255)
                                                                tGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCAccc                                                                                      >  2:1910721/1‑150 (MQ=255)
                                                                 gCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCg                                                                                     <  1:2254755/150‑1 (MQ=255)
                                                                    tcatcaGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATCCGTATCAGccg                                                                                                                                           >  2:1553438/1‑92 (MQ=255)
                                                                    tcatcaGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATCCGTATCAGccg                                                                                                                                           <  1:1553438/93‑2 (MQ=255)
                                                                       tcaGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAAcagca                                                                                            <  2:1483333/137‑1 (MQ=255)
                                                                       tcaGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAAcagca                                                                                            >  1:1483333/1‑137 (MQ=255)
                                                                          ggTTGTCGACCGAGAGAATGATATCGTCGCTCATGGTGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGcc                                                                                                    >  1:2242329/1‑126 (MQ=255)
                                                                          ggTTGTCGACCGAGAGAATGATATCGTCGCTCATGGTGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGcc                                                                                                    <  2:2242329/126‑1 (MQ=255)
                                                                            ttGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCg                                                                              <  2:2299866/146‑1 (MQ=255)
                                                                            ttGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACCCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCg                                                                              >  1:2299866/1‑146 (MQ=255)
                                                                               tCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGc                                                                       <  2:660861/150‑1 (MQ=255)
                                                                                   ccGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTgcg                                                                   >  1:1033882/1‑150 (MQ=255)
                                                                                    cGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCg                                                                                     >  2:316892/1‑131 (MQ=255)
                                                                                    cGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCg                                                                                     <  1:316892/131‑1 (MQ=255)
                                                                                     gagagaATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCGCCCGGTACTCGGAGAAGCTGCgca                                                                 <  1:2254287/150‑1 (MQ=255)
                                                                                     gagagaATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCTGATAAGCTGCgca                                                                 >  1:1696982/1‑150 (MQ=255)
                                                                                                atCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACGGCAGCACCCGGTACTCGTAGAAGCTGCGCAGCAGTTCCGgc                                                      <  2:691919/150‑1 (MQ=255)
                                                                                                   gtcgCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATcatcagcaccat                                                                                                             >  2:1163383/1‑92 (MQ=255)
                                                                                                   gtcgCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATcatcagcaccat                                                                                                             <  1:1163383/92‑1 (MQ=255)
                                                                                                   gtcgCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGcc                                                   >  1:2020676/1‑150 (MQ=255)
                                                                                                    tcgCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATcatcagca                                                                                                                 >  2:714045/1‑87 (MQ=255)
                                                                                                    tcgCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATcatcagca                                                                                                                 <  1:714045/87‑1 (MQ=255)
                                                                                                    tcgCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCg                                                                                     >  2:1945408/1‑115 (MQ=255)
                                                                                                    tcgCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCg                                                                                     <  1:1945408/115‑1 (MQ=255)
                                                                                                      gCTCATGGCGCCACTCCTCTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACACCGAACAGCAGCACCCGATACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACg                                                >  1:689317/1‑150 (MQ=255)
                                                                                                                 cACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGaa                                     <  1:101769/150‑1 (MQ=255)
                                                                                                                      cTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGga                                                                            <  1:1125097/106‑1 (MQ=255)
                                                                                                                      cTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGga                                                                            >  2:1125097/1‑106 (MQ=255)
                                                                                                                      cTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAAcgccg                                >  2:1085624/1‑150 (MQ=255)
                                                                                                                       tttACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAAcgccgc                               >  1:856587/1‑150 (MQ=255)
                                                                                                                            gTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCgca                                                                 >  1:1063111/1‑111 (MQ=255)
                                                                                                                            gTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCgca                                                                 <  2:1063111/111‑1 (MQ=255)
                                                                                                                               gggTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACca                       >  2:1671191/1‑150 (MQ=255)
                                                                                                                                ggTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACcac                      >  1:1975420/1‑150 (MQ=255)
                                                                                                                                    acacCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCg                                                                                     <  1:1409139/83‑1 (MQ=255)
                                                                                                                                    acacCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCg                                                                                     >  2:1409139/1‑83 (MQ=255)
                                                                                                                                      acCGGTACGGCTGATACGGATCAGCGCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCAc                >  2:2035438/1‑150 (MQ=255)
                                                                                                                                      acCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGAc                         >  2:608318/1‑141 (MQ=255)
                                                                                                                                      acCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGAc                         <  1:608318/141‑1 (MQ=255)
                                                                                                                                          gTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGc                                                                                                     <  2:679277/61‑1 (MQ=255)
                                                                                                                                          gTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGc                                                                                                     >  1:679277/1‑61 (MQ=255)
                                                                                                                                            tgggCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGa          >  1:948387/3‑150 (MQ=255)
                                                                                                                                             gggCTGATCCGTATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGt             >  2:1929567/2‑146 (MQ=255)
                                                                                                                                             gggCTGATCCGTATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGt             <  1:1929567/145‑1 (MQ=255)
                                                                                                                                              ggCTGAAACCGCACAGCACACACGATCACCACATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGACGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGAGc        >  2:149847/1‑150 (MQ=255)
                                                                                                                                               gCTGATCCGTATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGa                          >  2:1923569/1‑131 (MQ=255)
                                                                                                                                               gCTGATCCGTATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGa                          <  1:1923569/131‑1 (MQ=255)
                                                                                                                                               gCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCGCCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGa                          <  1:837848/131‑1 (MQ=255)
                                                                                                                                               gCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCGCCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGa                          >  2:837848/1‑131 (MQ=255)
                                                                                                                                                 tGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGAGCCCa     >  1:1979476/1‑150 (MQ=255)
                                                                                                                                                 tGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGAGCCCa     <  2:948387/150‑1 (MQ=255)
                                                                                                                                                  gATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGAGCCCAt    <  1:1974621/150‑1 (MQ=255)
                                                                                                                                                   aTACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGAGCCCATg   >  2:579464/1‑150 (MQ=255)
                                                                                                                                                    tACGGATCAGCCCACGCGGTCGCCAGAtcatcatca                                                                                                                    <  2:2374913/36‑1 (MQ=255)
                                                                                                                                                    tACGGATCAGCCCACGCGGTCGCCAGAtcatcatca                                                                                                                    >  1:2374913/1‑36 (MQ=255)
                                                                                                                                                    tACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACATCCTCGGTCGAGCCCATGc  <  2:2020676/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCAAATATCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGAGCCCATGC  >  NC_002947/5530909‑5531206

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