Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,510,121 +G intergenic (‑358/+151) yhfZ ← / ← trpS hypothetical protein/tryptophan tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,510,1201.G95.2% 507.6 / 17.1 167intergenic (‑357/+152)yhfZ/trpShypothetical protein/tryptophan tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base . (7/1);  new base G (100/59);  total (107/60)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 191 total aligned reads displayed.

AAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAG  >  NC_000913/3510028‑3510208
                                                                                                |                                                                                        
aaGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTcccc                                                                                           >  1:3742307/1‑93 (MQ=255)
aaGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTcccc                                                                                           <  2:617776/93‑1 (MQ=255)
 aGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAACAAACCGGGTTTCCCCg                                                                                         >  2:4659420/1‑93 (MQ=255)
 aGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCg                                                                                         >  1:1445649/1‑93 (MQ=255)
  gATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCgg                                                                                        >  1:1631916/1‑92 (MQ=255)
  gATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCgg                                                                                        >  2:4043287/1‑92 (MQ=255)
   aTAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCggt                                                                                       <  1:4228466/93‑3 (MQ=255)
    tAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGgtt                                                                                      >  2:1333739/1‑93 (MQ=255)
      aGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGgtttg                                                                                    >  1:2852224/1‑91 (MQ=255)
      aGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGgtttc                                                                                    >  1:4613694/1‑91 (MQ=255)
      aGGAGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGgtttc                                                                                    >  2:1257017/1‑91 (MQ=255)
          gcAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCagag                                                                                 <  1:4321652/93‑1 (MQ=255)
           cAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGt                                                                                <  2:435376/93‑1 (MQ=255)
            gaCAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtg                                                                               >  1:4422301/2‑93 (MQ=255)
            aaCAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtg                                                                               >  1:3429429/1‑93 (MQ=255)
             aCAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCTCCGGTTTCAGAGtga                                                                              >  1:3287407/1‑93 (MQ=255)
             aCAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtga                                                                              <  1:1278796/93‑1 (MQ=255)
             aCAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGgga                                                                              >  2:416836/1‑93 (MQ=255)
              cAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtgat                                                                             >  1:1411761/1‑93 (MQ=255)
                gCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGtgatga                                                                           <  1:1877460/93‑1 (MQ=255)
                  gagaTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATa                                                                         <  2:3814235/93‑1 (MQ=255)
                   agaTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATaa                                                                        <  2:2673514/93‑1 (MQ=255)
                    gaTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATaaa                                                                       <  1:4047184/93‑1 (MQ=255)
                     aTGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATaaaa                                                                      <  2:3549905/93‑1 (MQ=255)
                        tgGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCa                                                                   <  1:774831/93‑1 (MQ=255)
                         ggAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCaa                                                                  <  2:2927728/93‑1 (MQ=255)
                          gAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGATTCAGAGTGATGATAAAAGCaaa                                                                 >  2:4027123/1‑93 (MQ=255)
                            aaCGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAAt                                                               <  2:4254318/93‑1 (MQ=255)
                             aCGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCCGTTTCAGAGTGATGATAAAAGCAAAAtt                                                              >  2:3671780/1‑93 (MQ=255)
                             aCGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAAtt                                                              <  2:794854/93‑1 (MQ=255)
                                  gAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTg                                                         >  2:2558014/1‑93 (MQ=255)
                                  gAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTg                                                         <  2:2524285/93‑1 (MQ=255)
                                  gAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTg                                                         >  2:2471470/1‑93 (MQ=255)
                                      cACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgc                                                     <  2:4415750/93‑1 (MQ=255)
                                        cTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCAGGTTTCAGAGTGATGATAAAAGCAAAGTTGCCTGATgcgc                                                   >  2:2002918/1‑93 (MQ=255)
                                        cTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgcgc                                                   >  2:3239490/1‑93 (MQ=255)
                                        cTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATgcgc                                                   >  1:541905/1‑93 (MQ=255)
                                         ttCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCt                                                  >  1:2508329/1‑93 (MQ=255)
                                         ttCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCt                                                  >  1:747369/1‑93 (MQ=255)
                                         ttCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCt                                                  >  1:3895652/1‑93 (MQ=255)
                                          tCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTa                                                 >  1:1790488/1‑93 (MQ=255)
                                          tCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTa                                                 >  2:278014/1‑93 (MQ=255)
                                           cacaTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTAc                                                >  2:68842/1‑93 (MQ=255)
                                           cacaTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTAc                                                >  2:2117913/1‑93 (MQ=255)
                                           cacaTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTAc                                                >  2:1576076/1‑93 (MQ=255)
                                            acaTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACg                                               >  1:106809/1‑93 (MQ=255)
                                            acaTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACg                                               >  2:4683950/1‑93 (MQ=255)
                                             caTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGc                                              <  1:3552897/93‑1 (MQ=255)
                                             caTTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGc                                              >  2:44125/1‑93 (MQ=255)
                                              aTTTTTTCGGATTATTGCAAATGGGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCt                                             >  2:3779875/1‑93 (MQ=255)
                                               ttttttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCtt                                            >  2:1957525/1‑93 (MQ=255)
                                                  tttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATc                                         <  1:20596/93‑1 (MQ=255)
                                                   ttCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCa                                        >  1:2143010/1‑93 (MQ=255)
                                                    tCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAg                                       <  1:4247933/93‑1 (MQ=255)
                                                    tCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAAATGCCTGATGCGCTACGCTTATCAg                                       <  1:2667971/93‑1 (MQ=255)
                                                     cGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAgg                                      >  1:1951852/1‑93 (MQ=255)
                                                     cGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAgg                                      >  1:1081741/1‑93 (MQ=255)
                                                     cGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAgg                                      >  2:5951/1‑93 (MQ=255)
                                                      ggATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGc                                     <  1:3906397/93‑1 (MQ=255)
                                                      ggATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGGTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGc                                     >  2:1428463/1‑93 (MQ=255)
                                                        attattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGACt                                   >  2:361344/1‑93 (MQ=255)
                                                         ttattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTa                                  >  1:3146185/1‑93 (MQ=255)
                                                          tattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTAc                                 <  2:1031564/93‑1 (MQ=255)
                                                          tattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTAc                                 >  1:3806934/1‑93 (MQ=255)
                                                          tattGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTAc                                 >  1:610204/1‑93 (MQ=255)
                                                           attGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACa                                <  1:2287327/93‑1 (MQ=255)
                                                           attGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACa                                <  2:3389311/93‑1 (MQ=255)
                                                            ttGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAt                               >  1:4584983/1‑93 (MQ=255)
                                                             tGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACAtt                              >  1:3413983/1‑93 (MQ=255)
                                                                aaaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTcc                           >  2:4341269/1‑93 (MQ=255)
                                                                aaaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTcc                           >  1:3716592/1‑93 (MQ=255)
                                                                 aaTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCt                          >  2:3049669/1‑93 (MQ=255)
                                                                  aTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCtt                         >  1:2365320/1‑93 (MQ=255)
                                                                  aTGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCtt                         <  2:4673556/93‑1 (MQ=255)
                                                                      ggATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCaa                     >  2:86089/1‑93 (MQ=255)
                                                                      ggATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCaa                     >  1:3045298/1‑93 (MQ=255)
                                                                        aTAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAata                   >  2:3753189/1‑93 (MQ=255)
                                                                         tAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAatat                  >  1:1230396/1‑93 (MQ=255)
                                                                            aaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATAtgtg               >  1:4483320/1‑93 (MQ=255)
                                                                              aaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCa             >  2:1805349/1‑93 (MQ=255)
                                                                               aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCGTTGCAATATGTGCAt            >  2:3277061/1‑93 (MQ=255)
                                                                               aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt            <  2:1148756/93‑1 (MQ=255)
                                                                               aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt            <  1:931127/93‑1 (MQ=255)
                                                                               aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt            >  1:2191597/1‑93 (MQ=255)
                                                                               aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt            <  1:4287201/93‑1 (MQ=255)
                                                                               aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt            >  2:525635/1‑93 (MQ=255)
                                                                               aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt            >  2:2094649/1‑93 (MQ=255)
                                                                                aCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtt           >  1:2177271/1‑93 (MQ=255)
                                                                                  cGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTAc         >  2:1055068/1‑93 (MQ=255)
                                                                                  cGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTAc         <  1:811268/93‑1 (MQ=255)
                                                                                  cGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTAc         <  1:27738/93‑1 (MQ=255)
                                                                                  cGGGTTTCCCCGGTTTCAGAGTGAGGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTAc         <  2:1190789/93‑1 (MQ=255)
                                                                                    ggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCCTGGCTACATTTCCTTGCAATATGTGCATTACtt       >  1:4036682/1‑93 (MQ=37)
                                                                                    ggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACtt       >  2:4084694/1‑93 (MQ=255)
                                                                                    ggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACtt       <  2:1463836/93‑1 (MQ=255)
                                                                                    ggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACtt       >  2:1859811/1‑93 (MQ=255)
                                                                                      tttCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTCTg     >  2:389285/1‑93 (MQ=255)
                                                                                        tccccgGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTa   <  2:4617614/91‑1 (MQ=255)
                                                                                        tccccgGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTa   <  2:3371697/91‑1 (MQ=255)
                                                                                         ccccgGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAg  <  2:1599537/92‑1 (MQ=255)
                                                                                                |                                                                                        
AAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAG  >  NC_000913/3510028‑3510208

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: