Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,531,058 A→C 7.4% 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.4% 191.8 / 5.3 67R421R (CGT→CGGbraEhigh‑affinity branched‑chain amino acid ABC transporter permease BraE
Reads supporting (aligned to +/- strand):  ref base A (31/31);  new base C (2/3);  total (33/34)
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.79e-01

GGCAAATATCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCAC  >  NC_002947/5530909‑5531192
                                                                                                                                                     |                                                                                                                                      
ggCAAATATCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATa                                                                                                                                        >  1:1161584/1‑150 (MQ=255)
ggCAAATATCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCAATCCTTTACGTGGGGTCACACCGGTACGGCTGATa                                                                                                                                        >  1:769592/1‑150 (MQ=255)
   aaaTATCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACgg                                                                                                                                     >  2:884775/1‑150 (MQ=255)
   aaaTATCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACgg                                                                                                                                     >  2:884611/1‑150 (MQ=255)
     atatCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGAt                                                                                                                                   <  2:1408597/150‑1 (MQ=255)
      tatCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATc                                                                                                                                  <  1:884611/150‑1 (MQ=255)
      tatCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATc                                                                                                                                  <  1:884775/150‑1 (MQ=255)
               gcgcCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATCCGTATCAGccgtacc                                                                                                                         <  1:1362946/150‑6 (MQ=255)
                cgcCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCAcgcg                                                                                                                        <  1:1007899/150‑1 (MQ=255)
                 gcCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGTTACGGCTGATACGGATCAGCCCACGCgg                                                                                                                       >  1:771824/1‑150 (MQ=255)
                     gACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGc                                                                                                                   <  2:1161584/150‑1 (MQ=255)
                            aGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTc                                                                                                                     <  2:634016/141‑1 (MQ=255)
                            aGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTc                                                                                                                     >  1:634016/1‑141 (MQ=255)
                              gCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGatcatc                                                                                                          >  1:1332462/1‑150 (MQ=255)
                               cTGACGGCGGTGAGCGCCTTGATGCCACCGAACGTCATCATCAGGTTGTCGACCGAGAGAATGATATCTTCGCGCATGGCGCCACTCCTTTACGTCGTGTCACACCGGTACTGCTGATACGGATCAGCCCACGCGGTCGCCAGatcatca                                                                                                         <  2:534179/150‑1 (MQ=255)
                                    gTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGTCCACGCGGTCGCCAGatc                                                                                                             <  2:238512/141‑1 (MQ=255)
                                    gTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGTCCACGCGGTCGCCAGatc                                                                                                             >  1:238512/1‑141 (MQ=255)
                                    gTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCAcgc                                                                                                                         >  2:554353/1‑129 (MQ=255)
                                    gTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCAcgc                                                                                                                         <  1:554353/129‑1 (MQ=255)
                                             cgcCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATcatcagcaccatcagc                                                                                           <  1:1469604/150‑1 (MQ=255)
                                              gcCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATcatcagcaccatcagca                                                                                          >  2:617902/1‑150 (MQ=255)
                                                cTTGATGCCACCGAACTGCATATTCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGccc                                                                                                                             >  2:1089933/1‑113 (MQ=255)
                                                cTTGATGCCACCGAACTGCATACTCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGccc                                                                                                                             <  1:1089933/113‑1 (MQ=255)
                                                     tGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCCCTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATcatcagcaccatcagc                                                                                           <  1:112177/142‑1 (MQ=255)
                                                     tGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGACACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATcatcagcaccatcagc                                                                                           >  2:112177/1‑142 (MQ=255)
                                                         aCCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGccc                                                                                                                             >  1:256338/1‑104 (MQ=255)
                                                         aCCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGccc                                                                                                                             <  2:256338/104‑1 (MQ=255)
                                                           cGAACTGCATCATCAGGTTGTCGACCGAGAAAGCGATGTCGTCGCTGGTTGCGGCACTCCTTTACGTGGGCTCAAACAGGTCGTACTGAGACAGTGGATCCCACGCGGTCGCCAGATCATCATCAGCACCATCACCACGCCGAAcagcag                                                                             >  1:372851/1‑150 (MQ=255)
                                                             aaCTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAgcagca                                                                           <  2:1332462/150‑1 (MQ=255)
                                                               cTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCAcc                                                                         >  1:353022/1‑150 (MQ=255)
                                                                  catcatCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCgg                                                                      >  1:872511/1‑150 (MQ=255)
                                                                       tcaGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTc                                                                 <  2:771824/150‑1 (MQ=255)
                                                                          ggTTTTCGACCCAGATAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGga                                                              <  2:769592/150‑1 (MQ=255)
                                                                           gTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGc                                                                                       >  2:1479990/1‑124 (MQ=255)
                                                                           gTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGc                                                                                       <  1:1479990/124‑1 (MQ=255)
                                                                                  aCCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGaa                                                                                   >  2:152952/1‑121 (MQ=255)
                                                                                  aCCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGaa                                                                                   <  1:152952/121‑1 (MQ=255)
                                                                                          aaTGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATCCgt                                                                                                                                     >  1:1566337/1‑62 (MQ=255)
                                                                                          aaTGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATCCgt                                                                                                                                     <  2:1566337/63‑2 (MQ=255)
                                                                                               tatCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCgg                                         >  1:230273/1‑150 (MQ=255)
                                                                                                          aTGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCa                              <  2:372851/150‑1 (MQ=255)
                                                                                                              cgcCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCAc                          <  1:1421964/150‑1 (MQ=255)
                                                                                                                 cccTCCTTTACGTGGGGTCCCCACGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGaa                       <  2:353022/148‑1 (MQ=255)
                                                                                                                  aCTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCAc                          >  2:1066412/1‑146 (MQ=255)
                                                                                                                  aCTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCAc                          <  1:1066412/146‑1 (MQ=255)
                                                                                                                    tCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAAcgc                    >  1:1034159/1‑150 (MQ=255)
                                                                                                                      cTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCgcagc                                                 >  1:625432/1‑119 (MQ=255)
                                                                                                                      cTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCgcagc                                                 <  2:625432/119‑1 (MQ=255)
                                                                                                                              ggggTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGAt              >  1:1555556/1‑146 (MQ=255)
                                                                                                                              ggggTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGAt              <  2:1555556/146‑1 (MQ=255)
                                                                                                                                 gTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACGATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACcaca       >  2:1345247/1‑150 (MQ=255)
                                                                                                                                    acacCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAg                                                          <  1:605416/96‑1 (MQ=255)
                                                                                                                                    acacCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAg                                                          >  2:605416/1‑96 (MQ=255)
                                                                                                                                    acacCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACAccc    >  1:126865/1‑150 (MQ=255)
                                                                                                                                    acacCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACAccc    >  1:119078/1‑150 (MQ=255)
                                                                                                                                      acCGGTACGGCTGATACGGATCAGTCCACGCGGTCGCCAGATCATcatcagcaccat                                                                                               >  2:213194/1‑57 (MQ=255)
                                                                                                                                      acCGGTACGGCTGATACGGATCAGTCCACGCGGTCGCCAGATCATcatcagcaccat                                                                                               <  1:213194/57‑1 (MQ=255)
                                                                                                                                      acCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAGCGCCGCGATGACCACACCCAc  >  1:1536843/1‑150 (MQ=255)
                                                                                                                                      acCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCAc  <  2:309435/150‑1 (MQ=255)
                                                                                                                                      acCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCAc  <  2:1034159/150‑1 (MQ=255)
                                                                                                                                      acCGGTACAGCTGATACGGATCAGCCCACGCGGACGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCAGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCAc  >  2:321592/1‑150 (MQ=255)
                                                                                                                                          gTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGAc           >  2:812856/1‑137 (MQ=255)
                                                                                                                                          gTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGAc           <  1:812856/137‑1 (MQ=255)
                                                                                                                                           tACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGa            >  2:452887/1‑135 (MQ=255)
                                                                                                                                           tACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGa            <  1:452887/135‑1 (MQ=255)
                                                                                                                                               gCTGATCCGTATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCgca                                                   <  2:459254/92‑1 (MQ=255)
                                                                                                                                               gCTGATCCGTATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCgca                                                   >  1:459254/1‑92 (MQ=255)
                                                                                                                                                 tGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAAca                                                                                 <  2:376100/60‑1 (MQ=255)
                                                                                                                                                 tGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAAca                                                                                 >  1:376100/1‑60 (MQ=255)
                                                                                                                                                     |                                                                                                                                      
GGCAAATATCTGGTTGCGCCGGACCTTCAGGCTGACGTCGCTGAGCGCCTTGATGCCACCGAACTGCATCATCAGGTTGTCGACCGAGAGAATGATATCGTCGCTCATGGCGCCACTCCTTTACGTGGGGTCACACCGGTACGGCTGATACGGATCAGCCCACGCGGTCGCCAGATCATCATCAGCACCATCAGCACGCCGAACAGCAGCACCCGGTACTCGGAGAAGCTGCGCAGCAGTTCCGGCGCCACGGTCAGCACGAACGCCGCGATGACCACACCCAC  >  NC_002947/5530909‑5531192

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