Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 547,835 +G pseudogene (4/1008 nt) ylbE → putative protein, C‑ter fragment (pseudogene)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913547,8311.G98.0% 1748.3 / 24.5 450pseudogene (251/252 nt)ylbEputative protein, C‑ter fragment (pseudogene)
Reads supporting (aligned to +/- strand):  ref base . (3/6);  new base G (241/200);  total (244/206)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 525 total aligned reads displayed.

TGGCAGGAGATGACCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTG  >  NC_000913/547731‑547930
                                                                                                     |                                                                                                     
  gCAGGAGATGACCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgg                                                                                                    <  1:4167451/101‑1 (MQ=255)
   nnngagaTGACCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggg                                                                                                  <  2:8384641/98‑1 (MQ=255)
   gaGGAGATGACCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCannn                                                                                                     >  2:9102003/2‑98 (MQ=255)
    aagagaTGACCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgnnn                                                                                                   >  2:7777032/3‑98 (MQ=255)
    aGGAGATGACCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgnnn                                                                                                   >  2:8489396/1‑98 (MQ=255)
    aGGAGATGACCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgnnn                                                                                                   >  2:12372210/1‑98 (MQ=255)
    aGGAGATGACCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggg                                                                                                 >  1:2094435/1‑101 (MQ=255)
      nnnaTGACCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGACCAggggga                                                                                               <  2:13349260/98‑3 (MQ=255)
      gagaTGACCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCAAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAggggga                                                                                               <  1:8385616/101‑3 (MQ=255)
       agaTGACCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggggaa                                                                                              <  2:1004590/101‑4 (MQ=255)
       agaTGACCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggggaa                                                                                              >  1:10007858/1‑98 (MQ=255)
        gaTGACCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggggnnn                                                                                             >  2:8045520/1‑97 (MQ=255)
        gaTGACCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggggaag                                                                                             >  1:6700061/1‑97 (MQ=255)
        gaTGAACGGACCCAAGAAATGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAgggggaag                                                                                             <  2:1953441/101‑5 (MQ=255)
            aCCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGaa                                                                                          >  1:10288649/1‑101 (MQ=255)
             ccGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGCCGCTGGCAATACTGGAGCAGGGGGAAGTGaat                                                                                         >  2:607297/1‑100 (MQ=255)
               ggACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGATGGGCGAAAGTTGAAGTGCAGGCGCTGGCAATACTGGAGCAGGTGTAAGTGaannn                                                                                       >  2:13700011/1‑98 (MQ=255)
                gtcccATGAAAGGGGCGGGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTc                                                                                      <  2:4658034/99‑1 (MQ=255)
                 nnncATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCa                                                                                     <  2:9463007/98‑1 (MQ=255)
                   nnntGACAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAtt                                                                                   <  2:13243719/98‑1 (MQ=255)
                   ccATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCAtt                                                                                   <  1:3438013/101‑1 (MQ=255)
                      tGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCt                                                                                <  1:8710331/101‑1 (MQ=255)
                       gAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                                                               <  1:8489396/101‑1 (MQ=255)
                       gAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCtt                                                                               <  1:3864014/101‑1 (MQ=255)
                        aaaGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTg                                                                              <  1:7873241/101‑1 (MQ=255)
                          aGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAACAGGGGGAAGTGAACTTCATTCCttnnn                                                                            >  2:10410906/1‑98 (MQ=255)
                               cgtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                                                                       >  2:175835/1‑101 (MQ=255)
                               cgtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                                                                       <  1:9332452/101‑1 (MQ=255)
                               cgtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAt                                                                       >  1:2364534/1‑101 (MQ=255)
                                gtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAATTTCATTCCTTGTCACCAtg                                                                      <  1:12454294/101‑1 (MQ=255)
                                gtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCGTTCCTTGTCACCAtg                                                                      >  1:13576292/1‑101 (MQ=255)
                                gtgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtg                                                                      <  1:14105469/101‑1 (MQ=255)
                                 tgcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgt                                                                     >  1:9650945/1‑101 (MQ=255)
                                  gcgtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtg                                                                    >  2:2393570/1‑101 (MQ=255)
                                    gtgGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCAtgtnnn                                                                  >  2:14097726/1‑98 (MQ=255)
                                        gagcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGTCGGGAACTTTATTCCTTGTCAGCATGTGAAtnnn                                                              >  2:11038146/3‑98 (MQ=255)
                                         cgcATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCg                                                             >  1:6258505/1‑101 (MQ=255)
                                           cATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTc                                                           >  2:4240551/1‑101 (MQ=255)
                                            nnntCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCg                                                          <  2:7388815/98‑1 (MQ=255)
                                             tGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCgg                                                         >  1:5116214/1‑101 (MQ=255)
                                               tCTGTTCGAAGGTTGGGCGAAAGATGAAGCTCAGGCGCTGGCAATACTGGATCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGc                                                       <  1:6433540/101‑1 (MQ=255)
                                                nnnttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc                                                      <  2:8615598/98‑1 (MQ=255)
                                                nnnttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGcc                                                      <  2:13966032/98‑1 (MQ=255)
                                                 tGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCa                                                     <  1:3638051/101‑1 (MQ=255)
                                                  gTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCaa                                                    >  1:8031656/1‑101 (MQ=255)
                                                  gTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCCTGTCACCATGTGAATGCCGTCGGGcnnn                                                    >  2:12936546/1‑98 (MQ=255)
                                                   ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGccnnn                                                   >  2:9716104/1‑98 (MQ=255)
                                                   ttCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAAt                                                   <  1:2725064/101‑1 (MQ=255)
                                                   nnngAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAAt                                                   <  2:11851300/98‑1 (MQ=255)
                                                   nnngAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAAt                                                   <  2:13628963/98‑1 (MQ=255)
                                                      gAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATggg                                                >  1:6912295/1‑101 (MQ=255)
                                                       aaGGTTGGGCGAAAGATTAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc                                               <  2:5740746/101‑1 (MQ=255)
                                                       aaGGTTGGGCGAAAGATGAAGCGCAGGGGGTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGAGCCAATgnnn                                               >  2:14507501/1‑98 (MQ=255)
                                                       aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc                                               >  2:4974439/1‑101 (MQ=255)
                                                       aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc                                               <  1:175835/101‑1 (MQ=255)
                                                       aaGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGc                                               <  2:704385/101‑1 (MQ=255)
                                                        aGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCg                                              >  1:11346871/1‑101 (MQ=255)
                                                         ggTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCgg                                             >  2:5179730/1‑101 (MQ=255)
                                                          gTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGt                                            >  1:12183895/1‑101 (MQ=255)
                                                           ttGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCggga                                           >  1:4349080/1‑99 (MQ=255)
                                                           ttGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTa                                           <  1:10410906/101‑1 (MQ=255)
                                                            tGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTAt                                          >  2:545459/1‑101 (MQ=255)
                                                              nnngAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTa                                        <  2:8832273/98‑1 (MQ=255)
                                                              ggCGAAAGATGAAGCGCAGGCGCTGGCAATACTGTAGCAGGGGGAAGGGAACTTCATTCGTTGTCACCATGTTAATGCCGTCGGGCCAGTGGGCGGTannn                                        >  2:10524696/1‑98 (MQ=255)
                                                               gCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc                                       >  1:5589150/1‑101 (MQ=255)
                                                               gCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTAc                                       >  1:11090360/1‑101 (MQ=255)
                                                                nnnaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACt                                      <  2:14597578/98‑1 (MQ=255)
                                                                nnnaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACt                                      <  2:7512335/98‑1 (MQ=255)
                                                                 gAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACtt                                     <  1:11752372/101‑1 (MQ=255)
                                                                  aaaGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTc                                    >  1:6982076/1‑101 (MQ=255)
                                                                    aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACttnnn                                  >  2:14378202/1‑98 (MQ=255)
                                                                    aGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGCGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTg                                  <  1:14538374/101‑1 (MQ=255)
                                                                        gtaGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGcnnn                              >  2:10894151/3‑98 (MQ=255)
                                                                        gAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAg                              <  1:1370247/101‑1 (MQ=255)
                                                                        gAAGCGCAGGCGCTGGCAAGACTGGAGCAGGGGGAAGTGAACTTCATTGCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAg                              >  1:13150598/1‑101 (MQ=255)
                                                                         aaGCGCAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGGCGGGCCAATGGGCGGTATTACTTCTGCCAgc                             >  2:3994269/1‑100 (MQ=255)
                                                                           gcgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTAt                           >  1:1582067/1‑101 (MQ=255)
                                                                            cgcAGGCGCTGGCAATACTGTAGCAGGGGGAAGTGAACTTCATTCCTTGTCCCCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATg                          <  1:7708970/101‑1 (MQ=255)
                                                                            cgcAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATg                          <  1:2424632/101‑1 (MQ=255)
                                                                             gaaGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGc                         >  2:4512935/3‑101 (MQ=255)
                                                                              cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGcc                        <  1:9068466/101‑1 (MQ=255)
                                                                              cAGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGcc                        <  2:3736090/101‑1 (MQ=255)
                                                                               aGGCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCg                       >  1:7432775/1‑101 (MQ=255)
                                                                               aGGCGCTGGCAATACTGGAGCAGTGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCg                       <  2:6521980/101‑1 (MQ=255)
                                                                                tgcgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTGCTGCCAGTATGCCGa                      >  1:5416133/2‑101 (MQ=255)
                                                                                ggCGCTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGa                      >  1:3124303/1‑101 (MQ=255)
                                                                                 nnncTGGCATTACTGGAGCAGAGGGTAGTGGACTTCATTCCTTGTCACCATGTGAATGTCGTCGGGCCAATGGGCGGTATTTCTTCTCCCAGTATGCCGAt                     <  2:13836468/98‑1 (MQ=255)
                                                                                 gcgcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGAt                     <  2:3485755/101‑1 (MQ=255)
                                                                                   gcTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGGATTACTTCTGCCAGTATGCCTATGc                   >  2:2683933/1‑101 (MQ=255)
                                                                                    cTGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGACGGGCCAATGGGCGGGATTACTTCTGCCAGTATGCCGATGCt                  >  1:5863002/1‑101 (MQ=255)
                                                                                     tGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCtg                 <  1:4564988/101‑1 (MQ=255)
                                                                                     tGGCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCtg                 <  1:9946221/101‑1 (MQ=255)
                                                                                      ggCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCtgg                >  1:4899664/1‑101 (MQ=255)
                                                                                      ggCAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCtgg                >  2:3033471/1‑101 (MQ=255)
                                                                                        cAATACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCtggtg              <  1:12757176/101‑1 (MQ=255)
                                                                                         aaTACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCtggtgg             >  1:7264069/1‑101 (MQ=255)
                                                                                           tACTGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGtt           <  2:6786876/101‑1 (MQ=255)
                                                                                           nnntGGAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCtggtggat           <  2:7623869/98‑3 (MQ=255)
                                                                                            nnnggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTg          <  2:11215683/98‑1 (MQ=255)
                                                                                               ggAGCAGGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTgaga       >  2:1551130/1‑101 (MQ=255)
                                                                                                     |                                                                                                     
TGGCAGGAGATGACCGGACCCATGAAAGGGGCGTGCGTGGGCGCATGTCTGTTCGAAGGTTGGGCGAAAGATGAAGCGCAGGCGCTGGCAATACTGGAGCAGGGGAAGTGAACTTCATTCCTTGTCACCATGTGAATGCCGTCGGGCCAATGGGCGGTATTACTTCTGCCAGTATGCCGATGCTGGTGGTTGAGAACGTG  >  NC_000913/547731‑547930

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