Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,797,391 A→C 16.9% L185R (CTG→CGG)  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,3910AC16.9% 213.4 / 39.0 95L185R (CTG→CGG) pheTphenylalanine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (42/36);  new base C (3/13);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 1.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTT  >  NC_000913/1797292‑1797490
                                                                                                   |                                                                                                   
ccAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCa                                                                                                     <  1:1245723/100‑1 (MQ=255)
  aaGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGc                                                                                                   <  2:771272/100‑1 (MQ=255)
  aaGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGc                                                                                                   >  1:2495911/1‑100 (MQ=255)
   aGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCt                                                                                                  >  2:1496219/1‑100 (MQ=255)
   aGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCt                                                                                                  >  2:1081348/1‑100 (MQ=255)
      tAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGt                                                                                               >  2:798829/1‑100 (MQ=255)
      tAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGt                                                                                               <  2:1576336/100‑1 (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:2115422/1‑100 (MQ=255)
         cgcgGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCa                                                                                            >  1:1194137/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:314789/100‑4 (MQ=255)
            gggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaacc                                                                                         <  2:1386981/100‑4 (MQ=255)
            gggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaacc                                                                                         <  2:2317181/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                                              >  2:760466/1‑100 (MQ=255)
                                                       cgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGc                                              <  1:1002570/100‑1 (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                      >  2:1426653/1‑100 (MQ=255)
                                                                               tttCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTc                      <  1:291925/100‑1 (MQ=255)
                                                                               tttCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAACCATGCCTAAGCAGTCGGCACGGTTTGGCGCCACGCTGATTTc                      <  1:292950/100‑1 (MQ=255)
                                                                                  ccGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGAt                   <  2:1760536/100‑1 (MQ=255)
                                                                                   cGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATg                  <  2:2059897/100‑1 (MQ=255)
                                                                                    ggTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATgg                 <  2:671130/100‑1 (MQ=255)
                                                                                    ggTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATgg                 >  1:667767/1‑100 (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           >  1:77119/1‑100 (MQ=255)
                                                                                          aCCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgttgt           >  2:2282757/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        >  1:1215577/1‑100 (MQ=255)
                                                                                             aGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGtcat        >  2:542851/1‑100 (MQ=255)
                                                                                                ggCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCaa     >  1:221327/1‑100 (MQ=255)
                                                                                                 gCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAg    <  2:1190171/100‑1 (MQ=255)
                                                                                                  cAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGt   <  2:69350/100‑1 (MQ=255)
                                                                                                   aGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGtt  <  2:1774242/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTT  >  NC_000913/1797292‑1797490

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: