Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,217,147 (G)5→6 coding (261/927 nt) yehX ← glycine betaine ABC transporter ATP binding subunit YehX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,217,1421.G99.0% 809.5 / 0.2 191coding (266/927 nt)yehXglycine betaine ABC transporter ATP binding subunit YehX
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (99/90);  total (99/92)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 231 total aligned reads displayed.

CCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATC  >  NC_000913/2217042‑2217237
                                                                                                     |                                                                                               
ccccAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGt                                                                                                  >  2:3226758/1‑101 (MQ=255)
ccccAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGt                                                                                                  >  1:3226757/1‑101 (MQ=255)
 cccAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTg                                                                                                <  1:3226228/101‑1 (MQ=255)
 cccAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTg                                                                                                >  1:3226760/1‑101 (MQ=255)
 cccAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTg                                                                                                <  1:3226417/101‑1 (MQ=255)
   cAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTggg                                                                                              >  1:3226763/1‑101 (MQ=255)
    aGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTgggg                                                                                             <  2:3226446/101‑1 (MQ=255)
    aGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTgggg                                                                                             >  1:3226765/1‑101 (MQ=255)
     gTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTggggg                                                                                            >  1:3226768/1‑101 (MQ=255)
     gTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTggggg                                                                                            <  1:3226353/101‑1 (MQ=255)
      tAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTgggggg                                                                                           <  2:3226206/101‑2 (MQ=255)
      tAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTgggggg                                                                                           <  2:3226154/101‑2 (MQ=255)
      tAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTgggggg                                                                                           >  2:3226770/1‑100 (MQ=255)
       aGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTgggggga                                                                                          <  1:3226471/101‑3 (MQ=255)
       aGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTgggggga                                                                                          >  1:3226772/1‑99 (MQ=255)
         cgcCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTggggggaac                                                                                        >  1:3226775/1‑97 (MQ=255)
          gcCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTggggggaaca                                                                                       <  2:3226451/101‑6 (MQ=255)
          gcCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTggggggaaca                                                                                       >  2:3226777/1‑96 (MQ=255)
           ccATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAg                                                                                       <  1:3226283/101‑1 (MQ=255)
           ccATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAg                                                                                       >  1:3226781/1‑101 (MQ=255)
           ccATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAg                                                                                       >  2:3226782/1‑101 (MQ=255)
            cATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAgg                                                                                      <  2:3226274/101‑1 (MQ=255)
             aTTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGc                                                                                     <  2:3226466/101‑1 (MQ=255)
              ttAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGcc                                                                                    >  1:3226786/1‑101 (MQ=255)
               tAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCa                                                                                   >  2:3226787/1‑101 (MQ=255)
                 aTTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAAt                                                                                 <  1:3226344/101‑1 (MQ=255)
                    cgtcgATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGa                                                                              <  1:3226561/101‑1 (MQ=255)
                     gtcgATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGAt                                                                             <  2:3226576/101‑1 (MQ=255)
                      tcgATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGAtt                                                                            <  2:3226304/101‑1 (MQ=255)
                       cgATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTg                                                                           <  2:3226505/101‑1 (MQ=255)
                        gATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGa                                                                          >  2:3226791/1‑101 (MQ=255)
                        gATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGa                                                                          >  1:3226792/1‑101 (MQ=255)
                           aCGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATg                                                                       <  1:3226325/101‑1 (MQ=255)
                           aCGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATg                                                                       >  1:3226794/1‑101 (MQ=255)
                               tCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCAt                                                                   >  1:3226799/1‑101 (MQ=255)
                               tCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCAt                                                                   <  2:3226610/101‑1 (MQ=255)
                               tCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCAt                                                                   >  2:3226798/1‑101 (MQ=255)
                                cGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATa                                                                  <  1:3226425/101‑1 (MQ=255)
                                 gTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAg                                                                 >  2:3226801/1‑101 (MQ=255)
                                 gTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAg                                                                 <  1:3226578/101‑1 (MQ=255)
                                  tCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGc                                                                <  2:3226430/101‑1 (MQ=255)
                                    aaTCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGccc                                                              >  2:3226805/1‑101 (MQ=255)
                                     aTCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCa                                                             >  1:3226806/1‑101 (MQ=255)
                                      tCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCAt                                                            <  1:3226552/101‑1 (MQ=255)
                                       ccGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATc                                                           <  2:3226183/101‑1 (MQ=255)
                                        cgcgCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAAt                                                                        >  2:3226808/1‑87 (MQ=255)
                                        cgcgCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATcc                                                          <  1:3226397/101‑1 (MQ=255)
                                          cgcCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCcgg                                                        >  2:3226810/1‑101 (MQ=255)
                                           gcCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCcggc                                                       >  2:3226812/1‑101 (MQ=255)
                                            cccGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCcggcg                                                      <  2:3226587/101‑1 (MQ=255)
                                            cccGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCcggcg                                                      >  1:3226814/1‑101 (MQ=255)
                                             ccGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCcggcgg                                                     <  1:3226486/101‑1 (MQ=255)
                                             ccGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCcggcgg                                                     <  1:3226104/101‑1 (MQ=255)
                                              cgcgACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCggcggc                                                    <  1:3226385/101‑1 (MQ=255)
                                                cgACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCAt                                                            >  2:3226819/1‑91 (MQ=255)
                                                 gACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCa                                                 >  1:3226821/1‑101 (MQ=255)
                                                  aCCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCaa                                                >  1:3226823/1‑101 (MQ=255)
                                                   ccATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAAc                                               >  2:3226824/1‑101 (MQ=255)
                                                    cATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACt                                              <  2:3226535/101‑1 (MQ=255)
                                                    cATTTGTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACt                                              >  1:3226826/1‑101 (MQ=255)
                                                     aTTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTc                                             <  1:3226522/101‑1 (MQ=255)
                                                     aTTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTc                                             <  2:3226427/101‑1 (MQ=255)
                                                      tttttGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTcc                                            <  1:3226593/101‑1 (MQ=255)
                                                       ttttGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCa                                           >  2:3226831/1‑101 (MQ=255)
                                                       ttttGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCa                                           >  1:3226830/1‑101 (MQ=255)
                                                        tttGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAg                                          >  1:3226834/1‑101 (MQ=255)
                                                         ttGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGt                                         >  1:3226837/1‑101 (MQ=255)
                                                           gTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTAc                                       <  1:3226434/101‑1 (MQ=255)
                                                            tAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACt                                      <  2:3226340/101‑1 (MQ=255)
                                                            tAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACt                                      <  1:3226612/101‑1 (MQ=255)
                                                            tAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACt                                      >  1:3226839/1‑101 (MQ=255)
                                                             aaTTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTg                                     <  1:3226464/101‑1 (MQ=255)
                                                             aaTTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTg                                     >  2:3226841/1‑101 (MQ=255)
                                                               ttGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGc                                   <  2:3226418/101‑1 (MQ=255)
                                                                tGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCa                                  >  2:3226843/1‑101 (MQ=255)
                                                                tGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCa                                  <  2:3226424/101‑1 (MQ=255)
                                                                     cACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCagcgag                             >  1:3226847/1‑101 (MQ=255)
                                                                     cACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCagcgag                             >  2:3226846/1‑101 (MQ=255)
                                                                     cACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCagcgag                             <  2:3226356/101‑1 (MQ=255)
                                                                       cGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCagcgagcg                           >  1:3226850/1‑101 (MQ=255)
                                                                        ggTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAgcgagcga                          <  2:3226401/101‑1 (MQ=255)
                                                                          tAGCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAAt                        >  2:3226853/1‑101 (MQ=255)
                                                                            gCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAAttt                      >  1:3226855/1‑101 (MQ=255)
                                                                            gCGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAAttt                      >  2:3226856/1‑101 (MQ=255)
                                                                             cGATGTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATttc                     >  2:3226858/1‑101 (MQ=255)
                                                                                 gTTTTGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTtcttct                 >  1:3226863/1‑101 (MQ=255)
                                                                                   tttGCGCCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTcc               <  1:3226367/101‑1 (MQ=255)
                                                                                       cgcCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCa           <  1:3226310/101‑1 (MQ=255)
                                                                                        gcCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCaa          >  2:3226865/1‑101 (MQ=255)
                                                                                        gcCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCaa          >  1:3226866/1‑101 (MQ=255)
                                                                                        gcCACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCaa          <  2:3226378/101‑1 (MQ=255)
                                                                                         ccACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATgg                                                                      >  1:3226870/1‑40 (MQ=255)
                                                                                         ccACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCagcgag                             >  2:3226869/1‑81 (MQ=255)
                                                                                         ccACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCaaa         <  1:3226500/101‑1 (MQ=255)
                                                                                          cACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAg        <  1:3226489/101‑1 (MQ=255)
                                                                                          cACGCTCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAg        >  2:3226872/1‑101 (MQ=255)
                                                                                               tCCAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGAt   <  2:3226408/101‑1 (MQ=255)
                                                                                                ccAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGc       <  1:3226659/96‑1 (MQ=255)
                                                                                                ccAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGAtc  >  1:3226879/1‑101 (MQ=255)
                                                                                                ccAGTGGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGAtc  >  1:3226877/1‑101 (MQ=255)
                                                                                                     |                                                                                               
CCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGCGAGCGAATTTCTTCTCCGGCAAAGCGGATC  >  NC_000913/2217042‑2217237

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