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,6530CT98.9% 638.3 / ‑5.2 181intergenic (‑191/+47)ybaL/fsrputative transport protein/fosmidomycin resistance protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (101/78);  minor base G (2/0);  total (103/78)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

ATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATC  >  NC_000913/502563‑502742
                                                                                          |                                                                                         
 ttattGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGAc                                                                                        >  1:3673243/1‑93 (MQ=255)
 ttattGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGAc                                                                                        >  1:967549/1‑93 (MQ=255)
  tattGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACt                                                                                       >  2:1402170/1‑93 (MQ=255)
   attGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTg                                                                                      <  2:578505/93‑1 (MQ=255)
    ttGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGt                                                                                     >  2:1080728/1‑93 (MQ=255)
    ttGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGt                                                                                     <  1:2177685/93‑1 (MQ=255)
     tGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTa                                                                                    >  1:2622140/1‑93 (MQ=255)
      gAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaa                                                                                   <  2:1148352/93‑1 (MQ=255)
        aTTTTACAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaaaa                                                                                 <  1:421322/93‑1 (MQ=255)
        aTTTGGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTaaaa                                                                                 <  1:4194260/93‑1 (MQ=255)
         ttttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAAc                                                                                >  1:629489/1‑93 (MQ=255)
         ttttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAAc                                                                                >  2:2972560/1‑93 (MQ=255)
         ttttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAAc                                                                                <  1:4636362/93‑1 (MQ=255)
           ttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGa                                                                              >  1:3389557/1‑93 (MQ=255)
           ttGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGa                                                                              <  2:3558206/93‑1 (MQ=255)
              cAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGc                                                                           <  1:1909225/93‑1 (MQ=255)
              cAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGc                                                                           <  1:2898973/93‑1 (MQ=255)
               aGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAATACGATGCa                                                                          <  2:4575621/93‑1 (MQ=255)
               aGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCa                                                                          <  1:3275637/93‑1 (MQ=255)
               aGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCa                                                                          >  1:3841292/1‑93 (MQ=255)
                gAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAg                                                                         <  1:452448/93‑1 (MQ=255)
                gAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAg                                                                         >  2:1938020/1‑93 (MQ=255)
                 aaaaaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGc                                                                        <  1:3923552/93‑1 (MQ=255)
                    aaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCaa                                                                     >  1:1350571/1‑93 (MQ=255)
                    aaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCaa                                                                     >  2:1668111/1‑93 (MQ=255)
                    aaTGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCaa                                                                     >  2:3096385/1‑93 (MQ=255)
                       ggAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGt                                                                  >  2:4071456/1‑93 (MQ=255)
                        gAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGtt                                                                 <  1:1688899/93‑1 (MQ=255)
                         aaTTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGttt                                                                <  1:3781950/93‑1 (MQ=255)
                         aaTTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGttt                                                                >  2:1236845/1‑93 (MQ=255)
                           ttAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGc                                                              <  2:1621731/93‑1 (MQ=255)
                           ttAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGc                                                              >  1:2404244/1‑93 (MQ=255)
                            tAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGcc                                                             <  2:3187637/93‑1 (MQ=255)
                              gCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCtt                                                           >  2:4356116/1‑93 (MQ=255)
                              gCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGATGACTGTAAAACGATGCAGCCAAAGTTTGCCtt                                                           >  2:3630603/1‑93 (MQ=255)
                               cTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCttt                                                          <  1:2001883/93‑1 (MQ=255)
                                tGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTg                                                         >  1:1428942/1‑93 (MQ=255)
                                 gCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTgg                                                        >  2:2079067/1‑93 (MQ=255)
                                 gCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTgg                                                        >  2:2471259/1‑93 (MQ=255)
                                 gCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTgg                                                        <  1:1826004/93‑1 (MQ=255)
                                 gCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTgg                                                        <  1:646785/93‑1 (MQ=255)
                                  cAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGc                                                       >  1:4641390/1‑93 (MQ=255)
                                   aaaaaaaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGATTGCCTTTGGCt                                                      >  2:771727/1‑93 (MQ=255)
                                        aaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAAt                                                 <  2:4291231/93‑1 (MQ=255)
                                        aaGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAAt                                                 <  2:500173/93‑1 (MQ=255)
                                          gCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGa                                               <  2:169643/93‑1 (MQ=255)
                                          gCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGa                                               >  2:240139/1‑93 (MQ=255)
                                           cACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGaa                                              >  1:3051592/1‑93 (MQ=255)
                                             cGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAAt                                            <  1:555411/93‑1 (MQ=255)
                                              ggATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATc                                           >  2:2615500/1‑93 (MQ=255)
                                              ggATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATc                                           <  1:1010225/93‑1 (MQ=255)
                                               gATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCa                                          <  2:4639493/93‑1 (MQ=255)
                                                  cGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTc                                       <  1:3221854/93‑1 (MQ=255)
                                                   ggAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCt                                      >  2:2846642/1‑93 (MQ=255)
                                                    gAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCtt                                     >  1:249402/1‑93 (MQ=255)
                                                    gAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCtt                                     <  2:637601/93‑1 (MQ=255)
                                                      aaaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTa                                   <  1:868726/93‑1 (MQ=255)
                                                      aaaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTa                                   <  1:4584299/93‑1 (MQ=255)
                                                       aaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTAt                                  >  1:769450/1‑93 (MQ=255)
                                                       aaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTAt                                  <  1:1570541/93‑1 (MQ=255)
                                                       aaCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTAt                                  >  2:1777265/1‑93 (MQ=255)
                                                        aCAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAACTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATg                                 <  1:3690149/93‑1 (MQ=255)
                                                         cAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGc                                >  2:4533818/1‑93 (MQ=255)
                                                         cAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGc                                >  1:1474752/1‑93 (MQ=255)
                                                         cAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGc                                <  1:165858/93‑1 (MQ=255)
                                                           gAAGGCGGTTTAACAGAGGAAAGGCTTATGTTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGCCTGCAATGAAATCAGTCTTTATGCCg                              <  1:333127/93‑1 (MQ=255)
                                                           gAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCg                              >  1:3213172/1‑93 (MQ=255)
                                                            aaGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCgg                             <  2:3868306/93‑1 (MQ=255)
                                                            aaGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCgg                             >  1:205952/1‑93 (MQ=255)
                                                             aGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGt                            >  1:3551384/1‑93 (MQ=255)
                                                              gggggTTTAAGAGAGAAAGGGCTTATGGGAACTGTAAAACGGTGCCCGAAAAGTTTGCCTTTGGCTGCAAGGAAATCGGTTTTTATCCCGGtt                           <  2:2826577/90‑1 (MQ=255)
                                                              ggCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGtt                           >  2:3365678/1‑93 (MQ=255)
                                                               gCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTa                          <  1:2965444/93‑1 (MQ=255)
                                                               gCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTa                          >  2:1869512/1‑93 (MQ=255)
                                                               gCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTa                          <  1:4167487/93‑1 (MQ=255)
                                                                 ggTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATc                        <  2:2253803/93‑1 (MQ=255)
                                                                  gTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATca                       >  1:2730821/1‑93 (MQ=255)
                                                                  gTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATca                       >  2:2115817/1‑93 (MQ=255)
                                                                  gTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATca                       <  1:56782/93‑1 (MQ=255)
                                                                      aaCAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggca                   <  1:401831/93‑1 (MQ=255)
                                                                        cAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATcaggcagg                 >  2:2948883/1‑93 (MQ=255)
                                                                         agagGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGa                >  1:2646039/1‑93 (MQ=255)
                                                                          gagGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGaa               >  1:2349667/1‑93 (MQ=255)
                                                                           agGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAat              <  1:2234193/93‑1 (MQ=255)
                                                                            ggAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAata             <  1:4140336/93‑1 (MQ=255)
                                                                             gAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAatat            <  2:932227/93‑1 (MQ=255)
                                                                             gAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAatat            <  2:3468774/93‑1 (MQ=255)
                                                                             gAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAatat            >  1:2424857/1‑93 (MQ=255)
                                                                              aaaGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATg           <  1:1546997/93‑1 (MQ=255)
                                                                                aGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGt         >  2:3015875/1‑93 (MQ=255)
                                                                                 ggCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTc        >  1:4307977/1‑93 (MQ=255)
                                                                                    ttATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAAc     >  2:3798789/1‑93 (MQ=255)
                                                                                     tATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACa    <  2:3277336/93‑1 (MQ=255)
                                                                                     tATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACa    <  1:428507/93‑1 (MQ=255)
                                                                                      aTGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACAt   >  1:930751/1‑93 (MQ=255)
                                                                                       tGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATc  >  2:3772255/1‑93 (MQ=255)
                                                                                       tGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATc  >  2:1087672/1‑93 (MQ=255)
                                                                                       tGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATc  <  1:994624/93‑1 (MQ=255)
                                                                                       tGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATc  >  1:1407500/1‑93 (MQ=255)
                                                                                       tGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATc  <  1:1119283/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATC  >  NC_000913/502563‑502742

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