Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,428,391 G→C 15.7% L112L (CTC→CTG insH1 ← IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,428,3910GC15.7% 250.1 / 36.5 89L112L (CTC→CTGinsH1IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base G (38/37);  new base C (9/5);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

TGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCA  >  NC_000913/1428294‑1428490
                                                                                                 |                                                                                                   
tggCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAg                                                                                                   <  2:214452/100‑1 (MQ=255)
 ggCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACCATTGACGGGCCAGTTGATGCTGCTCGAGc                                                                                                  >  1:1547519/1‑100 (MQ=255)
   cATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAg                                                                                                <  1:394948/100‑1 (MQ=255)
   cATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAg                                                                                                <  2:1990708/100‑1 (MQ=255)
    aTCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAgg                                                                                               <  1:496556/100‑1 (MQ=255)
    aTCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAgg                                                                                               <  2:493524/100‑1 (MQ=255)
     tCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGt                                                                                              <  1:906771/100‑1 (MQ=255)
     tCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGt                                                                                              <  2:1285506/100‑1 (MQ=255)
       caccaAAGTGCCTTGGGTCATCATGACGCCTGCTTCGTTCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTg                                                                                             >  1:1791162/1‑99 (MQ=255)
       caccaAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTgg                                                                                            >  1:1484683/1‑100 (MQ=255)
        accaAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGc                                                                                           >  2:447634/1‑100 (MQ=255)
          caAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCgg                                                                                         <  2:2506021/100‑1 (MQ=255)
              gTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAAt                                                                                     <  1:1064966/100‑1 (MQ=255)
              gTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAAt                                                                                     <  1:339757/100‑1 (MQ=40)
              gTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAAt                                                                                     <  1:2375343/100‑1 (MQ=40)
              gTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAAt                                                                                     <  2:450539/100‑1 (MQ=40)
              gTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAAt                                                                                     <  2:450510/100‑1 (MQ=255)
                gCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTc                                                                                   <  2:1585921/100‑1 (MQ=40)
                 ccTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCa                                                                                  <  1:1720600/100‑1 (MQ=255)
                    tGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCatga                                                                               >  2:612663/1‑100 (MQ=255)
                     gggTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCatgat                                                                              <  2:2782354/100‑1 (MQ=40)
                          atcatGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGAtggtgg                                                                         <  2:2715945/100‑1 (MQ=38)
                          atcatGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGAtggtgg                                                                         >  2:2739668/1‑100 (MQ=38)
                           tcatGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATggtggt                                                                        >  1:2096871/1‑100 (MQ=38)
                           tcatGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATggtggt                                                                        >  2:1639728/1‑100 (MQ=38)
                             atGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGc                                                                      <  2:462073/100‑1 (MQ=38)
                                aCGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGa                                                                    >  2:151568/1‑99 (MQ=40)
                                             ccagccagCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCt                                                      >  1:428104/1‑100 (MQ=255)
                                                gccagcGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGCTGGTGGTGCGATCCGGCAGGGCGCTATc                                                   >  1:201627/1‑100 (MQ=255)
                                                  cagcGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCa                                                 >  2:2295188/1‑100 (MQ=255)
                                                   agcGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAg                                                >  2:2424616/1‑100 (MQ=255)
                                                   agcGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAg                                                <  2:922263/100‑1 (MQ=255)
                                                   agcGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAg                                                <  2:305983/100‑1 (MQ=255)
                                                    gcGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAgg                                               <  2:422165/100‑1 (MQ=255)
                                                    gcGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAgg                                               <  1:2064241/100‑1 (MQ=255)
                                                    gcGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAgg                                               <  2:2353041/100‑1 (MQ=255)
                                                        ttGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATa                                           <  1:1186047/100‑1 (MQ=255)
                                                           aTGGTCTTGAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATc                                        >  1:2108748/1‑100 (MQ=9)
                                                           aTGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATc                                        <  2:2180602/100‑1 (MQ=255)
                                                           aTGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATc                                        >  2:565843/1‑100 (MQ=255)
                                                            tGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCg                                       <  2:2510801/100‑1 (MQ=255)
                                                            tGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCg                                       >  2:2799209/1‑100 (MQ=255)
                                                             ggTCTTGAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCgg                                      >  1:1344474/1‑100 (MQ=255)
                                                              gTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCggg                                     >  2:2721830/1‑100 (MQ=255)
                                                              gTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCggg                                     >  2:501139/1‑100 (MQ=255)
                                                              gTCTTGAACAATTGACGGGCCAGTTGATGCTGCTAGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGAACCGGCAGGGCGCTATCCAGGGATAATCggg                                     >  2:2162188/1‑100 (MQ=255)
                                                               tCTTGAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCggg                                     >  2:2368131/1‑99 (MQ=12)
                                                                cTTGAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCggg                                     >  1:2552167/1‑98 (MQ=12)
                                                                cTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCa                                   <  1:1807375/100‑1 (MQ=255)
                                                                cTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCa                                   >  1:973275/1‑100 (MQ=255)
                                                                cTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCa                                   <  1:1600728/100‑1 (MQ=255)
                                                                cTTGAACAATTGACGGGCCAGTTGATGCGGCTCGAGCGGGGGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCa                                   <  1:638199/100‑1 (MQ=255)
                                                                 ttGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCaa                                  <  1:1082721/100‑1 (MQ=255)
                                                                 ttGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCaa                                  >  1:2111522/1‑100 (MQ=255)
                                                                   gAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCCAAc                                >  1:2286446/1‑100 (MQ=11)
                                                                   gAACAATTGGCGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAAc                                >  2:1322164/1‑100 (MQ=255)
                                                                   gAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAGATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAAc                                >  2:11769/1‑100 (MQ=255)
                                                                    aaCAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACa                               >  2:2347234/1‑100 (MQ=255)
                                                                    aaCAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACa                               >  1:1004691/1‑100 (MQ=255)
                                                                      cAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGGAAACAgg                             >  2:16061/1‑100 (MQ=255)
                                                                      cAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAgg                             >  2:1931517/1‑100 (MQ=255)
                                                                      cAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAgg                             >  2:37312/1‑100 (MQ=255)
                                                                      cAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAgg                             <  1:538062/100‑1 (MQ=255)
                                                                      cAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAgg                             <  2:183410/100‑1 (MQ=255)
                                                                         ttggcGGGCCAGTTGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGgcgc                          <  2:670586/96‑1 (MQ=255)
                                                                         ttGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGgcgc                          <  2:102525/100‑1 (MQ=255)
                                                                            gcGGGCCAGTTGATGCTGCTCCCGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATg                       >  1:1690108/2‑100 (MQ=255)
                                                                             cGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATgg                      >  1:505387/1‑100 (MQ=255)
                                                                             cGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATgg                      >  1:225364/1‑100 (MQ=255)
                                                                               ggCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAg                    >  2:1535324/1‑100 (MQ=255)
                                                                               ggCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAg                    >  2:937072/1‑100 (MQ=255)
                                                                                gCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAgg                   >  2:2191884/1‑100 (MQ=255)
                                                                                   aGTTGATGCTGCTCGAGCAGGTGGCGGAAATTTATGATGTTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGa                <  1:2046788/100‑1 (MQ=255)
                                                                                      tGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGAttt             >  2:977326/1‑100 (MQ=255)
                                                                                      tGATGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGAttt             <  1:1352162/100‑1 (MQ=255)
                                                                                       gATGCTGCTCGCGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTc            >  2:2837408/1‑100 (MQ=255)
                                                                                        aTGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCg           <  2:505390/100‑1 (MQ=255)
                                                                                        aTGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCg           <  1:1660045/100‑1 (MQ=255)
                                                                                        aTGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCg           <  2:251704/100‑1 (MQ=255)
                                                                                        aTGCTGCTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCg           <  1:457525/100‑1 (MQ=255)
                                                                                         tgctgcTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCg           <  1:76919/99‑1 (MQ=255)
                                                                                         tgctgcTCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGt          >  1:944791/1‑100 (MQ=255)
                                                                                         tgctgcTAGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGt          >  2:1526648/1‑100 (MQ=255)
                                                                                          gctgctCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTa         >  1:812572/1‑100 (MQ=255)
                                                                                           ctgctCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATGTCGTAc        >  2:1415753/1‑100 (MQ=255)
                                                                                            tgctCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACa       <  2:277303/100‑1 (MQ=255)
                                                                                             gctCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTCCAg      >  1:2707784/1‑100 (MQ=255)
                                                                                             gctCCAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATGACCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAg      >  1:1353063/1‑100 (MQ=255)
                                                                                                cGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGc   >  1:39218/1‑100 (MQ=255)
                                                                                                 gAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCa  >  1:2818319/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
TGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGATGGTGGTGCGATCCGGCAGGGCGCTATCCAGGGATAATCGGGCAAACAGGCGCATGGAGGCGATTTCGTACAGGGCA  >  NC_000913/1428294‑1428490

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