Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,797,388 G→T 16.4% P186Q (CCG→CAG)  pheT ← phenylalanine tRNA synthetase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,797,3880GT16.4% 240.2 / 29.9 94P186Q (CCG→CAG) pheTphenylalanine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (41/37);  new base T (3/13);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 1.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAA  >  NC_000913/1797289‑1797487
                                                                                                   |                                                                                                   
cGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCg                                                                                                     <  2:2132227/100‑1 (MQ=255)
  gCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGc                                                                                                   <  2:24740/100‑1 (MQ=255)
   ccAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCa                                                                                                  <  1:1245723/100‑1 (MQ=255)
     aaGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGc                                                                                                >  1:2495911/1‑100 (MQ=255)
     aaGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGc                                                                                                <  2:771272/100‑1 (MQ=255)
      aGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCt                                                                                               >  2:1081348/1‑100 (MQ=255)
      aGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCt                                                                                               >  2:1496219/1‑100 (MQ=255)
         tAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGt                                                                                            <  2:1576336/100‑1 (MQ=255)
         tAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGt                                                                                            >  2:798829/1‑100 (MQ=255)
         tAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGt                                                                                            <  2:1300664/100‑1 (MQ=255)
          aaCGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGt                                                                                            >  1:1625632/1‑99 (MQ=255)
          aaCGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGCTGAACCAGCTGCCGCTGGtt                                                                                           <  2:54046/100‑1 (MQ=255)
           aCGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTc                                                                                          >  1:1348268/1‑100 (MQ=255)
           aCGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTc                                                                                          >  2:1881802/1‑100 (MQ=255)
            cgcgGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGCACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCa                                                                                         >  1:2498037/1‑100 (MQ=255)
            cgcgGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCa                                                                                         >  1:1213235/1‑100 (MQ=255)
            cgcgGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCa                                                                                         >  1:1194137/1‑100 (MQ=255)
            cgcgGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCa                                                                                         >  1:2115422/1‑100 (MQ=255)
              cgGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGATCGATTTCCGGTTGAACCAGCTGCCGCTGGtt                                                                                           <  1:2577870/96‑1 (MQ=255)
              cgGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaac                                                                                       <  2:1469451/100‑3 (MQ=255)
              cgGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaac                                                                                       <  2:1653609/100‑3 (MQ=255)
               gggCCGGCTTCCGGGGCTTCGGCTGTTATCGGCGGCGTGTCGTTGATGGTCGCGCCACCCGGAACGATTTCCGGTTGAACCCGCTGCCGCTGGTTCaac                                                                                       <  1:1055640/99‑3 (MQ=255)
               gggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGACCGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCa                                                                                      <  2:1822836/100‑1 (MQ=255)
               gggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaacc                                                                                      <  2:2317181/100‑4 (MQ=255)
               gggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaacc                                                                                      <  2:1386981/100‑4 (MQ=255)
               gggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaacc                                                                                      <  2:314789/100‑4 (MQ=255)
               gggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGCTGAACCAGCTGCCGCTGGTTCaacc                                                                                      <  2:1235609/100‑4 (MQ=38)
               gggCAGGCTTACGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaacc                                                                                      <  2:84300/100‑4 (MQ=255)
               gggCAGGCCTCCGGCGCTGCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaacc                                                                                      <  1:727812/100‑4 (MQ=255)
                ggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaacc                                                                                      <  1:722184/99‑4 (MQ=255)
                  cAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACgg                                                                                   >  2:1002837/1‑100 (MQ=255)
                  cAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACgg                                                                                   >  2:864421/1‑100 (MQ=255)
                     gCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCaa                                                                                >  1:2008836/1‑100 (MQ=255)
                        ttccgcgcTTCGACTGTAGTCTGCAGCGTGTCGTCGATGGTCGCACGAACCGGAACGATGTCCGTTTGAACCAGCTGCCGCTGGTTCaacc                                                                                      <  1:57318/87‑4 (MQ=25)
                         ccGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAGCGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACg                                                                              >  2:520464/1‑98 (MQ=255)
                           ggCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCAc                                                                          <  2:1594986/100‑1 (MQ=255)
                                ttCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCa                                                                     >  2:1279764/1‑100 (MQ=255)
                                ttCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCa                                                                     <  2:2008847/100‑1 (MQ=255)
                                    aCTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAacac                                                                 <  2:892570/100‑1 (MQ=255)
                                     cTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACAcc                                                                >  2:150322/1‑100 (MQ=255)
                                      tGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAGCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCATGGCAACGTCACGCGCAacacn                                                                <  1:166723/99‑2 (MQ=255)
                                       gcaaTCGGCAGCGTGTCGTCGATCGTCGCACCAAGCGGAACGATTGCCGGGTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAactc                                                                 <  1:1037503/95‑3 (MQ=255)
                                        tAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                             >  2:94980/1‑100 (MQ=255)
                                        tAATCGGCAGCGTGTCGTCGATGGTCGCACCAACAGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCaa                                                              >  1:909327/1‑99 (MQ=255)
                                           tCGGCAGCGTGTCGTCGATGGGCGCACCAAGCGGAAGGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatga                                                           <  1:1749842/99‑1 (MQ=255)
                                            cGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatgatg                                                         <  1:1071967/100‑1 (MQ=255)
                                            cGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatgatg                                                         >  2:942428/1‑100 (MQ=255)
                                             ggCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGc                                                        >  1:868782/1‑100 (MQ=255)
                                                 gCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTa                                                     <  1:337584/99‑1 (MQ=255)
                                                   gtgtCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGc                                                  >  1:1066716/1‑100 (MQ=255)
                                                   gtgtCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGc                                                  <  2:1569512/100‑1 (MQ=255)
                                                   gtgtCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGc                                                  <  2:1499107/100‑1 (MQ=255)
                                                     gtcgtcGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAg                                                <  2:650443/100‑1 (MQ=255)
                                                         tcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCgg                                            <  1:673390/100‑1 (MQ=255)
                                                          cgATGGTCGCAGCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGc                                           >  1:1938479/1‑100 (MQ=255)
                                                          cgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGc                                           <  1:1002570/100‑1 (MQ=255)
                                                          cgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGc                                           >  2:760466/1‑100 (MQ=255)
                                                                    aCCAACCGGAACGATTTCCGGTGGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTGACGCGCAACACGAATGATGCCTAAGCAGTCCGCACGGttt                                    >  1:637124/1‑97 (MQ=255)
                                                                     ccAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCg                                >  1:2332838/1‑100 (MQ=255)
                                                                      cAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGt                               >  1:2174851/1‑100 (MQ=255)
                                                                       aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTc                              <  2:785533/100‑1 (MQ=255)
                                                                       aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGGCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCg                                <  1:850743/98‑1 (MQ=255)
                                                                         ccGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCAc                            >  2:312597/1‑100 (MQ=255)
                                                                         ccGGAACGATTTCCGGTTGAACCAGCGGCAGCGGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCAc                            <  2:1838323/100‑1 (MQ=255)
                                                                           ggaaatcgttCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGc                          >  1:1415915/9‑100 (MQ=255)
                                                                           ggAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCCACACCCATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGc                          >  1:2421162/1‑100 (MQ=255)
                                                                           ggAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGc                          >  2:2005376/1‑100 (MQ=255)
                                                                           ggAACGATTCCCGTTGGAACCGTCGGCAGCTGGTTCAGCACGGCAACGACACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACgg                          <  1:2450239/100‑2 (MQ=255)
                                                                           gcaaCGATTTCCGGTTGAACCCGCGGCAGCTGCTTCAGCACGGCAACGTCACGCGCAACCCCAATGATGCCTAAGCAGTCGGCACGGTTTGGc                                 <  1:680642/91‑1 (MQ=255)
                                                                             aaCGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACATCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTg                        <  1:1410301/100‑1 (MQ=255)
                                                                             aaCGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTg                        <  2:868785/100‑1 (MQ=255)
                                                                                 aTTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGAttt                    >  2:890147/1‑100 (MQ=255)
                                                                                  tttCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTc                   <  1:291925/100‑1 (MQ=255)
                                                                                  tttCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTc                   >  2:1426653/1‑100 (MQ=255)
                                                                                  tttCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAACCATGCCTAAGCAGTCGGCACGGTTTGGCGCCACGCTGATTTc                   <  1:292950/100‑1 (MQ=255)
                                                                                     ccGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGAt                <  2:1760536/100‑1 (MQ=255)
                                                                                      cGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATg               <  2:2059897/100‑1 (MQ=255)
                                                                                       ggTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATgg              >  1:667767/1‑100 (MQ=255)
                                                                                       ggTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATgg              <  2:671130/100‑1 (MQ=255)
                                                                                        gTTGAACCAGCGGCGGCTGGGTCAGCACGACAGCGTCACGCGCAACCCCAAGGATGGCTAAGCAGGCGGCACGGTTTGGCGTCACGCTGATTTAGAt                <  1:789034/97‑1 (MQ=255)
                                                                                         ttGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgtg            <  2:1086250/100‑1 (MQ=255)
                                                                                         ctGAACCAGCTGCCGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgtg            >  2:1083808/2‑100 (MQ=255)
                                                                                         ccgAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgtg            >  2:1408568/3‑100 (MQ=255)
                                                                                          tGAACCAGCTGCCGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgtgt           >  2:2386618/1‑100 (MQ=255)
                                                                                          tGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgtgt           <  1:696504/100‑1 (MQ=255)
                                                                                           gAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgtt          <  2:2155124/100‑1 (MQ=255)
                                                                                            aaCCAGCTGCCGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgttg         >  2:984169/1‑100 (MQ=255)
                                                                                            aaCCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgttg         >  2:658823/1‑100 (MQ=255)
                                                                                             aCCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgttgt        >  2:2282757/1‑100 (MQ=255)
                                                                                             aCCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgttgt        >  1:77119/1‑100 (MQ=255)
                                                                                             aCCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgttgt        >  1:1190167/1‑100 (MQ=255)
                                                                                               cAGCGGTAGCTGGTTCAGCACGGCAACGTTACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGtca      <  1:740127/100‑1 (MQ=255)
                                                                                               cAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGtca      <  1:1881490/100‑1 (MQ=255)
                                                                                                aGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGtcat     >  2:542851/1‑100 (MQ=255)
                                                                                                aGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGtcat     >  1:1215577/1‑100 (MQ=255)
                                                                                                   ggCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCaa  >  1:221327/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAA  >  NC_000913/1797289‑1797487

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