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.7% 1184.7 / ‑5.4 367intergenic (+1292/+81)ygcE/ygcFputative kinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (235/129);  total (237/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 378 total aligned reads displayed.

TACCGCAATAATTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAAT  >  NC_000913/2902598‑2902777
                                                                                          |                                                                                         
   cGCAATAATTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGtttt                                                                                      <  2:6867199/93‑1 (MQ=255)
   cGCAATAATTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGtttt                                                                                      >  2:5059714/1‑93 (MQ=255)
    gCAATAATTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTc                                                                                     >  1:6846018/1‑93 (MQ=255)
    gCAATAATTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTc                                                                                     >  1:757627/1‑93 (MQ=255)
    gCAATAATTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTc                                                                                     >  2:2076465/1‑93 (MQ=255)
     cAATAATTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCa                                                                                    <  2:6416211/93‑1 (MQ=255)
     cAATAATTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCa                                                                                    >  1:4250635/1‑93 (MQ=255)
      aataatTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCag                                                                                   >  1:9832179/1‑93 (MQ=255)
      aataatTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCag                                                                                   <  2:9141505/93‑1 (MQ=255)
       ataatTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCaga                                                                                  >  1:4973872/1‑93 (MQ=255)
       ataatTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCaga                                                                                  >  2:2747646/1‑93 (MQ=255)
          atTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGc                                                                               >  1:8102185/1‑93 (MQ=255)
           ttttCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCg                                                                              >  2:5633630/1‑93 (MQ=255)
           ttttCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCg                                                                              >  1:1724791/1‑93 (MQ=255)
           ttttCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCTTGATGGGTAAAATGGAGGGGTTTTCCGAGGCGa                                                                             >  2:1690408/1‑93 (MQ=255)
               cACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGacaa                                                                          >  2:9230196/1‑93 (MQ=255)
               cACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGacaa                                                                          <  2:9612670/93‑1 (MQ=255)
                aCTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGacaaa                                                                         >  1:6975682/1‑93 (MQ=255)
                aCTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGacaaa                                                                         >  1:9295752/1‑93 (MQ=255)
                aCTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGacaaa                                                                         >  1:4893022/1‑93 (MQ=255)
                 cTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGacaaac                                                                        <  1:2709861/93‑1 (MQ=255)
                  tCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGacaaaca                                                                       >  2:5888781/1‑93 (MQ=255)
                   ccAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGacaaacaa                                                                      <  1:1897620/93‑1 (MQ=255)
                   ccAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGacaaacaa                                                                      <  2:5772178/93‑1 (MQ=255)
                    cAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAat                                                                     <  2:1515502/93‑1 (MQ=255)
                     aGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAata                                                                    >  1:1315391/1‑93 (MQ=255)
                     aGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAata                                                                    >  1:8803126/1‑93 (MQ=255)
                     aGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAata                                                                    <  1:1041231/93‑1 (MQ=255)
                      gCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAatat                                                                   >  2:6750685/1‑93 (MQ=255)
                        gAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATaa                                                                 <  1:2460929/93‑1 (MQ=255)
                        gAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGACGGGCTTTCAGAGCAGACAACCAATATaa                                                                 >  1:6513168/1‑93 (MQ=255)
                         aaaaaTTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATaaa                                                                <  1:6002877/93‑1 (MQ=255)
                         aaaaaTTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCCGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATaaa                                                                >  2:3513033/1‑93 (MQ=255)
                            aaTTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGa                                                             >  2:1605049/1‑93 (MQ=255)
                            aaTTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGa                                                             <  1:7840968/93‑1 (MQ=255)
                             aTTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAg                                                            >  1:2084009/1‑93 (MQ=255)
                              ttAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGATTTCAGAGGCGACAAACAATATAAATGAGt                                                           >  2:1941813/1‑93 (MQ=255)
                              tgaaTGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGt                                                           <  1:5057027/91‑1 (MQ=255)
                               tAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTa                                                          >  2:6526495/1‑93 (MQ=255)
                               tAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTa                                                          <  2:7067393/93‑1 (MQ=255)
                                aaTGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAg                                                         >  1:5730071/1‑93 (MQ=255)
                                  tGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGAGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGaa                                                       <  1:5231923/93‑1 (MQ=255)
                                      acaGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCt                                                   >  1:9016974/1‑93 (MQ=255)
                                         gAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTa                                                >  1:831641/1‑93 (MQ=255)
                                          aaTTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTAc                                               <  1:4757893/93‑1 (MQ=255)
                                            tttGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAg                                             >  1:5685243/1‑93 (MQ=255)
                                            tttGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAg                                             >  2:22385/1‑93 (MQ=255)
                                            attGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAg                                             <  2:8197516/92‑1 (MQ=255)
                                              tGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGt                                           <  2:600413/93‑1 (MQ=255)
                                               gTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGtt                                          >  2:7964783/1‑93 (MQ=255)
                                               gTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGtt                                          >  2:7951019/1‑93 (MQ=255)
                                                tAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTc                                         >  1:1001753/1‑93 (MQ=255)
                                                tAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTc                                         >  1:3347565/1‑93 (MQ=255)
                                                tAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTc                                         >  1:6380426/1‑93 (MQ=255)
                                                   agaATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGGGTAGAAGTCTTTACAGGTTCCAc                                      >  2:8816736/1‑93 (MQ=255)
                                                     aaTGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCAccc                                    <  1:1618663/93‑1 (MQ=255)
                                                      aTGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCa                                   >  1:8558481/1‑93 (MQ=255)
                                                       tGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCAt                                  >  1:9760455/1‑93 (MQ=255)
                                                       tGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCAt                                  >  2:8682600/1‑93 (MQ=255)
                                                          tAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTta                               >  1:4018923/1‑93 (MQ=255)
                                                          tAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTta                               >  2:4249238/1‑93 (MQ=255)
                                                          tAATGATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTta                               >  2:6248739/1‑93 (MQ=255)
                                                              gATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAAc                           <  2:2528007/93‑1 (MQ=255)
                                                              gATTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAAc                           <  2:1178902/93‑1 (MQ=255)
                                                               aTTAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACg                          >  1:8392664/1‑93 (MQ=255)
                                                                ttAGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGc                         >  1:7095871/1‑93 (MQ=255)
                                                                  aGCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCtt                       >  2:3932563/1‑93 (MQ=255)
                                                                   gCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTta                      <  2:8494362/93‑1 (MQ=255)
                                                                   gCCCTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTta                      >  1:9473360/1‑93 (MQ=255)
                                                                     ccTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACGCATTTATAACGCTtata                    <  2:3223756/93‑1 (MQ=255)
                                                                      cTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATaa                   >  2:9520931/1‑93 (MQ=255)
                                                                      cTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATaa                   >  2:9442397/1‑93 (MQ=255)
                                                                      cTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATaa                   >  1:5908245/1‑93 (MQ=255)
                                                                      cTGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATaa                   >  1:6229439/1‑93 (MQ=255)
                                                                      cTGGATGGGTAAAATGGAGGGGATTTCAGCGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACACATTTATAACGCTTATaa                   >  2:9764671/1‑93 (MQ=255)
                                                                       tGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATaaa                  >  1:3913371/1‑93 (MQ=255)
                                                                       tGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATaaa                  <  1:9662983/93‑1 (MQ=255)
                                                                       tGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATaaa                  >  1:5171756/1‑93 (MQ=255)
                                                                       tGGATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATaaa                  >  2:6657872/1‑93 (MQ=255)
                                                                        ggATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAAGGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAAt                 <  2:5446896/93‑1 (MQ=255)
                                                                         gATGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATg                <  2:3412031/93‑1 (MQ=255)
                                                                           tGGGTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGtt              <  1:324760/93‑1 (MQ=255)
                                                                            gggTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGttt             <  1:7407849/93‑1 (MQ=255)
                                                                             ggTAAAGTGGAGAGGTTTTCAGAGGCGACAAACTATATATATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGtttc            <  2:8984631/93‑2 (MQ=255)
                                                                             ggTAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTa            >  2:3072527/1‑93 (MQ=255)
                                                                               tAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAAt          >  2:5267278/1‑93 (MQ=255)
                                                                               tAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAAt          <  2:7798988/93‑1 (MQ=255)
                                                                               tAAAATGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAAt          >  1:7371836/1‑93 (MQ=255)
                                                                                  aaTGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAg       >  2:4658138/1‑93 (MQ=255)
                                                                                  aaTGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAg       <  2:4119985/93‑1 (MQ=255)
                                                                                   aTGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAgg      >  2:5284569/1‑93 (MQ=255)
                                                                                   aTGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAgg      >  1:3625859/1‑93 (MQ=255)
                                                                                    tGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAgga     >  2:2608607/1‑92 (MQ=255)
                                                                                    tGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAGGc     >  1:9152224/1‑93 (MQ=255)
                                                                                    tGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAGGc     <  1:5444953/93‑1 (MQ=255)
                                                                                    tGGAGGGGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAGGc     <  1:6126492/93‑1 (MQ=255)
                                                                                    aggagggGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAGGc     <  1:2221086/92‑1 (MQ=255)
                                                                                      gaggggTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAGGCaa   >  1:4367157/1‑93 (MQ=255)
                                                                                      gaggggTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAGGCaa   <  1:2741241/93‑1 (MQ=255)
                                                                                      gaggggGTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAGGCaa   >  2:4778380/1‑93 (MQ=255)
                                                                                          |                                                                                         
TACCGCAATAATTTTCACTCCAGCGAAAAATTAATGCCACAGAATTTGTAGAGAATGCTAATGATTAGCCCTGGATGGGTAAAATGGAGGAGTTTTCAGAGGCGACAAACAATATAAATGAGTAGAAGTCTTTACAGGTTCCACCCATTTATAACGCTTATAAATGTTTAATCAGGCAAT  >  NC_000913/2902598‑2902777

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