Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 308,649 A→G V201A (GTG→GCG)  yagY ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913308,6490AG99.4% 1018.6 / ‑4.9 321V201A (GTG→GCG) yagYhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (163/155);  total (165/156)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 344 total aligned reads displayed.

TCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAG  >  NC_000913/308559‑308738
                                                                                          |                                                                                         
  gTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCcgcgcg                                                                                       >  1:5012181/1‑93 (MQ=255)
   ttATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCgg                                                                                      >  1:4771625/1‑93 (MQ=255)
   ttATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCgg                                                                                      >  2:7409917/1‑93 (MQ=255)
   ttATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCgg                                                                                      >  1:7304174/1‑93 (MQ=255)
     aTCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGt                                                                                    >  1:141240/1‑93 (MQ=255)
     aTCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGt                                                                                    >  1:8943147/1‑93 (MQ=255)
      tCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTa                                                                                   <  1:9887333/93‑1 (MQ=255)
      tCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTa                                                                                   >  2:3210843/1‑93 (MQ=255)
      tCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTa                                                                                   >  1:5409405/1‑93 (MQ=255)
       ccATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTaa                                                                                  <  2:2722754/93‑1 (MQ=255)
       ccATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTaa                                                                                  >  2:7617978/1‑93 (MQ=255)
       ccATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTaa                                                                                  >  1:4893060/1‑93 (MQ=255)
         aTCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTaaaa                                                                                <  2:9526990/93‑1 (MQ=255)
          tCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGTAACCCGTCCTTTGTTATCCGCAGTGTCCGTGCGGGTAATac                                                                               >  2:8951781/1‑93 (MQ=255)
          tCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAac                                                                               >  1:102834/1‑93 (MQ=255)
          tCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAac                                                                               >  2:8848628/1‑93 (MQ=255)
          tCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAac                                                                               >  2:4300539/1‑93 (MQ=255)
           cTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGAGCGGGTAAAacg                                                                              >  1:4705081/1‑93 (MQ=255)
            tGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAacga                                                                             <  2:1652003/93‑1 (MQ=255)
            tGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAacga                                                                             <  1:3697705/93‑1 (MQ=255)
            ctcccaGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTTCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGGAAAacga                                                                             <  1:4426938/89‑1 (MQ=255)
             gCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAacgac                                                                            <  1:614291/93‑1 (MQ=255)
               acaGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGAcgc                                                                          >  2:3766005/1‑93 (MQ=255)
                  gCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGAc                                                                       >  2:1774931/1‑93 (MQ=255)
                  gCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGAc                                                                       <  2:4200239/93‑1 (MQ=255)
                  gCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGAc                                                                       <  1:7019295/93‑1 (MQ=255)
                   ctatctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACt                                                                      <  2:6509265/93‑1 (MQ=255)
                   ctatctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACt                                                                      >  1:8185005/1‑93 (MQ=255)
                   ctatctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACt                                                                      <  1:1091575/93‑1 (MQ=255)
                   ctatctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGCCCGCGCGGGTAAAACGACGCGCCt                                                                      >  2:7686340/1‑93 (MQ=255)
                    tatctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACtt                                                                     >  1:6481542/1‑93 (MQ=255)
                    tatctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACtt                                                                     >  2:4651936/1‑93 (MQ=255)
                     atctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAGCGACGCGACTTg                                                                    >  1:1933137/1‑93 (MQ=255)
                     atctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTg                                                                    >  1:5738329/1‑93 (MQ=255)
                      tctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGc                                                                   >  1:5197867/1‑93 (MQ=255)
                       ctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTTcc                                                                  >  2:9263550/1‑93 (MQ=255)
                       ctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGcc                                                                  <  2:5255945/93‑1 (MQ=255)
                       ctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGcc                                                                  <  2:615162/93‑1 (MQ=255)
                        tatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGccc                                                                 >  1:9642129/1‑93 (MQ=255)
                          tttCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCgg                                                               >  1:452368/1‑93 (MQ=255)
                          tttCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCgg                                                               >  1:794963/1‑93 (MQ=255)
                           ttCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGc                                                              >  1:9969046/1‑93 (MQ=255)
                            tCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCa                                                             >  1:7370971/1‑93 (MQ=255)
                            tCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCa                                                             <  1:10018205/93‑1 (MQ=255)
                            tCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCa                                                             <  1:2385071/93‑1 (MQ=255)
                              aCGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATc                                                           <  1:1825250/93‑1 (MQ=255)
                              aCGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATc                                                           >  2:6823248/1‑93 (MQ=255)
                                 ggAATGAACTTATCGCCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAgg                                                        <  1:6585486/93‑1 (MQ=255)
                                 ggAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAgg                                                        >  2:1872023/1‑93 (MQ=255)
                                  gAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGCAAAACGACGCGACTTGCCCGGCATCAGgt                                                       >  1:9236432/1‑93 (MQ=255)
                                     tGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgtagt                                                    >  2:2930692/1‑93 (MQ=255)
                                       aaCTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTaa                                                  <  1:3293362/93‑1 (MQ=255)
                                          ttATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtc                                               <  2:9912724/93‑1 (MQ=255)
                                          ttATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtc                                               >  1:7645408/1‑93 (MQ=255)
                                            aTCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtctc                                             >  1:3041123/1‑93 (MQ=255)
                                              cACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCtt                                           >  2:7575781/1‑93 (MQ=255)
                                                cccTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCNGTGTCCGCGCGNGTAAAACGACGCGACNTGCCCNGCATCAGGTAGTAATTCTNTTTa                                         >  1:4156010/1‑93 (MQ=255)
                                                cccTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACTACGCGACTTGCCCGGCATCAGGGAGTAATTCTCTTTa                                         >  2:8699585/1‑93 (MQ=255)
                                                   tGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTacac                                      <  2:8223783/93‑1 (MQ=255)
                                                    gCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACt                                     <  1:4360833/93‑1 (MQ=255)
                                                    gCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACt                                     <  2:5506017/93‑1 (MQ=255)
                                                    gCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACt                                     <  2:3054431/93‑1 (MQ=255)
                                                      cAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTcc                                   >  2:6401431/1‑93 (MQ=255)
                                                      cAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTcc                                   <  1:8038393/93‑1 (MQ=255)
                                                      cAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTcc                                   >  2:4116863/1‑93 (MQ=255)
                                                       ataGTGCAGCCCGTCGTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCt                                  <  2:3703084/91‑1 (MQ=255)
                                                        aaGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCtt                                 >  2:7106027/1‑93 (MQ=255)
                                                         aGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTa                                >  2:7336825/1‑93 (MQ=255)
                                                         aGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTa                                <  1:2088047/93‑1 (MQ=255)
                                                           tGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTAcc                              <  1:7485613/93‑1 (MQ=255)
                                                           tGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTAcc                              >  2:1678158/1‑93 (MQ=255)
                                                            gCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCg                             >  1:684888/1‑93 (MQ=255)
                                                            gCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCg                             >  2:502355/1‑93 (MQ=255)
                                                              aaCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGtt                           <  2:2084677/93‑1 (MQ=255)
                                                              aaCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGtt                           >  2:1167021/1‑93 (MQ=255)
                                                                cccGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTgg                         >  1:846559/1‑93 (MQ=255)
                                                                cccGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTgg                         >  2:9684783/1‑93 (MQ=255)
                                                                  cGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggcg                       >  2:2636317/1‑93 (MQ=255)
                                                                  cGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggcg                       <  2:6653753/93‑1 (MQ=255)
                                                                   gTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggcgg                      >  1:3835707/1‑93 (MQ=255)
                                                                    tCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggcggc                     >  1:2820238/1‑93 (MQ=255)
                                                                     ccTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCt                    <  2:1405766/93‑1 (MQ=255)
                                                                       tttGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCttt                  <  2:8487291/93‑1 (MQ=255)
                                                                        ttGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTc                 >  2:6074859/1‑93 (MQ=255)
                                                                        ttGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTc                 <  1:5351211/93‑1 (MQ=255)
                                                                           ttATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCGGTTGGTAATTCTTTTTCCACTCCTTACGGTTGCGGGTTTTcggc              >  1:3457069/1‑90 (MQ=255)
                                                                           ttACCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTTTTTACACTCCTTACCGTTGGCGGCTTTCAgg              <  2:7794956/93‑1 (MQ=255)
                                                                             aTCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCa            <  2:5903945/93‑1 (MQ=255)
                                                                             aTCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCa            >  2:5734107/1‑93 (MQ=255)
                                                                             aTCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCa            >  1:6369355/1‑93 (MQ=255)
                                                                             aTCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCa            >  1:7542644/1‑93 (MQ=255)
                                                                              tCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTTCACTCCTTTCCGTTGGCGGCTTTCAGGCaa           <  2:9411708/93‑1 (MQ=255)
                                                                               ccGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAg          <  1:2210879/93‑1 (MQ=255)
                                                                               ccGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAg          <  1:967502/93‑1 (MQ=255)
                                                                                cGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAgg         <  2:3274541/93‑1 (MQ=255)
                                                                                cGCAGTGGCGACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTAGCGTTGGCGGCTTTCAGGCAAgg         <  2:3998967/93‑1 (MQ=255)
                                                                                  cAGTGTCGGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTc       <  1:1492112/93‑1 (MQ=255)
                                                                                  cAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTc       <  2:1237877/93‑1 (MQ=255)
                                                                                    gtgtCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCg     >  1:4574618/1‑93 (MQ=255)
                                                                                    gtgtCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCg     >  2:9179940/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAG  >  NC_000913/308559‑308738

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