Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,572,928 C→G 15.3% A759P (GCA→CCA)  pqqL ← putative periplasmic M16 family zinc metalloendopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,572,9280CG15.3% 366.5 / 58.2 144A759P (GCA→CCA) pqqLputative periplasmic M16 family zinc metalloendopeptidase
Reads supporting (aligned to +/- strand):  ref base C (69/52);  new base G (5/17);  total (75/69)
Fisher's exact test for biased strand distribution p-value = 4.60e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.60e-01

Only 100 of 152 total aligned reads displayed.

GATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGA  >  NC_000913/1572829‑1573027
                                                                                                   |                                                                                                   
gATCTTTTGCCAGTGCGACGTTAAACCCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGc                                                                                                     >  1:1712169/1‑100 (MQ=255)
gATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGc                                                                                                     >  2:1515748/1‑100 (MQ=255)
gATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGc                                                                                                     >  1:2338273/1‑100 (MQ=255)
 aTCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGcc                                                                                                    >  2:575624/1‑100 (MQ=255)
 aTCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGcc                                                                                                    >  2:2826866/1‑100 (MQ=255)
 aTCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGcc                                                                                                    >  2:2001398/1‑100 (MQ=255)
   cTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCcac                                                                                                  >  2:826697/1‑100 (MQ=255)
   cTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCcac                                                                                                  >  1:326160/1‑100 (MQ=255)
        gCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCcacaag                                                                                               >  2:3984405/1‑96 (MQ=255)
           aGTGCGACGTTAAATGTATTGAGTGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTCTGAAACCTGTGGCACAGGtttt                                                                                           <  1:107080/99‑3 (MQ=255)
           aGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttca                                                                                          <  2:449769/100‑4 (MQ=255)
            gTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcac                                                                                         <  2:979465/100‑5 (MQ=255)
            gTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcac                                                                                         <  2:2653616/100‑5 (MQ=255)
            gTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcac                                                                                         <  1:533217/100‑5 (MQ=255)
            gTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcac                                                                                         <  1:166701/100‑5 (MQ=255)
            gTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtt                                                                                         >  1:1372521/1‑100 (MQ=255)
             tGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcaca                                                                                        <  1:686681/100‑6 (MQ=255)
             tGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcaca                                                                                        <  2:1374574/100‑6 (MQ=255)
             tGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAttt                                                                                        <  2:511199/100‑1 (MQ=255)
              gCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCCGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcacag                                                                                       <  1:3328832/100‑7 (MQ=255)
              gCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcacag                                                                                       <  1:265850/100‑7 (MQ=255)
              gCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtttt                                                                                       <  2:22266/100‑1 (MQ=255)
              gCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtttt                                                                                       <  1:3287076/100‑1 (MQ=255)
              gCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCAtttt                                                                                       >  1:4146160/1‑100 (MQ=255)
               cGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTg                                                                                      >  1:2698702/1‑100 (MQ=255)
                gACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcacagtg                                                                                     <  1:2362367/100‑9 (MQ=255)
                gACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGGCACAGGtttcacagtg                                                                                     <  1:1160451/100‑9 (MQ=255)
                  cGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTc                                                                                   <  1:1308089/100‑1 (MQ=255)
                   tttAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCt                                                                                  <  2:603232/99‑1 (MQ=255)
                    ttAAAAGCATCGAGCGCCATACGCGTCTGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCt                                                                                  <  1:2095346/99‑1 (MQ=255)
                    ttAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCta                                                                                 >  1:603230/1‑99 (MQ=255)
                     tAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCttt                                                                                >  2:3427552/1‑100 (MQ=255)
                      aaaaGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCtttt                                                                               <  2:285606/100‑1 (MQ=255)
                       aaaGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTa                                                                              >  1:690396/1‑100 (MQ=255)
                       aaaGCATCGAGCGCCATACGCGACGGCAGATTAACAGGTGTCCGGGAAGCAGAACGCGTCCACTGTGAAACCAGTGCCACAGATACATGTTGTTCTTTTa                                                                              >  1:3322419/1‑100 (MQ=255)
                        aaGGATCGAGCGCCATACGCGTCGTCACATTAACATGTGTCCGGTAATCATAACGCTTCCTCTGTGAAACCTGTGCCACAGGTTCATTTTGTTCtttt                                                                               <  1:3390628/98‑1 (MQ=255)
                        aaGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTAc                                                                             >  2:4384320/1‑100 (MQ=255)
                        aaGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTAc                                                                             >  1:1837947/1‑100 (MQ=255)
                         aGCATCGAGCGCCATACGCGTCGGCAGATTAAAAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACa                                                                            <  2:1147425/100‑1 (MQ=255)
                          gCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAg                                                                           >  1:2225841/1‑100 (MQ=255)
                           cATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACTGGTTCATTTTGTTCTTTTACatt                                                                          >  1:1325351/1‑98 (MQ=255)
                           cATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGt                                                                          >  2:1050040/1‑100 (MQ=255)
                            aTCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTa                                                                         >  1:2030974/1‑100 (MQ=255)
                             tCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTaa                                                                        <  1:4206985/100‑1 (MQ=255)
                              cGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAAc                                                                       >  2:1736866/1‑100 (MQ=255)
                                aGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCg                                                                     >  2:1949617/1‑100 (MQ=255)
                                    ccATACGCGTAGGCAGATTAACAGGTGTCCGGGAATCATACCGCTTACACTGTGAAACTTGTGACAAAGGTTCATTTTGTTCTTTTAAAGTAACCgtt                                                                   >  1:1573085/1‑96 (MQ=255)
                                       tACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGtt                                                              <  1:2851621/100‑1 (MQ=255)
                                       tACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGtt                                                              >  1:1176309/1‑100 (MQ=255)
                                          gcgTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTc                                                           <  2:3254782/100‑1 (MQ=255)
                                          gcgTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTc                                                           <  2:728981/100‑1 (MQ=255)
                                              cGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTc                                                       >  2:3143543/1‑100 (MQ=255)
                                               ggCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcg                                                      <  2:2082686/100‑1 (MQ=255)
                                               ggCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcg                                                      <  1:1349011/100‑1 (MQ=255)
                                                gCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcgc                                                     >  1:3960077/1‑100 (MQ=255)
                                                  agattccccggTGTACGGAAATCATAACGCTTCCAATATACCACCTCTGCCACAGGTTCATTCTGTGCTTTCACAGTAACCGAAGCGTTGTCCGTcgcga                                                   >  1:3499554/10‑99 (MQ=255)
                                                  aGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcgcgc                                                   >  2:1386513/1‑100 (MQ=255)
                                                   gATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcgcgcg                                                  <  2:223374/100‑1 (MQ=255)
                                                   gATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTcgcgcg                                                  >  2:3283296/1‑100 (MQ=255)
                                                     ttAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAg                                                <  2:3935013/100‑1 (MQ=255)
                                                     ttAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAg                                                <  1:1644195/100‑1 (MQ=255)
                                                        aCAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTa                                             <  1:1010508/100‑1 (MQ=255)
                                                           ggTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATg                                          >  2:718695/1‑100 (MQ=255)
                                                           ggTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATg                                          >  2:3674827/1‑100 (MQ=255)
                                                            gtgtCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATgg                                         <  1:3067769/100‑1 (MQ=255)
                                                             tgtCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGATCTTTGACAGTAACCGATGCGTTGTCCGTCGCGCGAGGTAATGGt                                        <  1:2043846/100‑1 (MQ=255)
                                                               tCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGttt                                      >  2:3624255/1‑100 (MQ=255)
                                                                ccGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTa                                     >  1:4088461/1‑100 (MQ=255)
                                                                  gggAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTAcc                                   <  1:4239800/100‑1 (MQ=255)
                                                                  gggAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTAcc                                   <  1:3849377/100‑1 (MQ=255)
                                                                  gggAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTAcc                                   >  1:1489446/1‑100 (MQ=255)
                                                                        cATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCgg                              >  1:1487823/1‑99 (MQ=255)
                                                                        cATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGc                             >  2:4064334/1‑100 (MQ=255)
                                                                         aTAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCt                            >  2:2592928/1‑100 (MQ=255)
                                                                            aCGCTTCCACGGTGAAACCTGTGCCACAGTTTTATTTTGTTCTTTTACAGTAACCGATGCTTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTaa                          <  1:1226284/99‑1 (MQ=255)
                                                                             cGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATg                        >  2:1541879/1‑100 (MQ=255)
                                                                             cGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATg                        >  2:387593/1‑100 (MQ=255)
                                                                             cGCTTCCACTATGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATg                        >  1:3000337/1‑100 (MQ=255)
                                                                               cTTCCACTGTGAAACCTGTGCCACAGTGTCATTTTGTACTTTTCCAGTAATCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGc                      <  1:3507484/100‑1 (MQ=255)
                                                                                   cACTTTGAAACCTGTGCCCCAGGTTGATCTTGTTCGTTTACAGTAACCGATGCGGTGTTCGCCGCGCGGCTTAAGGCTTTACCTTCGGCTAATGTCGAAt                  <  1:3930303/100‑1 (MQ=255)
                                                                                   cACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAAt                  >  2:3452432/1‑100 (MQ=255)
                                                                                   acaaaatgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAAt                  >  2:2998820/9‑100 (MQ=255)
                                                                                    aCTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATc                 <  2:4006083/100‑1 (MQ=255)
                                                                                     cTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCa                >  1:3738342/1‑100 (MQ=255)
                                                                                     aaaatgaaCCTGTGGCACAGGTTCAGTTTGTTCTTTTACAGTAACCTATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCa                >  1:1970378/7‑100 (MQ=255)
                                                                                      tgtgAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCag               >  2:117081/1‑100 (MQ=255)
                                                                                       gtgAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACGGATGCGTTGTCCGTCGGGCGAGTTAGTGGTTTACCTGCGGCTAATGGCGAATCaga              >  1:3089611/1‑100 (MQ=255)
                                                                                       gtgAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGaagatcg              >  1:1846164/1‑95 (MQ=255)
                                                                                       gtgAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGGCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCag               <  1:1214355/99‑1 (MQ=255)
                                                                                        atgaaCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTCATGGCGAGTccg               >  1:4225599/4‑96 (MQ=255)
                                                                                         gAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAgt            >  2:3654060/1‑100 (MQ=255)
                                                                                          aaaCCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAgtg           <  2:531708/100‑1 (MQ=255)
                                                                                          aaaCCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATCCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAGTCagcg             >  1:2129419/1‑96 (MQ=255)
                                                                                              cTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttg       >  1:1902691/1‑100 (MQ=255)
                                                                                              cTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttg       <  1:3311396/100‑1 (MQ=255)
                                                                                              cTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttg       >  2:1672033/1‑100 (MQ=255)
                                                                                              cTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttg       <  2:302631/100‑1 (MQ=255)
                                                                                               tgtgCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttga      >  1:3635085/1‑100 (MQ=255)
                                                                                                 tgCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgatn    <  2:838131/100‑2 (MQ=255)
                                                                                                  gCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTttgattg   >  2:2294965/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGA  >  NC_000913/1572829‑1573027

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: