Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 531,655 C→A A189A (GCG→GCT allS ← allD operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913531,6550CA97.7% 298.0 / ‑0.5 86A189A (GCG→GCTallSallD operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (36/48);  total (36/50)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

AATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTA  >  NC_000913/531559‑531751
                                                                                                |                                                                                                
aataatCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCt                                                                                               >  2:1876767/1‑100 (MQ=255)
aataatCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCt                                                                                               >  1:1668256/1‑100 (MQ=255)
    atCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCTTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTg                                                                                           <  2:294749/100‑1 (MQ=255)
    atCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTg                                                                                           <  2:420215/100‑1 (MQ=255)
      ctctTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGtt                                                                                         >  2:1624992/1‑100 (MQ=255)
      ctctTGTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCCTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTgct                                                                                         <  1:1580635/100‑3 (MQ=255)
       tctTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTa                                                                                        >  1:931524/1‑100 (MQ=255)
       tctTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTa                                                                                        >  1:492844/1‑100 (MQ=255)
       tctTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTa                                                                                        >  1:2158568/1‑100 (MQ=255)
       tctTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTa                                                                                        >  1:2105631/1‑100 (MQ=255)
        ctTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAg                                                                                       >  2:2034181/1‑100 (MQ=255)
         tttttGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGc                                                                                      <  2:928536/100‑1 (MQ=255)
          ttttGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCg                                                                                     >  1:2100862/1‑100 (MQ=255)
          ttttGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCg                                                                                     >  1:101905/1‑100 (MQ=255)
          ttttGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCg                                                                                     <  2:1499887/100‑1 (MQ=255)
          ttttGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCg                                                                                     <  1:1752073/100‑1 (MQ=255)
           tttGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCgg                                                                                    <  2:2105624/100‑1 (MQ=255)
           tttGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCgg                                                                                    <  2:931522/100‑1 (MQ=255)
            ttGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGc                                                                                   <  1:213048/100‑1 (MQ=255)
              gCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGctc                                                                                 >  1:919113/1‑100 (MQ=255)
              gCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGctc                                                                                 >  1:29545/1‑100 (MQ=255)
                ccTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCtt                                                                               <  2:2264434/100‑1 (MQ=255)
                ccTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCtt                                                                               >  1:883746/1‑100 (MQ=255)
                 cTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGAGCGTCGGAAATCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTc                                                                              >  2:1641918/1‑100 (MQ=255)
                 cTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGAGCGACGGAAATCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTc                                                                              >  2:1641911/1‑100 (MQ=255)
                 cTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTc                                                                              <  2:411744/100‑1 (MQ=255)
                  tGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCa                                                                             >  1:1947967/1‑100 (MQ=255)
                    gCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAAc                                                                           >  2:678534/1‑100 (MQ=255)
                     cgattgCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACg                                                                          >  1:446084/6‑100 (MQ=255)
                     cAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACg                                                                          >  1:1254304/1‑100 (MQ=255)
                       aTCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGtt                                                                        <  2:466544/100‑1 (MQ=255)
                        tCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTc                                                                       <  1:1546003/100‑1 (MQ=255)
                              ggCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGc                                                                 >  1:260251/1‑100 (MQ=255)
                               gCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCg                                                                <  1:1642636/100‑1 (MQ=255)
                                cGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCgg                                                               <  2:477784/100‑1 (MQ=255)
                                cGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCgg                                                               <  2:1937535/100‑1 (MQ=255)
                                cGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCgg                                                               <  2:2058479/100‑1 (MQ=255)
                                 gACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGga                                                              <  2:1158515/100‑1 (MQ=255)
                                  aCGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGgat                                                             >  1:1866003/1‑100 (MQ=255)
                                   tctgTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGgatg                                                            <  2:942414/97‑1 (MQ=255)
                                   cgcgTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGgatg                                                            <  2:101903/100‑1 (MQ=255)
                                     cgTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGgatgat                                                          <  2:1172019/100‑1 (MQ=255)
                                      gTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGgatgaga                                                         >  1:942416/1‑98 (MQ=255)
                                      gTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATc                                                         <  2:559727/100‑1 (MQ=255)
                                      gTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATc                                                         <  1:861038/100‑1 (MQ=255)
                                          ggTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATccgcc                                                     >  1:2278983/1‑100 (MQ=255)
                                           gTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATccgccg                                                    <  2:1947960/100‑1 (MQ=255)
                                           gTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATccgccg                                                    <  2:373123/100‑1 (MQ=255)
                                            ttAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCcgccgc                                                   >  2:1139661/1‑100 (MQ=255)
                                             tAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCgccgcc                                                  >  2:596330/1‑100 (MQ=255)
                                             tAAGGTGCGGGCGCAGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCgccgcc                                                  <  2:1591760/100‑1 (MQ=255)
                                              aaGGTGCGGGCGCTGTCTTCAATATTGACCGCCGTAAAGCGCCGCAACTGAGCTTCTGTGAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCgccgcc                                                  <  1:302800/99‑1 (MQ=255)
                                              aaGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCa                                                 >  2:1618073/1‑100 (MQ=255)
                                              aaGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCa                                                 <  2:501312/100‑1 (MQ=255)
                                              aaGGTGCGGGCGCTGTCTTCAAGATTGACCGACGGAACGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCa                                                 <  2:1576746/100‑1 (MQ=255)
                                               aGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCGGCGGATCATccgctggcga                                                <  1:1663153/100‑7 (MQ=255)
                                                ggTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATg                                               >  1:994428/1‑100 (MQ=255)
                                                ggTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATg                                               <  1:421077/100‑1 (MQ=255)
                                                  tGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGAc                                             <  1:1968873/100‑1 (MQ=255)
                                                          cTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGcc                                     >  2:510338/1‑100 (MQ=255)
                                                           tGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCa                                    <  1:1641916/100‑1 (MQ=255)
                                                           tGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTCAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCa                                    <  1:1641923/100‑1 (MQ=255)
                                                            gTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCAt                                   >  1:304025/1‑100 (MQ=255)
                                                             tCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCAtt                                  <  1:2078590/100‑1 (MQ=255)
                                                              cTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTg                                 >  2:1541180/1‑100 (MQ=255)
                                                              cTTCAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTg                                 >  1:1499884/1‑100 (MQ=255)
                                                                 cAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCAc                              <  2:883744/100‑1 (MQ=255)
                                                                 cAATATTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCAc                              <  1:1733897/100‑1 (MQ=255)
                                                                   atatTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCAcc                             <  2:1984885/99‑1 (MQ=255)
                                                                   atatTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCg                            >  2:2137295/1‑100 (MQ=255)
                                                                     atTGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGAt                          >  2:602313/1‑100 (MQ=255)
                                                                       tGACCGCCGGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATcc                        <  2:2100855/100‑1 (MQ=255)
                                                                             ccgGAAAGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAggg                  <  2:994426/100‑1 (MQ=255)
                                                                                 aaaGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCa              <  2:29544/100‑1 (MQ=255)
                                                                                 aaaGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCa              <  1:2034188/100‑1 (MQ=255)
                                                                                 aaaGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCa              <  2:919111/100‑1 (MQ=255)
                                                                                  aaGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCaa             <  1:1876773/100‑1 (MQ=255)
                                                                                   aGCGCCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAg            <  2:260249/100‑1 (MQ=255)
                                                                                    gcgcCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGa           >  1:505501/1‑100 (MQ=255)
                                                                                    gcgcCGCAACTGAGCTTCTGTTAGCGGCTCTTCAACGGTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGa           <  1:86381/100‑1 (MQ=255)
                                                                                      gccgcAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGAc          <  2:1865997/99‑1 (MQ=255)
                                                                                       ccgcAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTa        >  2:1098882/1‑100 (MQ=255)
                                                                                          cAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAg     <  1:1893331/100‑1 (MQ=255)
                                                                                          cAACTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAg     <  2:2158561/100‑1 (MQ=255)
                                                                                           aaCTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAgg    <  2:492842/100‑1 (MQ=255)
                                                                                             cTGAGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTa  >  1:320967/1‑100 (MQ=255)
                                                                                                |                                                                                                
AATAATCTCTTTTTGCCCTGGCAATCGCCAGGCGACGCGTTTGGTTAAGGTGCGGGCGCTGTCTTCAATATTGACCGCCGGAAAGCGCCGCAACTGCGCTTCTGTTAGCGGCTCTTCAACGTTCGCCAGCGGATGATCCGCCGCCATGACAAAGCGCCATTGCACCGATCCTAAGGGATCAAGACTAAAGGTA  >  NC_000913/531559‑531751

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