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.0% 511.2 / ‑1.9 151G239G (GGC→GGTmviNputative virulence factor
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base T (66/82);  total (69/82)
Fisher's exact test for biased strand distribution p-value = 9.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 157 total aligned reads displayed.

AAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGCCG  >  NC_000913/1127689‑1127866
                                                                                         |                                                                                        
 agaTCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGtct                                                                                      >  2:806911/1‑93 (MQ=255)
 agaTCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGtct                                                                                      >  2:3394197/1‑93 (MQ=255)
  gaTCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGGAATGCGCGGGGGGTAACAGATGGGACCGGCGATCCTTGGTGtctc                                                                                     >  1:1016990/1‑93 (MQ=255)
  gaTCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGtctc                                                                                     >  1:1027595/1‑93 (MQ=255)
  gaTCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGtctc                                                                                     >  1:4165163/1‑93 (MQ=255)
     cGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCtgt                                                                                  <  1:236737/93‑1 (MQ=255)
      ggCATGCTGGTCGTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCtgtg                                                                                 <  1:1466389/93‑1 (MQ=255)
      ggCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCtgtg                                                                                 <  2:1524143/93‑1 (MQ=255)
       gCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGa                                                                                <  1:1245169/93‑1 (MQ=255)
       gCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGa                                                                                <  1:2597325/93‑1 (MQ=255)
       gCAGGCTGGTCCTGCCGCGGATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGa                                                                                <  1:3567520/93‑1 (MQ=255)
         aTGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGc                                                                              >  2:881422/1‑93 (MQ=255)
         aTGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGc                                                                              >  1:4099606/1‑93 (MQ=255)
          tGCTGGTCCTGCCGCGCATTAACTTCGACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGcc                                                                             <  2:1307408/93‑1 (MQ=255)
          tGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGcc                                                                             <  2:2012541/93‑1 (MQ=255)
             tGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATTGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGa                                                                          <  1:3492412/93‑1 (MQ=255)
              ggTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGAt                                                                         <  2:1807211/93‑1 (MQ=255)
                  cTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTcc                                                                     <  1:3278166/93‑1 (MQ=255)
                   tGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCt                                                                    >  1:2288864/1‑93 (MQ=255)
                   tGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCt                                                                    <  1:1801234/93‑1 (MQ=255)
                    gCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCtt                                                                   >  2:4480199/1‑93 (MQ=255)
                    gCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCtt                                                                   <  2:855534/93‑1 (MQ=255)
                    gCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCtt                                                                   <  1:3398426/93‑1 (MQ=255)
                       gcgcATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATTCTTGGTGTCTCTGGGAGTCAGATCTCCTTAat                                                                >  2:1639809/1‑93 (MQ=255)
                        cgcATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAatc                                                               >  1:4669234/1‑93 (MQ=255)
                           aTTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAatcatc                                                            >  1:3316450/1‑93 (MQ=255)
                             tAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCaa                                                          >  2:3372110/1‑93 (MQ=255)
                             tAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCaa                                                          <  1:1510489/93‑1 (MQ=255)
                             tAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCaa                                                          >  1:1588035/1‑93 (MQ=255)
                              aaCTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAac                                                         >  1:1953905/1‑93 (MQ=255)
                              aaCTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAac                                                         <  1:1338505/93‑1 (MQ=255)
                                cTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAacac                                                       <  1:1389447/93‑1 (MQ=255)
                                 ttCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACAcc                                                      <  2:2908438/93‑1 (MQ=255)
                                 ttCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACAcc                                                      >  1:4565660/1‑93 (MQ=255)
                                 ttCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACAcc                                                      >  1:358497/1‑93 (MQ=255)
                                  tCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCa                                                     >  1:987158/1‑93 (MQ=255)
                                   ccACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAt                                                    <  1:929498/93‑1 (MQ=255)
                                   ccACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAt                                                    >  2:472771/1‑93 (MQ=255)
                                     aCGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAttt                                                  >  1:2037933/1‑93 (MQ=255)
                                     aCGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGAACGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAttt                                                  >  1:2656215/1‑93 (MQ=255)
                                      cGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAtttt                                                 <  1:2895869/93‑1 (MQ=255)
                                       gATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAttttt                                                >  1:34265/1‑93 (MQ=255)
                                       gATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCAttttt                                                <  1:381606/93‑1 (MQ=255)
                                        aTGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTg                                               >  1:1955787/1‑93 (MQ=255)
                                           ccGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCt                                            >  2:841715/1‑93 (MQ=255)
                                            cGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTc                                           >  2:4354926/1‑93 (MQ=255)
                                            cGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTCTGCCTc                                           >  1:866938/1‑93 (MQ=255)
                                            cGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTc                                           >  1:4012857/1‑93 (MQ=255)
                                              gAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGt                                         <  1:17361/93‑1 (MQ=255)
                                              gAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATACTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGt                                         >  2:736033/1‑93 (MQ=255)
                                               aGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGtt                                        <  1:2758086/93‑1 (MQ=255)
                                               aGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGtt                                        <  1:2902171/93‑1 (MQ=255)
                                               aGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGtt                                        >  2:3957960/1‑93 (MQ=255)
                                               aGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGtt                                        >  1:3442725/1‑93 (MQ=255)
                                                gCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGttt                                       >  1:2690328/1‑93 (MQ=255)
                                                gCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGttt                                       >  1:3198418/1‑93 (MQ=255)
                                                gCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGttt                                       <  1:2719796/93‑1 (MQ=255)
                                                 cAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTc                                      <  2:1372023/93‑1 (MQ=255)
                                                 cAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTc                                      <  2:3506454/93‑1 (MQ=255)
                                                  aaTGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCt                                     >  1:2924088/1‑93 (MQ=255)
                                                   aTGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCtt                                    >  1:4720431/1‑93 (MQ=255)
                                                   aTGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCtt                                    <  1:2221338/93‑1 (MQ=255)
                                                     gcgcGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGc                                  <  2:3673186/93‑1 (MQ=255)
                                                     gcgcGTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGc                                  <  1:1748946/93‑1 (MQ=255)
                                                      cgcgTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCt                                 <  2:2984518/93‑1 (MQ=255)
                                                       gggtggtgAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGAGAGCCAGATCTCCTTAATCATCAACACCATGTTTGCCTCGTTTCTTGCtg                                >  2:2339284/3‑92 (MQ=255)
                                                       gcgTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCtt                                >  1:2675770/1‑93 (MQ=255)
                                                        cgTGGTGAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTc                               <  1:3560127/93‑1 (MQ=255)
                                                          tggtgAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATAATCAACACCATTTTTGCCTCGTTTCTTGCTTCCg                             <  1:1805127/93‑1 (MQ=255)
                                                           ggtgAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCgg                            >  2:4157618/1‑93 (MQ=255)
                                                            gtgAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGt                           <  1:3624669/93‑1 (MQ=255)
                                                             tgAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGtt                          <  2:573391/93‑1 (MQ=255)
                                                             ggAAACAGATGGAACCGGCGATACTTGGTGTCCCTGTGAGCCAGATCTCCTTACTCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGtt                          <  2:3579531/92‑1 (MQ=255)
                                                              gAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAATACCATTTTTGCCTCGTTTCTTGCTTCCGGTTc                         <  1:4292625/93‑1 (MQ=255)
                                                              gAAACAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTc                         >  2:3028851/1‑93 (MQ=255)
                                                               aaaCAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCg                        <  1:2844207/93‑1 (MQ=255)
                                                               aaaCAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCg                        <  1:2111818/93‑1 (MQ=255)
                                                               aaaCAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCg                        <  2:1084364/93‑1 (MQ=255)
                                                                aaCAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCgg                       >  2:3750426/1‑93 (MQ=255)
                                                                aaCAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCgg                       <  2:3276651/93‑1 (MQ=255)
                                                                 aCAGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGgt                      >  1:225307/1‑93 (MQ=255)
                                                                   aGATGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGgtgt                    <  1:2109403/93‑1 (MQ=255)
                                                                      tGGGACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCtt                 >  1:4301814/1‑93 (MQ=255)
                                                                       gggACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTg                >  1:3814761/1‑93 (MQ=255)
                                                                         gACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGa              <  1:1413428/93‑1 (MQ=255)
                                                                         gACCGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGa              <  1:3114974/93‑1 (MQ=255)
                                                                            cGGCGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGt           >  1:826985/1‑93 (MQ=255)
                                                                               cGATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTAtt        <  1:4372518/93‑1 (MQ=255)
                                                                               cGATCCTTGGTGGCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTAtt        <  1:3405455/93‑1 (MQ=255)
                                                                                gATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACGCCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTa       <  1:589757/93‑1 (MQ=255)
                                                                                gATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTa       <  1:589756/93‑1 (MQ=255)
                                                                                gATCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTa       <  2:2752479/93‑1 (MQ=255)
                                                                                 aTCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTAc      <  1:4640781/93‑1 (MQ=255)
                                                                                 aTCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTAc      >  2:1640755/1‑93 (MQ=255)
                                                                                 aTCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTAc      >  2:2733238/1‑93 (MQ=255)
                                                                                 aTCCTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTAc      >  2:2057894/1‑93 (MQ=255)
                                                                                  tCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTTCTTACAGTTCGGTGTCTTGGATGTATTACg     >  2:3115959/1‑93 (MQ=255)
                                                                                   ccTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGc    <  1:2762434/93‑1 (MQ=255)
                                                                                    cTTGGTGTGTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGcc   <  2:2394535/93‑1 (MQ=255)
                                                                                    cTTGGTGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGcc   >  2:133141/1‑93 (MQ=255)
                                                                                         |                                                                                        
AAGATCGGCATGCTGGTCCTGCCGCGCATTAACTTCCACGATGCCGGAGCAATGCGCGTGGTGAAACAGATGGGACCGGCGATCCTTGGCGTCTCTGTGAGCCAGATCTCCTTAATCATCAACACCATTTTTGCCTCGTTTCTTGCTTCCGGTTCGGTGTCTTGGATGTATTACGCCG  >  NC_000913/1127689‑1127866

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