Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 502,653 C→T intergenic (‑191/+47) ybaL ← / ← fsr putative transport protein/fosmidomycin resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913502,6530CT99.5% 765.9 / NA 209intergenic (‑191/+47)ybaL/fsrputative transport protein/fosmidomycin resistance protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (103/105);  total (104/105)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

ATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATC  >  NC_000913/502563‑502742
                                                                                          |                                                                                         
attattGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGa                                                                                         >  2:5153533/1‑93 (MQ=255)
  tattGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACt                                                                                       >  2:8333/1‑93 (MQ=255)
  tattGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACt                                                                                       >  2:3089026/1‑93 (MQ=255)
  tattGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACt                                                                                       >  2:3666375/1‑93 (MQ=255)
    ttGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGt                                                                                     <  1:4927760/93‑1 (MQ=255)
    ttGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGt                                                                                     >  1:1508495/1‑93 (MQ=255)
     tGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGACGAAAGGCTTATGGTGACTGTa                                                                                    <  1:4012893/93‑1 (MQ=255)
       aaTTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaaa                                                                                  >  2:433165/1‑93 (MQ=255)
        aTTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaaaa                                                                                 >  2:354652/1‑93 (MQ=255)
         ttttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTATAAAAc                                                                                <  1:4304146/93‑1 (MQ=255)
          tttGCAGAAAAATGGAATTAGCTGCGAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACg                                                                               >  1:237333/1‑93 (MQ=255)
           ttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGa                                                                              >  1:2293092/1‑93 (MQ=255)
           ttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTCTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGa                                                                              >  2:2951045/1‑93 (MQ=255)
                gAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAg                                                                         >  1:4650127/1‑93 (MQ=255)
                  aaaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGcc                                                                       >  2:5776156/1‑93 (MQ=255)
                  aaaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGcc                                                                       <  2:946/93‑1 (MQ=255)
                   aaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCa                                                                      >  1:1175337/1‑93 (MQ=255)
                   aaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCa                                                                      >  2:4792189/1‑93 (MQ=255)
                   aaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCa                                                                      >  1:1407978/1‑93 (MQ=255)
                    aaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCaa                                                                     >  2:505352/1‑93 (MQ=255)
                    aaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCaa                                                                     >  2:2583087/1‑93 (MQ=255)
                     aTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCaaa                                                                    >  1:3130327/1‑93 (MQ=255)
                     aTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCaaa                                                                    >  1:3008115/1‑93 (MQ=255)
                      tGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAg                                                                   >  2:1652679/1‑93 (MQ=255)
                        gAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGtt                                                                 <  2:390481/93‑1 (MQ=255)
                        gAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGtt                                                                 >  1:4735790/1‑93 (MQ=255)
                          aTTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGTAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTg                                                               <  1:1521983/93‑1 (MQ=255)
                          aTTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTg                                                               >  1:2235804/1‑93 (MQ=255)
                          aTTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTg                                                               >  1:5570310/1‑93 (MQ=255)
                              gCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCtt                                                           <  2:1976518/93‑1 (MQ=255)
                              gCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCtt                                                           <  1:708596/93‑1 (MQ=255)
                              gCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCtt                                                           >  1:4207821/1‑93 (MQ=255)
                               cTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCttt                                                          <  2:4128823/93‑1 (MQ=255)
                                tGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAGAACGATGCAGACAAAGTTTGCCTTTg                                                         >  1:3271911/1‑93 (MQ=255)
                                  gcaaaaaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGAAGACAAAGTTTGCCTTTGGc                                                       >  2:5221495/3‑93 (MQ=255)
                                  cAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGc                                                       >  1:3857066/1‑93 (MQ=255)
                                  cAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGc                                                       <  1:5278526/93‑1 (MQ=255)
                                  cAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGc                                                       <  1:4258448/93‑1 (MQ=255)
                                  cAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGc                                                       <  2:502316/93‑1 (MQ=255)
                                    aaaaaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTg                                                     >  1:6119523/1‑93 (MQ=255)
                                      aaaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCa                                                   <  2:2003103/93‑1 (MQ=255)
                                       aaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCaa                                                  <  1:3407598/93‑1 (MQ=255)
                                        aaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAAt                                                 <  1:270207/93‑1 (MQ=255)
                                        aaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAAt                                                 >  1:3905845/1‑93 (MQ=255)
                                          gCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGa                                               <  1:6070254/93‑1 (MQ=255)
                                          gCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGa                                               <  1:897291/93‑1 (MQ=255)
                                           cACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGaa                                              <  1:635079/93‑1 (MQ=255)
                                           cACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGaa                                              >  2:3250815/1‑93 (MQ=255)
                                           cACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGaa                                              >  2:1385252/1‑93 (MQ=255)
                                             cGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAAt                                            <  1:2886197/93‑1 (MQ=255)
                                              ggATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATc                                           <  1:2819302/93‑1 (MQ=255)
                                               gATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCa                                          <  1:788972/93‑1 (MQ=255)
                                                 tCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGt                                        >  2:3121072/1‑93 (MQ=255)
                                                 tCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGt                                        >  1:3719854/1‑93 (MQ=255)
                                                   ggAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGGAATCAGTCt                                      >  1:406177/1‑93 (MQ=255)
                                                   ggAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGAGACTGTAAAACGCTGCAGCCAAAGTCTGCCTTTGGCTGCAATGAAATCAGTCt                                      >  2:408993/1‑93 (MQ=255)
                                                    gAAAACAGAAGGCGGTTTTTCAAAGGAAAAGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGTCTGCAATGAAATCAGTCtt                                     <  2:1090135/93‑1 (MQ=255)
                                                      aaaCAGAAGGCGGTTTAGCAGAGGAAAGTCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTa                                   <  1:177687/93‑1 (MQ=255)
                                                      aaaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTa                                   <  2:107733/93‑1 (MQ=255)
                                                      aaaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTa                                   <  2:2166989/93‑1 (MQ=255)
                                                       aaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTGAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTAt                                  >  2:4659861/1‑93 (MQ=255)
                                                       aaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTAt                                  >  2:5134636/1‑93 (MQ=255)
                                                        aCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATg                                 >  1:5057211/1‑93 (MQ=255)
                                                          aGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGcc                               <  1:1001668/93‑1 (MQ=255)
                                                             aGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGt                            >  1:1619140/1‑93 (MQ=255)
                                                             aGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGt                            >  2:1204894/1‑93 (MQ=255)
                                                              ggCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGtt                           <  2:2802357/93‑1 (MQ=255)
                                                                cGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTAt                         <  2:2697574/93‑1 (MQ=255)
                                                                cGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTAt                         <  2:5397348/93‑1 (MQ=255)
                                                                 ggTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATc                        >  2:5596086/1‑93 (MQ=255)
                                                                 ggTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATc                        >  2:1490322/1‑93 (MQ=255)
                                                                  gTTTACCAGAGTAAACGCTGATTGTGACTGTAAAACGATGCAGCCAAAGTGTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATca                       <  1:5517555/93‑1 (MQ=255)
                                                                  gTTTAACAGAGGAAAGGCTTATGGTGAATGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATca                       >  2:641085/1‑93 (MQ=255)
                                                                   tttAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcag                      <  2:5710167/93‑1 (MQ=255)
                                                                     tAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggc                    >  1:3882038/1‑93 (MQ=255)
                                                                      aaCAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggca                   <  2:4243381/93‑1 (MQ=255)
                                                                      aaCAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggca                   >  2:1040538/1‑93 (MQ=255)
                                                                       aCAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggcag                  >  2:3850430/1‑93 (MQ=255)
                                                                       aCAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggcag                  <  1:5649583/93‑1 (MQ=255)
                                                                         agagGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGa                >  2:5479082/1‑93 (MQ=255)
                                                                         agagGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGa                <  2:3326350/93‑1 (MQ=255)
                                                                          gagGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGaa               >  1:2432215/1‑93 (MQ=255)
                                                                          gagGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGCTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGaa               >  1:5421070/1‑93 (MQ=255)
                                                                            ggAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAata             <  2:3775063/93‑1 (MQ=255)
                                                                            ggAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAata             >  1:2020488/1‑93 (MQ=255)
                                                                            ggAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAata             >  2:2517633/1‑93 (MQ=255)
                                                                             gAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAatat            <  1:2349859/93‑1 (MQ=255)
                                                                             gAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAatat            <  1:2578881/93‑1 (MQ=255)
                                                                              aaaGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATg           <  1:5752980/93‑1 (MQ=255)
                                                                              aaaGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATg           <  1:5886899/93‑1 (MQ=255)
                                                                              aaaGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATg           <  1:653177/93‑1 (MQ=255)
                                                                               aaGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCGTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATgg          >  2:3422233/1‑93 (MQ=255)
                                                                                    ttATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAAc     >  2:3372023/1‑93 (MQ=255)
                                                                                     tATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGTCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACa    <  2:2253796/93‑1 (MQ=255)
                                                                                     tATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACa    <  1:4404400/93‑1 (MQ=255)
                                                                                     tATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACa    >  1:5566998/1‑93 (MQ=255)
                                                                                      aTGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACac   >  2:3764776/1‑92 (MQ=255)
                                                                                       tGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATc  <  2:195904/93‑1 (MQ=255)
                                                                                       tGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATc  >  2:5187733/1‑93 (MQ=255)
                                                                                       tGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATc  >  1:1277468/1‑93 (MQ=255)
                                                                                          |                                                                                         
ATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATC  >  NC_000913/502563‑502742

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