Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,718,305 T→C intergenic (‑76/‑4) yiaZ ← / → yi5A putative membrane permeability altering protein/IS150 hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,718,3050TC84.2% 483.5 / 84.7 196intergenic (‑76/‑4)yiaZ/yi5Aputative membrane permeability altering protein/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (16/15);  new base C (77/88);  total (93/103)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

Only 100 of 203 total aligned reads displayed.

ATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTC  >  NC_000913/3718215‑3718394
                                                                                          |                                                                                         
atatTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCaa                                                                                         >  2:2376361/1‑93 (MQ=255)
 tatTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAat                                                                                        >  2:549325/1‑93 (MQ=255)
  atTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAata                                                                                       >  1:436055/1‑93 (MQ=255)
  atTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAata                                                                                       >  1:2507021/1‑93 (MQ=255)
   tttCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAatat                                                                                      >  1:3633408/1‑93 (MQ=255)
    ttCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATg                                                                                     <  2:2639279/93‑1 (MQ=255)
     tCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGt                                                                                    >  2:2110345/1‑93 (MQ=255)
     tCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGt                                                                                    <  1:2662274/93‑1 (MQ=255)
      cTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTc                                                                                   <  1:2403051/93‑1 (MQ=255)
      cTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTc                                                                                   <  2:5423842/93‑1 (MQ=255)
       tGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCa                                                                                  >  1:5252978/1‑93 (MQ=255)
       tGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCa                                                                                  >  2:4409026/1‑93 (MQ=255)
        gCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCaa                                                                                 >  1:335907/1‑93 (MQ=255)
        gCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCaa                                                                                 >  2:3452114/1‑93 (MQ=255)
        gCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCaa                                                                                 >  2:5766571/1‑93 (MQ=255)
         cGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCaaa                                                                                <  1:1256887/93‑1 (MQ=255)
         cGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCaaa                                                                                <  2:6102223/93‑1 (MQ=255)
         cGGCATTCCAAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCaaa                                                                                <  1:3597400/93‑1 (MQ=255)
          ggCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAg                                                                               <  1:2226061/93‑1 (MQ=255)
           gCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGc                                                                              >  2:840837/1‑93 (MQ=255)
            cATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGcc                                                                             <  2:3531379/93‑1 (MQ=255)
              ttCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCaa                                                                           >  2:1225428/1‑93 (MQ=255)
               tCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCaaa                                                                          <  1:207557/93‑1 (MQ=255)
                cccAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCaaaa                                                                         >  1:1899650/1‑93 (MQ=255)
                 ccAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAAt                                                                        >  2:5855223/1‑93 (MQ=255)
                  cAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATa                                                                       >  2:5950639/1‑93 (MQ=255)
                  cAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATa                                                                       >  1:3468473/1‑93 (MQ=255)
                  cAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATa                                                                       >  1:1462755/1‑93 (MQ=255)
                   aGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATAc                                                                      >  2:4641958/1‑93 (MQ=255)
                   aGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATAc                                                                      >  2:4274064/1‑93 (MQ=255)
                    gCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATAcc                                                                     <  2:244458/93‑1 (MQ=255)
                    gCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATAcc                                                                     <  1:6199312/93‑1 (MQ=255)
                    gCCTCCTGAGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATAcc                                                                     >  1:6502536/1‑93 (MQ=255)
                     cctcctGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATAccc                                                                    <  1:3277568/93‑1 (MQ=255)
                      ctcctGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCt                                                                   >  2:877060/1‑93 (MQ=255)
                      ctcctGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCt                                                                   >  1:3778250/1‑93 (MQ=255)
                      ctcctGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCt                                                                   <  1:2306318/93‑1 (MQ=255)
                       tcctGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCtt                                                                  >  1:5026623/1‑93 (MQ=255)
                       tcctGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCtt                                                                  <  2:2855479/93‑1 (MQ=255)
                       tcctGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCtt                                                                  >  1:756966/1‑93 (MQ=255)
                         ctGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCtttt                                                                >  1:4126210/1‑93 (MQ=255)
                         ctGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCtttt                                                                >  1:3531080/1‑93 (MQ=255)
                         ctGTGTTGATTTCCAAAGAGTATGTACTGCACCCATTTGGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCtttt                                                                <  2:4838872/93‑1 (MQ=255)
                          tgtgTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTg                                                               >  2:5002554/1‑93 (MQ=255)
                          cgtgtTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTg                                                               <  2:4314904/92‑1 (MQ=255)
                           gtgtTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGa                                                              >  1:5676826/1‑93 (MQ=255)
                             gtTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGaaa                                                            >  1:5071432/1‑93 (MQ=255)
                              ttGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGaaaa                                                           >  2:4581855/1‑93 (MQ=255)
                              ttGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGaaaa                                                           >  1:3129963/1‑93 (MQ=255)
                                gATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAgc                                                         >  1:5836976/1‑93 (MQ=255)
                                gATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAgc                                                         <  1:5960777/93‑1 (MQ=255)
                                gATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAgc                                                         >  1:5306360/1‑93 (MQ=255)
                                  tttCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAgcgc                                                       >  1:291761/1‑93 (MQ=255)
                                   ttCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGcc                                                      <  2:1451205/93‑1 (MQ=255)
                                    tCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCt                                                     <  2:6634223/93‑1 (MQ=255)
                                       aaCGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGa                                                  >  1:5992372/1‑93 (MQ=255)
                                       aaCGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGa                                                  >  2:5977732/1‑93 (MQ=255)
                                          gAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGt                                               <  1:3157967/93‑1 (MQ=255)
                                          gAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGt                                               <  1:5389760/93‑1 (MQ=255)
                                                 tACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAAt                                        <  1:3076749/93‑1 (MQ=255)
                                                  aCTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATc                                       <  1:1384915/93‑1 (MQ=255)
                                                   cTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCa                                      <  2:3043489/93‑1 (MQ=255)
                                                   cTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCa                                      <  2:1530958/93‑1 (MQ=255)
                                                    tGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCac                                     <  2:4366025/93‑1 (MQ=255)
                                                     gCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCact                                    <  1:4952618/93‑1 (MQ=255)
                                                     gCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCact                                    <  2:5239142/93‑1 (MQ=255)
                                                       aCCCATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCactac                                  <  1:943582/93‑1 (MQ=255)
                                                         ccATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACtt                                <  1:4332632/93‑1 (MQ=255)
                                                         ccATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACtt                                >  1:853247/1‑93 (MQ=255)
                                                         ccATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACtt                                <  1:641584/93‑1 (MQ=255)
                                                          cATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTc                               >  1:1004365/1‑93 (MQ=255)
                                                          cATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTc                               >  1:6241719/1‑93 (MQ=255)
                                                          cATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTc                               <  2:3290220/93‑1 (MQ=255)
                                                          cATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTc                               <  2:2697014/93‑1 (MQ=255)
                                                          cATTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGCCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTc                               >  2:4790848/1‑93 (MQ=255)
                                                           aTTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTca                              >  1:4611329/1‑93 (MQ=255)
                                                           aTTTTGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTca                              <  2:3758040/93‑1 (MQ=255)
                                                             tttGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATCCCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTcaca                            <  2:4092300/93‑1 (MQ=255)
                                                             tttGTTGGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATTCCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTcaca                            <  2:2953886/93‑1 (MQ=255)
                                                              ttgttgGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACaa                           >  1:4861313/1‑93 (MQ=255)
                                                                gttgGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACt                         >  1:5741067/1‑93 (MQ=255)
                                                                  tgGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACtga                       <  2:4126598/93‑1 (MQ=255)
                                                                  tgGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACtga                       >  2:646716/1‑93 (MQ=255)
                                                                  tgGACGATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACtga                       <  1:388480/93‑1 (MQ=255)
                                                                       gATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGatgatg                  <  2:3901411/93‑1 (MQ=255)
                                                                       gATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGatgatg                  >  2:1976779/1‑93 (MQ=255)
                                                                       gATGAAATGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGatgatg                  <  2:3869561/93‑1 (MQ=255)
                                                                           aaaTGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTa              >  1:3676824/1‑93 (MQ=255)
                                                                            aaTGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTAc             >  2:5112847/1‑93 (MQ=255)
                                                                             aTGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACa            >  1:2465265/1‑93 (MQ=255)
                                                                             aTGGAATAGCCCCCAATATGTCAAAGCCAAAACACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACa            <  2:714520/93‑1 (MQ=255)
                                                                              tGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAg           <  1:6311967/93‑1 (MQ=255)
                                                                              tGGAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAg           >  1:3858922/1‑93 (MQ=255)
                                                                               ggAATAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAgg          <  1:5510851/93‑1 (MQ=255)
                                                                                 aaTAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGat        >  1:4007486/1‑93 (MQ=255)
                                                                                 aaTAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGat        <  2:4326293/93‑1 (MQ=255)
                                                                                 aaTAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGat        >  2:813416/1‑93 (MQ=255)
                                                                                   tAGCCCCCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGatca      <  1:4322032/93‑1 (MQ=255)
                                                                                      ccccCAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCAtct   <  1:4742104/93‑1 (MQ=255)
                                                                                       ccccAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCAtctc  <  1:4665686/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATATTTCTGCGGCATTCCCAGCCTCCTGTGTTGATTTCCAACGAGTATGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTC  >  NC_000913/3718215‑3718394

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