Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,127,778 C→T G239G (GGC→GGT mviN → putative virulence factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,127,7780CT98.7% 792.0 / ‑3.5 227G239G (GGC→GGTmviNputative virulence factor
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (118/106);  total (120/107)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-01

Only 100 of 235 total aligned reads displayed.

GAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGCCGA  >  NC_000913/1127688‑1127867
                                                                                          |                                                                                         
   gaTCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGtctc                                                                                      >  2:9816559/1‑93 (MQ=255)
     tCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCtg                                                                                    <  2:2655523/93‑1 (MQ=255)
      cGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCtgt                                                                                   <  2:4707363/93‑1 (MQ=255)
       ggCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCtgtg                                                                                  <  1:6827533/93‑1 (MQ=255)
          aTGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGc                                                                               <  2:3090142/93‑1 (MQ=255)
           tGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGcc                                                                              >  1:1038166/1‑93 (MQ=255)
           tGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGcc                                                                              >  1:2173298/1‑93 (MQ=255)
           tGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGcc                                                                              >  1:3889148/1‑93 (MQ=255)
           tGCTGGTCCGGCCGCGCATTAACTTCTACGATGCCGTAGGGCTTCGCGTGGTGTAACAGATGGGACCTTCGATCCTTGGTGGCTCTGTGAgtc                                                                              >  2:6027840/1‑91 (MQ=255)
              tGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGa                                                                           >  2:5631065/1‑93 (MQ=255)
              tGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGa                                                                           <  2:645333/93‑1 (MQ=255)
               ggTCCTGCTGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGAt                                                                          <  1:4593948/93‑1 (MQ=255)
               ggTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGAt                                                                          <  2:7605440/93‑1 (MQ=255)
               ggTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGAt                                                                          <  2:7842193/93‑1 (MQ=255)
                gTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGAtc                                                                         >  1:4308403/1‑93 (MQ=255)
                gTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGAtc                                                                         <  1:6347508/93‑1 (MQ=255)
                   cTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTcc                                                                      <  1:2148993/93‑1 (MQ=255)
                    tGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCt                                                                     >  1:7209796/1‑93 (MQ=255)
                    tGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCt                                                                     <  2:7041647/93‑1 (MQ=255)
                     gCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCtt                                                                    >  1:7522935/1‑93 (MQ=255)
                     gCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCtt                                                                    <  2:9419662/93‑1 (MQ=255)
                     gCCGCGCATTAACTTCCACGATGCCGGAGCAATGCACGTGGTGAAACGGATGGGACCGGCGATCCTTGGTGTCCCTGTGAGCCAGATCTCCtt                                                                    <  2:3436670/93‑1 (MQ=255)
                      ccGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTa                                                                   >  1:5866375/1‑93 (MQ=255)
                        gcgcATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAat                                                                 >  2:1806013/1‑93 (MQ=255)
                        gcgcATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAat                                                                 >  1:8719392/1‑93 (MQ=255)
                           cATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAatcat                                                              >  2:6356019/1‑93 (MQ=255)
                           cATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAatcat                                                              <  1:8981974/93‑1 (MQ=255)
                             ttAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAAtcatca                                                            >  1:3479375/1‑93 (MQ=255)
                             ttAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAAtcatca                                                            >  1:9601869/1‑93 (MQ=255)
                             ttAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAAtcatca                                                            >  1:3182441/1‑93 (MQ=255)
                              tAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCaa                                                           >  2:4441446/1‑93 (MQ=255)
                               aaCTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAac                                                          <  1:9811142/93‑1 (MQ=255)
                                aCTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAaca                                                         >  2:2440123/1‑93 (MQ=255)
                                aCTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAaca                                                         <  1:8290056/93‑1 (MQ=255)
                                aCTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAaca                                                         <  2:5933328/93‑1 (MQ=255)
                                aCTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAaca                                                         <  1:7165489/93‑1 (MQ=255)
                                  ttCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACAcc                                                       <  2:9457011/93‑1 (MQ=255)
                                   tCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCa                                                      >  1:9173679/1‑93 (MQ=255)
                                   tCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCa                                                      >  1:5692128/1‑93 (MQ=255)
                                     cACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAtt                                                    >  2:9747046/1‑93 (MQ=255)
                                     cACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAtt                                                    >  2:5438782/1‑93 (MQ=255)
                                      aCGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAttt                                                   >  1:3334987/1‑93 (MQ=255)
                                      aCGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAttt                                                   <  1:9409092/93‑1 (MQ=255)
                                      aCGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAttt                                                   >  1:2389024/1‑93 (MQ=255)
                                      aCGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAttt                                                   <  2:8029917/93‑1 (MQ=255)
                                        gATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAttttt                                                 <  2:3550908/93‑1 (MQ=255)
                                        gATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGCGAGCCAGATCTCCTTAATCATCAACACCAttttt                                                 >  1:6215150/1‑93 (MQ=255)
                                          tGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGc                                               >  2:6855478/1‑93 (MQ=255)
                                           gCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGcc                                              >  2:8840866/1‑93 (MQ=255)
                                           gCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGcc                                              <  1:6299060/93‑1 (MQ=255)
                                           gCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGcc                                              <  1:8633754/93‑1 (MQ=255)
                                            ccGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCt                                             <  1:4740048/93‑1 (MQ=255)
                                               gAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGt                                          <  1:8954906/93‑1 (MQ=255)
                                               gAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGt                                          >  2:7171680/1‑93 (MQ=255)
                                                 gCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGttt                                        <  2:3700177/93‑1 (MQ=255)
                                                  cAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTc                                       <  1:1962766/93‑1 (MQ=255)
                                                  cAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTc                                       <  1:6423293/93‑1 (MQ=255)
                                                   aaTGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCt                                      >  2:9466277/1‑93 (MQ=255)
                                                   aaTGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCt                                      >  2:5476577/1‑93 (MQ=255)
                                                    aTGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCtt                                     >  2:3280278/1‑93 (MQ=255)
                                                    aTGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCtt                                     >  2:6529943/1‑93 (MQ=255)
                                                    aTGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCtt                                     >  1:9336681/1‑93 (MQ=255)
                                                     tGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTg                                    >  2:9556669/1‑93 (MQ=255)
                                                      gcgcGTGGTGGAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACTCCATTTTTGCCTCGTTTCTTGc                                   <  2:4862317/93‑1 (MQ=255)
                                                      gcgcGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGc                                   <  1:2739105/93‑1 (MQ=255)
                                                        gcgTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCtt                                 >  1:9800583/1‑93 (MQ=255)
                                                         cgTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTc                                >  1:6018015/1‑93 (MQ=255)
                                                         cgTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTc                                <  1:5024781/93‑1 (MQ=255)
                                                         cgTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTc                                <  2:7595252/93‑1 (MQ=255)
                                                           tggtgAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCg                              <  2:5145863/93‑1 (MQ=255)
                                                           tggtgAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCg                              >  1:4879716/1‑93 (MQ=255)
                                                             gtgAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGt                            >  2:4609761/1‑93 (MQ=255)
                                                               gAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTc                          <  1:2943569/93‑1 (MQ=255)
                                                               gAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTc                          >  2:7412649/1‑93 (MQ=255)
                                                               gAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTc                          >  1:1604760/1‑93 (MQ=255)
                                                               gAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTc                          <  2:198627/93‑1 (MQ=255)
                                                                aaaCAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCg                         >  1:6965909/1‑93 (MQ=255)
                                                                  aCAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGgt                       <  2:4906960/93‑1 (MQ=255)
                                                                  aCAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGgt                       >  2:2251972/1‑93 (MQ=255)
                                                                   cAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGgtg                      >  2:5145459/1‑93 (MQ=255)
                                                                    aGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGgtgt                     <  2:31587/93‑1 (MQ=255)
                                                                       tGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCtt                  <  2:5765157/93‑1 (MQ=255)
                                                                        gggACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTg                 >  1:4210308/1‑93 (MQ=255)
                                                                          gACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGa               >  2:3712569/1‑93 (MQ=255)
                                                                          gACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCACCTTAATCATCAACCCCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGa               <  2:1688554/93‑1 (MQ=255)
                                                                           aCCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGAt              >  2:2927209/1‑93 (MQ=255)
                                                                            ccGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATg             >  2:9298632/1‑93 (MQ=255)
                                                                              ggCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTa           >  1:6776986/1‑93 (MQ=255)
                                                                              ggCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGCTCGGTGTCTTGGCTGTa           >  1:3986711/1‑93 (MQ=255)
                                                                               gCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTAt          <  2:5765645/93‑1 (MQ=255)
                                                                               gCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTAt          <  1:6455877/93‑1 (MQ=255)
                                                                                 gATCCTTTGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTa        <  2:934989/93‑1 (MQ=255)
                                                                                 gATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTa        <  2:8627095/93‑1 (MQ=255)
                                                                                  aTCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTAc       <  1:6167114/93‑1 (MQ=255)
                                                                                  aTCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTAc       >  2:5010402/1‑93 (MQ=255)
                                                                                   tCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACg      >  1:5425567/1‑93 (MQ=255)
                                                                                    ccTTGGTGTCTCTCTCAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTGGCTTCCGGTTCGGTGTCTTGGATGTATTACGc     <  2:3982771/93‑1 (MQ=255)
                                                                                      ttGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGCCg   >  2:8415278/1‑93 (MQ=255)
                                                                                      ttGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGCCg   >  2:8687802/1‑93 (MQ=255)
                                                                                      ttGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACAATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGCCg   >  2:7291033/1‑93 (MQ=255)
                                                                                          |                                                                                         
GAAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGCCGA  >  NC_000913/1127688‑1127867

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