Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,921,059 T→G P279P (CCT→CCG rsmF → 16S rRNA m(5)C1407 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,921,0590TG97.6% 415.9 / 0.9 124P279P (CCT→CCGrsmF16S rRNA m(5)C1407 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (63/57);  total (65/59)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

Only 100 of 126 total aligned reads displayed.

GTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCTGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACT  >  NC_000913/1920962‑1921155
                                                                                                 |                                                                                                
gTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCggg                                                                                                <  1:456199/100‑1 (MQ=255)
 tACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGt                                                                                               <  1:237781/100‑1 (MQ=255)
 tACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGt                                                                                               >  1:3319576/1‑100 (MQ=255)
  aCCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTg                                                                                              <  2:374673/100‑1 (MQ=255)
  aCCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTg                                                                                              >  1:1809824/1‑100 (MQ=255)
  aCCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTg                                                                                              >  2:787544/1‑100 (MQ=255)
  aCCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTg                                                                                              >  1:3716605/1‑100 (MQ=255)
   ccTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGc                                                                                             <  2:2111457/100‑1 (MQ=255)
   ccTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGc                                                                                             >  1:1155030/1‑100 (MQ=255)
   ccTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGc                                                                                             >  1:2458312/1‑100 (MQ=255)
   ccTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGGTTTTACCACTTGGCGATCTCTTCCCGGGTGc                                                                                             >  1:347011/1‑100 (MQ=255)
    cTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGca                                                                                            >  2:408999/1‑100 (MQ=255)
      tAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGcaaa                                                                                          >  1:2026800/1‑100 (MQ=255)
       aaaCCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGcaaac                                                                                         <  1:3697024/100‑1 (MQ=255)
       aaaCCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGcaaac                                                                                         <  2:1988481/100‑1 (MQ=255)
       aaaCCACGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACGTACCTCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGgaaac                                                                                         <  2:3764854/100‑1 (MQ=255)
        aatcgggAAGAAAAATAAGCCGCTTGCCTGTGGTTGAAAGAGACTTACCTCGACGCAGTAGAGTTCTTACCACTTGGCGATCTCTTCCCGCGTTcaaaca                                                                                        >  2:3089114/6‑100 (MQ=255)
          ccAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCTGGTGcaaacaaa                                                                                      >  1:993386/1‑100 (MQ=255)
           cAGGAAGAAAACGAAGCCGTTTGTCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAg                                                                                     >  2:3587363/1‑100 (MQ=255)
           cAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAg                                                                                     >  2:2378799/1‑100 (MQ=255)
            aGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAgc                                                                                    >  1:287989/1‑100 (MQ=255)
              gaagaaAACGAAGCCGTTTGCCTTTGGCTGAAAGAGACTTACCCCAACGCAGTAGATTTTTTACCACTTGGCGATCTCTTCCCGGCTGCACACAAAgcgc                                                                                  <  2:117268/100‑1 (MQ=255)
               aagaaAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCt                                                                                 >  1:1612299/1‑100 (MQ=255)
               aagaaAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCt                                                                                 >  1:4016335/1‑100 (MQ=255)
                  aaaaCGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGAc                                                                              <  1:1023279/100‑1 (MQ=255)
                  aaaaCGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGAc                                                                              <  1:564104/100‑1 (MQ=255)
                   aaaCGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGAcc                                                                             <  2:2207034/100‑1 (MQ=255)
                     aCGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCga                                                                           <  1:2018947/100‑1 (MQ=255)
                       gAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCgaag                                                                         <  1:124893/100‑1 (MQ=255)
                        aaGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCgaaga                                                                        >  1:3283263/1‑100 (MQ=255)
                           gcGTTTGCCCGTGTCTGAAAGCGACTTACACCGACGCAGTAGAGTATTTGCCAATTGGCGCCCTCTTCCCGGGTGCAAACAAAGCGCTGACCCaagaagc                                                                     <  2:546871/99‑2 (MQ=255)
                                tGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCtttt                                                                <  1:2404303/100‑1 (MQ=255)
                                 gCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCttttt                                                               >  1:3285886/1‑100 (MQ=255)
                                 gCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGAGCAAACAAAGCGCTCACCGAAGAATGCttttt                                                               >  1:3744706/1‑100 (MQ=255)
                                  ccTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTAGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTg                                                              >  2:3641576/1‑100 (MQ=255)
                                    tgtgGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCa                                                            >  1:3493042/1‑100 (MQ=255)
                                     gtgGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCAt                                                           <  1:1227099/100‑1 (MQ=255)
                                     gtgGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCAt                                                           >  1:2944377/1‑100 (MQ=255)
                                       ggCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAGGCGCTGACCGAAGCAGGCTTTTTGCATGt                                                         >  2:3966276/1‑100 (MQ=255)
                                        gCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGtt                                                        <  2:766711/100‑1 (MQ=255)
                                         cTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGttt                                                       >  1:1744552/1‑100 (MQ=255)
                                         cTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGttt                                                       >  1:1747212/1‑100 (MQ=255)
                                         cTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGttt                                                       >  2:102460/1‑100 (MQ=255)
                                         cTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGttt                                                       >  1:873351/1‑100 (MQ=255)
                                          tGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAGGAAGGCTTTGTGCATGtttt                                                      >  2:875454/1‑100 (MQ=255)
                                          tGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGtttt                                                      <  2:3190578/100‑1 (MQ=255)
                                             aaGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTccc                                                   <  2:1249969/100‑1 (MQ=255)
                                             aaGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTccc                                                   <  2:2082929/100‑1 (MQ=255)
                                               gagaCTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCcac                                                 >  2:3677029/1‑100 (MQ=255)
                                               gagaCTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCcac                                                 <  1:24837/100‑1 (MQ=255)
                                               gagaCTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCcac                                                 <  1:1479863/100‑1 (MQ=255)
                                                 gaCTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACaa                                               >  1:214412/1‑100 (MQ=255)
                                                 gaCTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACaa                                               >  2:3699980/1‑100 (MQ=255)
                                                   cTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAAt                                             >  1:3093849/1‑100 (MQ=255)
                                                   cTTACCCCGACGCAGTAGAGTGTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAAt                                             >  2:2405202/1‑100 (MQ=255)
                                                    ttACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAAtt                                            <  2:3051140/100‑1 (MQ=255)
                                                       ncccGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTAc                                         >  1:113379/2‑100 (MQ=255)
                                                       ccccGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTAc                                         >  2:875906/1‑100 (MQ=255)
                                                       ccccGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTAc                                         >  2:665557/1‑100 (MQ=255)
                                                       ccccGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTAc                                         <  1:1797689/100‑1 (MQ=255)
                                                       ccccGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTAc                                         >  1:3294849/1‑100 (MQ=255)
                                                       ccccGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTAc                                         >  2:3525865/1‑100 (MQ=255)
                                                            aCGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTg                                    <  1:313346/100‑1 (MQ=255)
                                                              gCAGTCGAGTTTTTACCACTTGGCCATCTATTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAAGTTACGACTGCg                                  <  2:622384/100‑1 (MQ=255)
                                                                aGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGaa                                >  1:4020777/1‑100 (MQ=255)
                                                                 gTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAg                               <  2:503354/100‑1 (MQ=255)
                                                                 gTAGAGTTTTTACCACTTGGCGATCTCTTCCCAGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAg                               >  2:3595544/1‑100 (MQ=255)
                                                                  tAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAgg                              <  2:636076/100‑1 (MQ=255)
                                                                  tAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAgg                              <  2:2978075/100‑1 (MQ=255)
                                                                   agagTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGc                             >  1:492924/1‑100 (MQ=255)
                                                                   agagTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGc                             <  1:1670663/100‑1 (MQ=255)
                                                                    gagTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACATAGCGCTGACCGAAGAAGGCTTTTTGCATGTTATCCCACAAATTTACGACTGCGAAGGct                            >  1:3535511/1‑100 (MQ=255)
                                                                    gagTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGct                            <  2:2974160/100‑1 (MQ=255)
                                                                     agTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGctt                           <  1:1654813/100‑1 (MQ=255)
                                                                      gTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGcttc                          >  2:155731/1‑100 (MQ=255)
                                                                      gTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGcttc                          >  1:3992968/1‑100 (MQ=255)
                                                                      gTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATCTACGACTGCGAAGGcttc                          >  1:3829665/1‑100 (MQ=255)
                                                                       tttttACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGcttct                         <  1:804467/100‑1 (MQ=255)
                                                                         tttACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCttcttc                       >  1:610353/1‑100 (MQ=255)
                                                                         tttACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCttcttc                       <  2:2644459/100‑1 (MQ=255)
                                                                         tttACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCttcttc                       <  1:66220/100‑1 (MQ=255)
                                                                          ttACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCg                      <  2:773657/100‑1 (MQ=255)
                                                                           tACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGt                     <  2:907709/100‑1 (MQ=255)
                                                                            aCCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGtt                    <  1:1138222/100‑1 (MQ=255)
                                                                              cACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGc                  <  1:2738325/100‑1 (MQ=255)
                                                                              cACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGc                  <  2:3235/100‑1 (MQ=255)
                                                                               aCTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCt                 <  1:3425368/100‑1 (MQ=255)
                                                                               aCTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCt                 <  1:1311390/100‑1 (MQ=255)
                                                                               aCTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCt                 >  1:687737/1‑100 (MQ=255)
                                                                                 ttGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCg               <  1:667568/100‑1 (MQ=255)
                                                                                     cGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTg           >  2:3100265/1‑100 (MQ=255)
                                                                                       aTCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCg         <  1:3056528/100‑1 (MQ=255)
                                                                                        tctcTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGt        <  1:1650792/100‑1 (MQ=255)
                                                                                         ctctTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTTTTTTTTTTTCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTa       >  2:754039/1‑100 (MQ=255)
                                                                                         ctctTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTa       >  1:18749/1‑100 (MQ=255)
                                                                                         ctctTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTAAGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTa       >  2:4694/1‑100 (MQ=255)
                                                                                         ctctTCCCGGATGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAACTTTACGACTGCGAAGGCATCTTCGTTGCGCGTCTGCgga       <  2:1555722/100‑3 (MQ=255)
                                                                                          tctTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTaa      <  1:4031670/100‑1 (MQ=255)
                                                                                           ctTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTaaa     >  1:3189771/1‑100 (MQ=255)
                                                                                             tCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAAc   <  2:2354703/100‑1 (MQ=255)
                                                                                                 |                                                                                                
GTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCTGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACT  >  NC_000913/1920962‑1921155

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