Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,902,688 A→G intergenic (+1292/+81) ygcE → / ← ygcF putative kinase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,902,6880AG99.2% 363.4 / ‑4.6 117intergenic (+1292/+81)ygcE/ygcFputative kinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (92/24);  total (93/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

Only 100 of 120 total aligned reads displayed.

TACCGCAATAATTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAGG  >  NC_000913/2902598‑2902773
                                                                                          |                                                                                     
tACCGCAATAATTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGt                                                                                     >  2:2401301/1‑93 (MQ=255)
   cGCAATAATTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGtttt                                                                                  <  1:4572896/93‑1 (MQ=255)
    gCAATAATTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTc                                                                                 >  1:1368987/1‑93 (MQ=255)
      aataatTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCag                                                                               >  2:961656/1‑93 (MQ=255)
       ataatTTTCACTCCAGCGATAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCaga                                                                              >  2:4567323/1‑93 (MQ=255)
       ataatTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCaga                                                                              >  2:2914925/1‑93 (MQ=255)
       ataatTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCaga                                                                              >  1:4195688/1‑93 (MQ=255)
        taatTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCagag                                                                             >  1:4531365/1‑93 (MQ=255)
        taatTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCagag                                                                             >  1:2602549/1‑93 (MQ=255)
        taatTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCagag                                                                             >  1:4263437/1‑93 (MQ=255)
        taatTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCagag                                                                             >  1:4489999/1‑93 (MQ=255)
          atTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGc                                                                           >  1:4023557/1‑93 (MQ=255)
             ttCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGac                                                                        >  2:741091/1‑93 (MQ=255)
              tCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGaca                                                                       <  2:205905/93‑1 (MQ=255)
              tCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGaca                                                                       >  2:2536277/1‑93 (MQ=255)
              tCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGaca                                                                       >  2:3047432/1‑93 (MQ=255)
                aCTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGacaaa                                                                     <  2:1288858/93‑1 (MQ=255)
                aCTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGacaaa                                                                     >  1:25406/1‑93 (MQ=255)
                aCTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGacaaa                                                                     >  1:3749090/1‑93 (MQ=255)
                   ccAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGacaaacaa                                                                  >  2:3649040/1‑93 (MQ=255)
                     aGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAata                                                                >  1:4282237/1‑93 (MQ=255)
                     aGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAata                                                                >  1:1826209/1‑93 (MQ=255)
                      gCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAatat                                                               >  1:527902/1‑93 (MQ=255)
                      gCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAatat                                                               <  1:3363844/93‑1 (MQ=255)
                      gCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAatat                                                               >  1:2695689/1‑93 (MQ=255)
                       cGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAAtata                                                              <  2:2358808/93‑1 (MQ=255)
                         aaaaaTTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATaaa                                                            >  1:1369528/1‑93 (MQ=255)
                          aaaaTTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAAt                                                           >  2:1258602/1‑93 (MQ=255)
                            aaTTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGa                                                         >  2:1797798/1‑93 (MQ=255)
                               tAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTa                                                      >  1:1569643/1‑93 (MQ=255)
                                  tGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGaa                                                   >  1:2084241/1‑93 (MQ=255)
                                    ccACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGGCAAAAAAAATAAATAAGTAAAAga                                                 >  2:3364969/1‑92 (MQ=255)
                                     cacaGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTc                                                >  1:3113138/1‑93 (MQ=255)
                                     cacaGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTc                                                >  2:71629/1‑93 (MQ=255)
                                     cacaGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGCTTTCAGAGGCGACAAACAATATAAATGAGTAGAAggc                                                >  1:3785009/1‑91 (MQ=255)
                                     cacaGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGGAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTc                                                >  1:1795629/1‑93 (MQ=255)
                                       caGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCtt                                              >  1:396902/1‑93 (MQ=255)
                                       caGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCtt                                              >  2:822499/1‑93 (MQ=255)
                                       caGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCtt                                              >  2:3220210/1‑93 (MQ=255)
                                       caGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGACtt                                              >  2:2584447/1‑93 (MQ=255)
                                        aGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCttt                                             >  2:4412734/1‑93 (MQ=255)
                                        aGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCttt                                             <  1:2583586/93‑1 (MQ=255)
                                         gAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTa                                            >  2:3058397/1‑93 (MQ=255)
                                         gAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTa                                            >  1:1118295/1‑93 (MQ=255)
                                         gAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTa                                            >  2:4471617/1‑93 (MQ=255)
                                         gAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTa                                            <  2:4387327/93‑1 (MQ=255)
                                           aTTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACa                                          >  1:3170513/1‑93 (MQ=255)
                                           aTTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACa                                          >  1:4339293/1‑93 (MQ=255)
                                              tGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGt                                       >  1:201710/1‑93 (MQ=255)
                                               gTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGtt                                      >  1:2865834/1‑93 (MQ=255)
                                                tAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTc                                     <  2:2359871/93‑1 (MQ=255)
                                                 agagAATGCTAATGATTAGCCCTGGCTGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTcc                                    >  2:3437749/1‑93 (MQ=255)
                                                  gagaATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCa                                   >  2:3341951/1‑93 (MQ=255)
                                                  gagaATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCa                                   >  2:1233816/1‑93 (MQ=255)
                                                    gaATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGGTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGGTCCAcc                                 >  2:232769/1‑93 (MQ=255)
                                                      aTGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCa                               <  2:3690831/93‑1 (MQ=255)
                                                      aTGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCa                               <  2:3754634/93‑1 (MQ=255)
                                                       tGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCAt                              >  1:1924694/1‑93 (MQ=255)
                                                        gCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCAtt                             >  1:3157845/1‑93 (MQ=255)
                                                          tAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACCAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTta                           >  2:2358982/1‑93 (MQ=255)
                                                          tAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTta                           >  1:667373/1‑93 (MQ=255)
                                                          tAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTta                           <  1:3800607/93‑1 (MQ=255)
                                                          tAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTta                           >  2:4594998/1‑93 (MQ=255)
                                                           aaTGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTtat                          >  1:3557098/1‑93 (MQ=255)
                                                           aaTGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTtat                          >  1:2528339/1‑93 (MQ=255)
                                                            aTGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTtata                         >  1:673872/1‑93 (MQ=255)
                                                            aTGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTtata                         >  1:884717/1‑93 (MQ=255)
                                                             tGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATaa                        >  1:1829657/1‑93 (MQ=255)
                                                             tGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATaa                        >  2:2693090/1‑93 (MQ=255)
                                                             tGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATaa                        >  2:767100/1‑93 (MQ=255)
                                                              gATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAAc                       >  1:3264589/1‑93 (MQ=255)
                                                              gATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAAc                       >  1:997920/1‑93 (MQ=255)
                                                              gATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAAc                       <  2:458543/93‑1 (MQ=255)
                                                               aTTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACg                      >  1:2224234/1‑93 (MQ=255)
                                                                ttAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGc                     <  2:2234197/93‑1 (MQ=255)
                                                                ttAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGc                     >  2:2094460/1‑93 (MQ=255)
                                                                 tAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCt                    >  1:1372266/1‑93 (MQ=255)
                                                                    cccTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTtat                 <  1:3444734/93‑1 (MQ=255)
                                                                    cccTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTtat                 >  1:3045432/1‑93 (MQ=255)
                                                                      cTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATaa               >  1:4369117/1‑93 (MQ=255)
                                                                      cTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATaa               >  1:3589205/1‑93 (MQ=255)
                                                                       tGGATGGGTAAAATGTAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATaaa              >  2:3867773/1‑93 (MQ=255)
                                                                       tGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATaaa              >  2:2783068/1‑93 (MQ=255)
                                                                       tGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATaaa              <  2:1396471/93‑1 (MQ=255)
                                                                        ggATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAAt             >  2:1700674/1‑93 (MQ=255)
                                                                        ggATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAAt             >  2:4512097/1‑93 (MQ=255)
                                                                         gATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATg            >  2:44919/1‑93 (MQ=255)
                                                                          aTGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGt           <  1:2399658/93‑1 (MQ=255)
                                                                           tGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGtt          >  2:2155125/1‑93 (MQ=255)
                                                                            gggTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGttt         >  2:3319448/1‑93 (MQ=255)
                                                                             ggTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTa        >  1:1730531/1‑93 (MQ=255)
                                                                             ggTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTa        >  1:1791925/1‑93 (MQ=255)
                                                                             ggTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTa        >  2:3439642/1‑93 (MQ=255)
                                                                               tAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAAt      <  1:4392295/93‑1 (MQ=255)
                                                                                 aaaTGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCa    <  2:1516087/93‑1 (MQ=255)
                                                                                  aaTGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAg   <  2:4158230/93‑1 (MQ=255)
                                                                                  aaTGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAg   >  2:324677/1‑93 (MQ=255)
                                                                                  aaTGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAg   >  2:316617/1‑93 (MQ=255)
                                                                                  aaTGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAg   >  1:3961456/1‑93 (MQ=255)
                                                                                   aTGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAgg  >  2:3291720/1‑93 (MQ=255)
                                                                                          |                                                                                     
TACCGCAATAATTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAGG  >  NC_000913/2902598‑2902773

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: