Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,313,777 G→A 5.5% A143T (GCA→ACA)  rpiB → ribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,313,7770GA 5.5% 209.0 / 2.2 72A143T (GCA→ACA) rpiBribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (40/28);  new base A (2/2);  total (42/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

GCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCG  >  NC_000913/4313629‑4313925
                                                                                                                                                    |                                                                                                                                                    
gcAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACgg                                                                                                                                                      <  1:234831/149‑1 (MQ=255)
  aaCTTTGGCGGCAGCACAACGACACCAATGTGGTGCCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCa                                                                                                                                                    <  2:558173/149‑1 (MQ=255)
    cTTTCGCAGCAGCATAACGACACAAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCACCAACGCGTGGAGGCGATTACGGCAAt                                                                                                                                                  >  2:513706/1‑149 (MQ=255)
         gcgGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGAAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAAACGagc                                                                                                                                             >  2:70055/1‑149 (MQ=255)
           ggCAGCATAACGACACCAACGAGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTCGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGCGGCGATCACGGCAATAGcgcag                                                                                                                                           >  2:428457/1‑149 (MQ=255)
               gcaTAACGAAACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGga                                                                                                                                       <  2:715913/149‑1 (MQ=255)
                 aTAACTACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGgaga                                                                                                                                     >  2:466908/1‑149 (MQ=255)
                 aTAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGgaga                                                                                                                                     >  2:509972/1‑149 (MQ=255)
                 aTAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGGCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGgaga                                                                                                                                     >  2:641269/1‑149 (MQ=255)
                    aCGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCATGGCTGGGCGCCCCGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGGGAAAt                                                                                                                                  >  1:595249/1‑149 (MQ=255)
                        cacCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGCTTGTGGTTTTGTGGTTGGTCTCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTgag                                                                                                                              <  2:269220/149‑1 (MQ=255)
                               gTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATc                                                                                                                       <  1:17257/149‑1 (MQ=255)
                                     gCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCactac                                                                                                                 >  2:268725/1‑149 (MQ=255)
                                     gCTTTTGGTTCACGAGGGGTTGGCCTCGAACTGGCAAAAATGGTTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGATCAGCGGAGAAATTGGGATTCATCCgctac                                                                                                                 >  1:4327/1‑149 (MQ=255)
                                       ttttGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGGTGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGGAATTGAGATTGATCCACTACtt                                                                                                               >  2:458365/1‑149 (MQ=255)
                                       ttttGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACtt                                                                                                               >  2:654621/1‑149 (MQ=255)
                                       ttttGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTg                                                                                                              >  1:419062/1‑149 (MQ=255)
                                          ggggtCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCa                                                                                                            <  2:309194/145‑1 (MQ=255)
                                           ggTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCAt                                                                                                           <  1:297653/149‑1 (MQ=255)
                                             ttCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGAATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATgg                                                                                                         >  2:725180/1‑149 (MQ=255)
                                               cACGGGTGGTTGGCCTCGAACTGGCAACCATGATTGTGGATGCGTGGTTGGGCGCATCGTACGAAGGCGGTCGTCATCAATAACGCGTGTAAGCTATGTCCGCAATAGAGCGGCGGAGAGATTGAGATTCATCCACTACTTGCATGGAt                                                                                                       <  2:178924/149‑1 (MQ=255)
                                               cACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGACGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGAt                                                                                                       <  1:242298/149‑1 (MQ=255)
                                                 cGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGGAGTTGAGATTCATACACTACTTGCATGGATGa                                                                                                     >  2:384181/1‑149 (MQ=255)
                                                    gTGGTTGGCCTCGAACTGGCAAACATGATCGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAAAAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTa                                                                                                  <  2:473015/149‑1 (MQ=255)
                                                    gTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCTGTCGTCATCAACAACGCGTGGAGGTGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTa                                                                                                  <  2:321556/149‑1 (MQ=255)
                                                     tGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTaa                                                                                                 >  1:346312/1‑149 (MQ=255)
                                                          ggCCTCGAACTGGCAAAAAGGATTGTGGATGCGTGGCTGGGCGCACGGTACGAAGGCGGTCGTCATCAACCACGCGTGGAGGCGATTACGGCAATAGAACATCGGAGAAAGAGAGATTCATCCACTACTTGCATGCATGAGTACTGAtc                                                                                            >  2:112497/1‑148 (MQ=255)
                                                            ccTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGCAGGAGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGGAAGTGAGATTCATCCACTACGTGCATGGCTGGGAGAGGATTaa                                                                                          >  1:31413/1‑149 (MQ=255)
                                                            ccTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTaa                                                                                          >  2:23184/1‑149 (MQ=255)
                                                             cTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCGGCGGAGAAATTGAGATTCATCCACTACTTGCATGGAAGaggaaggacgaat                                                                                         >  2:421012/1‑144 (MQ=255)
                                                              tCGACCTGGCAAAAATGATTGTGGATGCGTGCGTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAAtg                                                                                        <  1:505009/149‑1 (MQ=255)
                                                              tCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATTCACTACTTGCATGGATGAGTAATGCTTAAtg                                                                                        >  2:44160/1‑149 (MQ=255)
                                                              tCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAGtg                                                                                        >  2:597646/1‑149 (MQ=255)
                                                              tCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAAtg                                                                                        >  2:85568/1‑149 (MQ=255)
                                                              tCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAAtg                                                                                        >  1:672315/1‑149 (MQ=255)
                                                               cGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAAtgt                                                                                       <  2:673737/149‑1 (MQ=255)
                                                                    tGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCGCAGTACGAAGGCGGTAGTAATCAACACCGCGTGGAGGGGATTGCGGCACAAGAGCAGAGGAGAAATTGAGTTTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGAtg                                                                                  >  2:192278/1‑148 (MQ=255)
                                                                    tGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTGCTGCACTACTTGCATTGATGATTAATGATTAATGTGGATa                                                                                  <  2:690603/149‑1 (MQ=255)
                                                                      gCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTGATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATGGAGATTTATCCACTACTTGCAAGGGTGAGTAATGATTAATGTGGATaga                                                                                >  2:440415/1‑149 (MQ=255)
                                                                      gCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGACAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATaga                                                                                >  1:426966/1‑149 (MQ=255)
                                                                      gCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGACAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATaga                                                                                >  1:427202/1‑149 (MQ=255)
                                                                             tGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCCACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGGAATGATTAATGTGGCTAGAGTTTCtt                                                                         >  2:233649/1‑149 (MQ=255)
                                                                                ttGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCttttt                                                                      <  1:654621/149‑1 (MQ=255)
                                                                                  gtgGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGTATAGAGTTTCTTTTTGa                                                                    <  1:78418/149‑1 (MQ=255)
                                                                                   tgGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAg                                                                   <  1:641339/149‑1 (MQ=255)
                                                                                     gATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGt                                                                 >  2:605398/1‑149 (MQ=255)
                                                                                          gTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGAAATAGCGAAGCGGCGAAATTGAGATTCATCCACTACTTGGCTGCATGTGTGATGGAACATGTGGAGAAAGTTACTTTTTGAggctggct                                                            >  2:248922/1‑149 (MQ=255)
                                                                                               tGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTCAtaa                                                       >  2:175322/1‑147 (MQ=255)
                                                                                                  gcgcACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGGTTCTTTTTGAGGTTGGCTAATGAAAc                                                    >  1:243524/1‑149 (MQ=255)
                                                                                                     cacaGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCt                                                 <  1:429860/149‑1 (MQ=255)
                                                                                                          tACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGGTGCTTGTTGAGGTCGGCTAATGAAACggcttcat                                            >  2:371995/1‑142 (MQ=255)
                                                                                                          tACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGACAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCg                                            <  2:426966/149‑1 (MQ=255)
                                                                                                          tACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGACAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCg                                            <  2:427202/149‑1 (MQ=255)
                                                                                                            cGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGGAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCt                                          >  1:159920/1‑149 (MQ=255)
                                                                                                             gAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTa                                         >  1:602763/1‑149 (MQ=255)
                                                                                                                 gCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTc                                     <  2:318917/149‑1 (MQ=255)
                                                                                                                    gtcgtcATCAACAACGCGTGGTGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAGACGCTTTCCGCTATTTCttt                                  >  1:533047/1‑149 (MQ=255)
                                                                                                                    gtcgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCttt                                  <  1:562150/149‑1 (MQ=255)
                                                                                                                     tcgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGGTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCtttt                                 >  1:575022/1‑149 (MQ=255)
                                                                                                                     tcgtcAGCAACAACCCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCtttt                                 <  2:694004/149‑1 (MQ=255)
                                                                                                                       gtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTAt                               <  1:314515/149‑1 (MQ=255)
                                                                                                                          atcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGAGAGAGGTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTGTTTCTTTTATTCa                            >  1:78955/1‑149 (MQ=255)
                                                                                                                          accaacaaCGCGTGGAGGGGATTACGGCAATTGAGCAGCGGAGAAATTGAGATTCATCCCCTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTCGGCTGATGAAACGCTTTCCGCTATTTCTTTTATTCa                            <  2:583407/147‑1 (MQ=255)
                                                                                                                           tcaACAACGCGGGGAGGCGATTACGGCAATATGGCAGCGGAGTAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTACTTCTTTTATTCAc                           <  2:458159/149‑1 (MQ=255)
                                                                                                                             aacaacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCAccc                         >  1:112529/1‑149 (MQ=255)
                                                                                                                              acaacGCGTGGAGGCGATTACGGCAATAGAGCAGACGCGTAAATGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAAGGCTTTCCTCTATTTCTTTTATTCACCCt                        >  1:180699/1‑149 (MQ=255)
                                                                                                                                 gcgcgTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCt                     <  1:156356/148‑1 (MQ=255)
                                                                                                                                 gcgcgTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCt                     <  2:729156/148‑1 (MQ=255)
                                                                                                                                 acGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCt                     >  1:32292/1‑149 (MQ=255)
                                                                                                                                          ggCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTc            <  1:557659/149‑1 (MQ=255)
                                                                                                                                           gCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTcc           <  1:725180/149‑1 (MQ=255)
                                                                                                                                            cGATTACGGCAATCGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCa          <  2:224095/149‑1 (MQ=255)
                                                                                                                                             gATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTcctg         >  2:448366/1‑147 (MQ=255)
                                                                                                                                                   ggCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTcc   >  1:189358/1‑149 (MQ=255)
                                                                                                                                                    gCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTCCCg  >  2:364491/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCAACTTTCGCGGCAGCATAACGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCG  >  NC_000913/4313629‑4313925

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