Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,531,061 G→T 7.2% I420I (ATC→ATA 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,0610GT 7.2% 245.0 / 5.5 79I420I (ATC→ATAbraEhigh‑affinity branched‑chain amino acid ABC transporter permease BraE
Reads supporting (aligned to +/- strand):  ref base G (35/37);  new base T (3/3);  total (39/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 = 8.72e-01

ATCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGAGCCCATGCCAC  >  NC_002947/5530916‑5531209
                                                                                                                                                 |                                                                                                                                                    
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:2373704/150‑1 (MQ=255)
                                               gCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAAc                                                                                                   <  1:1612534/150‑1 (MQ=255)
                                                 cACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGCTACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAAcag                                                                                                 >  1:1836132/1‑150 (MQ=255)
                                                   ccGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCAcgc                                                                                                                                          <  1:19124/107‑1 (MQ=255)
                                                   ccGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCAcgc                                                                                                                                          >  2:19124/1‑107 (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                                                                                                                                              <  1:1553438/93‑2 (MQ=255)
                                                             tcatcaGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATCCGTATCAGccg                                                                                                                                              >  2:1553438/1‑92 (MQ=255)
                                                                tcaGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAAcagca                                                                                               >  1:1483333/1‑137 (MQ=255)
                                                                tcaGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAAcagca                                                                                               <  2:1483333/137‑1 (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                                                                                        <  1:316892/131‑1 (MQ=255)
                                                                             cGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCg                                                                                        >  2:316892/1‑131 (MQ=255)
                                                                              gagagaATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCGCCCGGTACTCGGAGAAGCTGCgca                                                                    <  1:2254287/150‑1 (MQ=255)
                                                                              gagagaATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCTGATAAGCTGCgca                                                                    >  1:1696982/1‑150 (MQ=255)
                                                                                         atCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACGGCAGCACCCGGTACTCGTAGAAGCTGCGCAGCAGTTCCGgc                                                         <  2:691919/150‑1 (MQ=255)
                                                                                            gtcgCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATcatcagcaccat                                                                                                                <  1:1163383/92‑1 (MQ=255)
                                                                                            gtcgCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATcatcagcaccat                                                                                                                >  2:1163383/1‑92 (MQ=255)
                                                                                            gtcgCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGcc                                                      >  1:2020676/1‑150 (MQ=255)
                                                                                             tcgCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATcatcagca                                                                                                                    <  1:714045/87‑1 (MQ=255)
                                                                                             tcgCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATcatcagca                                                                                                                    >  2:714045/1‑87 (MQ=255)
                                                                                             tcgCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCg                                                                                        <  1:1945408/115‑1 (MQ=255)
                                                                                             tcgCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCg                                                                                        >  2:1945408/1‑115 (MQ=255)
                                                                                               gCTCATGGCGCCACTCCTCTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACACCGAACAGCAGCACCCGATACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACg                                                   >  1:689317/1‑150 (MQ=255)
                                                                                                          cACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGaa                                        <  1:101769/150‑1 (MQ=255)
                                                                                                               cTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGga                                                                               >  2:1125097/1‑106 (MQ=255)
                                                                                                               cTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGga                                                                               <  1:1125097/106‑1 (MQ=255)
                                                                                                               cTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAAcgccg                                   >  2:1085624/1‑150 (MQ=255)
                                                                                                                tttACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAAcgccgc                                  >  1:856587/1‑150 (MQ=255)
                                                                                                                     gTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCgca                                                                    <  2:1063111/111‑1 (MQ=255)
                                                                                                                     gTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCgca                                                                    >  1:1063111/1‑111 (MQ=255)
                                                                                                                        gggTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACca                          >  2:1671191/1‑150 (MQ=255)
                                                                                                                         ggTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACcac                         >  1:1975420/1‑150 (MQ=255)
                                                                                                                             acacCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCg                                                                                        >  2:1409139/1‑83 (MQ=255)
                                                                                                                             acacCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCg                                                                                        <  1:1409139/83‑1 (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                <  1:1929567/145‑1 (MQ=255)
                                                                                                                                      gggCTGATCCGTATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGt                >  2:1929567/2‑146 (MQ=255)
                                                                                                                                       ggCTGAAACCGCACAGCACACACGATCACCACATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGACGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGAGc           >  2:149847/1‑150 (MQ=255)
                                                                                                                                        gCTGATCCGTATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGa                             <  1:1923569/131‑1 (MQ=255)
                                                                                                                                        gCTGATCCGTATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGa                             >  2:1923569/1‑131 (MQ=255)
                                                                                                                                        gCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCGCCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGa                             >  2:837848/1‑131 (MQ=255)
                                                                                                                                        gCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCGCCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGa                             <  1:837848/131‑1 (MQ=255)
                                                                                                                                          tGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGAGCCCa        <  2:948387/150‑1 (MQ=255)
                                                                                                                                          tGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGAGCCCa        >  1:1979476/1‑150 (MQ=255)
                                                                                                                                           gATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGAGCCCAt       <  1:1974621/150‑1 (MQ=255)
                                                                                                                                            aTACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGAGCCCATg      >  2:579464/1‑150 (MQ=255)
                                                                                                                                             tACGGATCAGCCCACGCGGTCGCCAGAtcatcatca                                                                                                                       >  1:2374913/1‑36 (MQ=255)
                                                                                                                                             tACGGATCAGCCCACGCGGTCGCCAGAtcatcatca                                                                                                                       <  2:2374913/36‑1 (MQ=255)
                                                                                                                                             tACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACATCCTCGGTCGAGCCCATGc     <  2:2020676/150‑1 (MQ=255)
                                                                                                                                                ggATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAAGAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGAGCCCATGCCAc  >  2:1418557/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ATCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCACGGTCGAGCCCATGCCAC  >  NC_002947/5530916‑5531209

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: