Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 88,571 G→A V182M (GTG→ATG)  cra → DNA‑binding transcriptional dual regulator Cra

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091388,5710GA97.0% 204.8 / ‑5.2 68V182M (GTG→ATG) craDNA‑binding transcriptional dual regulator Cra
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (39/26);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

AGAAGATCAGCCAGACAACGAAATGCGGTGCATTGAGCACCTTTTACAGCGTCAGGTTGATGCCATTATTGTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGGTGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCGCAGGCGCTGTTCACAACGTCGTTTGCGTTGTTGCAAGGAGTGATGGATGTCACGCTGCGTCGCGACGGCA  >  NC_000913/88324‑88814
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
agaagaTCAGCCAGACAACGAAATGCGGTGCATTGAGCACCTTTTACAGCGTCAGGTTGATGCCATTATTGTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGACGGTAGAGTTACGTAAGTTTCCCGCCGTGACGCTGc                                                                                                                                                                                                                                                    >  2:137410/1‑251 (MQ=255)
      tCAGCCAGACAACGAAATGCGGTGCATTGAGCACCTTTTACAGCGTCAGGTTGATGCCATTATTGTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGCAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATc                                                                                                                                                                                                                                              >  1:88578/1‑251 (MQ=255)
         gCCAGGCAACGAAATGCGGTGCGTTGAGCCCTTTTTACAGCGTCGGGGTGATGTCATTATTGTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTAACAACGCTGGGCTAACGACCCGTCCCCGTTTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCCCCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGAAAGTTTCCCGCCGAGACGATGCTTTATCTTg                                                                                                                                                                                                                                           <  1:10318/251‑1 (MQ=255)
          ccAGACAACGAAATGCGGTGCATTGAGCGACTTTTACAGCGTCAGGTTGCTGCCATTATTGTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTgg                                                                                                                                                                                                                                          >  1:196916/1‑251 (MQ=255)
           cAGACAACGAAATGCGGTGCATTGAGCACCTTTTACAGCGTCAGGTTGATGCCATTATTGTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGt                                                                                                                                                                                                                                         >  1:203932/1‑251 (MQ=255)
           aaGACAACGAAAGGGGGTGCATTGAGCACCTTTTACAGCTGCAGTTTTATGCCTTTATTGTGTCGGCGGCGTTGCCTCCTGAGCATCTTTTTTCTCACCGCGGGGCAACCGACCCGTTCCCGATTGTCGCGCGGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGGGGTTGGTGCCGATCAGGATGAGGCCGAAATGCGGGCGGAAGAGTTAGGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGt                                                                                                                                                                                                                                         <  2:199662/250‑1 (MQ=255)
               aaaCGAAATGCGGTGCATTGAGCACCTTTTCCAGCGTCAGGTTGATGCCTTTATTGTTTCGACGTCGTTGCCCCCTGAGCATCCTTTTTACCACCGCGGGGCTAACGACCCGTCCCCGATTGTCGCGCTGTACCGCGCCCTCGATCGTGAACACTTCCCCAGCGGGGTTGGTGCCGATCAGGATGATGCCGAAATGTTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTgcgc                                                                                                                                                                                                                                     <  1:139688/250‑1 (MQ=255)
                     aaggCGGTGCATGGGGCACCTCTCACAGCGCCAGTTTGATGCATTTTTTTTTTCGACGTTGTTGCCCCAGGTGCTGTCTTTTTATCAACGCTGGGTTAACGACCTGGTCCCGATTGTGGCGGGGGACCGACCCCTAGATCGTGAACTCGTCTCTAGCGTGGTTGGGGCCGATCAGGATGATGCCGAAATGCTGGCGAAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCgg                                                                                                                                                                                                                               <  2:491149/248‑1 (MQ=255)
                            tGCATTGAGCACCTTTTACAGCGTCAGGTTGATGCCATTATTGTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTc                                                                                                                                                                                                                        >  1:118883/1‑251 (MQ=255)
                                ttGAGCACCTTTTACAGCGTCAGGTTGATGCCATTATTGTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTc                                                                                                                                                                                                                    >  1:278045/1‑251 (MQ=255)
                                  gAGCACCTTTTACAGCGTCAGGTTGATGCCATTATTGTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTc                                                                                                                                                                                                                    >  1:300928/1‑249 (MQ=255)
                                     ccctttttACAGCGTCAGGTTGATGCCATTATTGTTTCGCGGTCTTTGCCTCCTGAGCACCCTTTTTATCACCGCTGGGCTAACGCCCCGTTCCCGTTTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAATTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCtt                                                                                                                                                                                                               <  1:426041/247‑1 (MQ=255)
                                          tttACAGCGTCAGTTTGATGCCATTTTTGTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTACCGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCCCCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCt                                                                                                                                                                                                            <  1:438080/249‑1 (MQ=255)
                                              cAGCGTCAGGATGATGCCATTATTGTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGCAATGCTGGCGGACGAGTTACGAAAGTTTCCCGCCGAGACTATTCTTTTACTTGGTGCGCTCCCGGAGCTTTCTTTAAGTTTCCTTCGTg                                                                                                                                                                                                       >  2:552074/1‑250 (MQ=255)
                                                  gTCAGGTTGATGCCATTATTGTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGGGACGGTGCTTTATGTTGGGGCGATACCGGAGCTTTCTGTCAGCTTCCTGCgggca                                                                                                                                                                                                     >  2:204222/1‑244 (MQ=255)
                                                       gTTGATGCCATTATTGTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACAGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGttt                                                                                                                                                                                             >  1:142508/1‑251 (MQ=255)
                                                                  tattGTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCt                                                                                                                                                                                  >  1:322066/1‑251 (MQ=255)
                                                                      gTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGa                                                                                                                                                                               >  1:117441/1‑250 (MQ=255)
                                                                        ttCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGaa                                                                                                                                                                              >  1:287210/1‑249 (MQ=255)
                                                                         tCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGATGGAACGCGTCATCGATCGTGAACACTTCACCAGCGTGGTTGATGCCGATCAGGATGATGCCGAAATGCAGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTTCGCTACCGGAGCTTTCTGGCAGCTTCCTGCGTGAACAAGGTTTCCGTACTTCCTGGGAAgc                                                                                                                                                                           >  1:76133/1‑250 (MQ=255)
                                                                           gACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGCGCCGATGCTTTATCTTGGCGCGCTACCAGAGCTTTCTGTCAGCGTCCTGCGTGAACAAGGTTTCCGTACTTCCTGGAAAgctg                                                                                                                                                                         >  2:112112/1‑248 (MQ=255)
                                                                                ctttccctcctGAGCTCCCTTTTTACCACCGCTGGGCTAACTCCCCGTTCCCGTTTGTCCCGCGGCCCCGCGCCCTCGATCGTAACCCCTTCCCCAGCGTGTTTGGTGCCGATCAGGTTGATGCCAAAATCCTGGCGGAAGATTTACTTAAGTTCCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATc                                                                                                                                                                      <  1:253526/244‑1 (MQ=255)
                                                                                cgtTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCcg                                                                                                                                                                    >  2:569686/1‑251 (MQ=255)
                                                                                cgtTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCcg                                                                                                                                                                    >  1:310569/1‑251 (MQ=255)
                                                                                 gtTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCcgc                                                                                                                                                                   >  1:463512/1‑251 (MQ=255)
                                                                                           gcactccatttgtagcAACGCGGGGCGTACGCCCCTTTCCCGATTGTCGCGCTGGCCCGCGCCCTATATCGTGACCATTTCACCAGGGGGGTTGGGGCCGAACAGGAGTATGCCGAATGGCGGGCGGAAGAGCTACGTAATTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTg                                                                                                                                                           <  2:142508/234‑1 (MQ=255)
                                                                                             gcttcttttttATCACCGCGGGGTGAACGCCCCTTCCCCGTTTGGCGCGCTGGACCGCGCCCTCTATCGGGAACATTACACCAGGGTGTTGGGTGCCGATCAGGAGAATGCCGAAAGGCTGGCGGAAGAGTTACGTAATTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCAtt                                                                                                                                                       <  2:310569/245‑1 (MQ=255)
                                                                                             gcaacttttttATCAACGCGGGGCTAACGACCCGTTCCCGATGTTCGCGATGGACCGCGCCCTGGATCGGGAAAATTTCACAAGCGGGGTGGGGGCCGATCAGGAGGATGCCGAAATGTTGGGGGAAGAGTTACGAAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCAtt                                                                                                                                                       <  1:434824/245‑1 (MQ=255)
                                                                                               aTCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGAGGTTGGTGCCGATCAGGATGATGCCGCAATGCTGGCGGAAGAGTTAAGTAAGTTTCCCGCCGAGACGATTCTTTATCTTGGTGGGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGGAATAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGGGAAGTGCAtttt                                                                                                                                                     >  2:114397/1‑250 (MQ=255)
                                                                                                 ccTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCt                                                                                                                                                  >  1:179640/1‑251 (MQ=255)
                                                                                                    ttttATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCCCCAGCGTGGTTGGTGCCGATCAGGATGATGCCAAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTg                                                                                                                                                 <  1:377138/249‑1 (MQ=255)
                                                                                                        aTCAGCGCTGGGTTAACGTCCGGTTCCCGATTGGCGCGATGGACCGCACACTCTATCGCGACCCTTTCACCAGCGTTGTTGGTCCCGATCAGAATATGGCCGAAATGCTGGCGGAAGAGTAACGTAAGTTTTCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGTTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGc                                                                                                                                            <  1:241771/250‑1 (MQ=255)
                                                                                                        aTCAACGCTGGGCTAACGCCCCGTCCCCGATTGCCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGATTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGc                                                                                                                                            <  1:407381/250‑1 (MQ=255)
                                                                                                                 gggCTACCGCCCCGTTCCCGATTGTCGCGCTGGCCCGCGCCCTCGATCGTGAACACTTCACCGGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTAt                                                                                                                                  <  1:224790/251‑1 (MQ=255)
                                                                                                                 gggCTACCGCCCCGTCCCCGATGGTCGCGCTGGACCGCGCCCTCGCTCGTAAACATTTCACCAGCGGGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTAt                                                                                                                                  <  1:6456/251‑1 (MQ=255)
                                                                                                                     tAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGCTGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGGTTCCCGCCGAGCCGAAGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGACAGATGATCCGCGCGCAGTGCATTTCCTGTATGCCACCCGCTATGAGc                                                                                                                              >  2:29542/1‑251 (MQ=255)
                                                                                                                       aCGACCCTTTCCCTATTGCCGCGCGGGACCGCGCCCGCGATCGTGAACACTTCACCAGGGTGGTTGGGGCCGATCAGGAGGGTGCCGAAAGGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGAGGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGAGCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCgg                                                                                                                            <  1:534192/251‑1 (MQ=255)
                                                                                                                         gACCCGTTCCCGATTCTCGCGCTGGACCGCGCCCTAGATCGTGAACACTTCACCAGCGTGGTTGGAGCCGATCAGGATGATGCCGAACTGCTGGCGGAAGAGTTACGTAAATTTCCCGCCGAGACGGATCTTTATCCTGGTGCGCTACCGGCGCTTTCTGTAAGCTTCCTGCGTGATCAAGGTTTCCGTAATGCCTGGACAGATGATACGCGCGCAGTACATTTCCGTGTTGCCAACTGGTATGaagggg                                                                                                                           >  2:231773/1‑245 (MQ=255)
                                                                                                                            ccGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTAACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAggc                                                                                                                       >  1:431694/1‑251 (MQ=255)
                                                                                                                              gTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAggcgg                                                                                                                     >  1:179616/1‑251 (MQ=255)
                                                                                                                                         tCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCtt                                                                                                                                                                                                               >  1:368365/1‑151 (MQ=255)
                                                                                                                                           gcgcTGGACCGCCCCCTCGATCGTAAACACTTCACCAGCGTGGTTGGTGCCAAGCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCATCAGCTATGAGCGGGAGGCGGCTGCCCAGTTAtt                                                                                                        <  1:450810/251‑1 (MQ=255)
                                                                                                                                             gcTGGACCGCGCCCTCGAGCGTGAACATTTCCCCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCg                                                                                                      <  1:569686/251‑1 (MQ=255)
                                                                                                                                               tGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGa                                                                                                     >  1:15578/1‑250 (MQ=255)
                                                                                                                                                ggACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGaa                                                                                                    >  1:471435/1‑250 (MQ=255)
                                                                                                                                                 gACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGAGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCttt                                                                                                                                                                                                                         >  1:449182/1‑133 (MQ=255)
                                                                                                                                                                          tCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCg                                                                         >  1:224719/1‑251 (MQ=255)
                                                                                                                                                                            aCCAGCGTGGTGGGTGCCGAACAGGAGGTGGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCGCa                                                                       <  1:552074/251‑1 (MQ=255)
                                                                                                                                                                              cAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTg                                                                                                                                                 >  1:127213/1‑175 (MQ=255)
                                                                                                                                                                              cAGCGTGGTGGGGGCCGATCAGGAGGATGCCGAAATGCTGGCGGAAGAGTTCCGTAAGTTCCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTg                                                                                                                                                 <  2:127213/175‑1 (MQ=255)
                                                                                                                                                                               aGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTACGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCGCAGgc                                                                    >  1:77167/1‑251 (MQ=255)
                                                                                                                                                                                                  aGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCGCAGGCGCTGTTCACAACGTCGttt                                                >  1:298463/1‑251 (MQ=255)
                                                                                                                                                                                                              gAAATGCTGGCGGAAGAGTAACGTAAGTTCCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCGCAGGCGCTGTTCACAACGTCGTTTGCGTTGTTGCa                                     <  1:114397/251‑1 (MQ=255)
                                                                                                                                                                                                                     tGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGACGGCTGCCCAGTTATTCGAACAATGGCTGGAAACGCATCCGATGCCGCAGGCGCTGTTCACAACGTCGTTTGCGTTGTTGCAAGGGGTg                              >  2:265932/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaGAGTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACAGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCGTAGGCGCTGTTCACAACGTCGTTTGCGTTGTTGCAAGGAGTgatgggt                        >  1:364283/1‑249 (MQ=255)
                                                                                                                                                                                                                              agTTACGTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCCTCCGATGCCGCAGGCGCTGTTCACAACGCCGTTTGGGGTGTTGCAAGGAGTgcgggaggt                      >  2:7747/1‑242 (MQ=255)
                                                                                                                                                                                                                                 tACGGAAGTCTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGTCAACAGCTATGATCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCGCAGGCGCTGTACACAACGTCGTTTGCGTTGTTGCAAGGAGTGATGGATGTCAc                   <  1:40705/250‑1 (MQ=255)
                                                                                                                                                                                                                                   cGTAAGTTTACAGCCGCGACGATGATTTATCTTGGTGCGCTACCGGAGCTTCCTGTCAGCTTCAAGCGAGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCAACGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCACAGTTATTCGAAAAATAGCTGTACACGAATCCGCTGCCGCAAGCGCTGTTCCAACCGTCGTTGTCGTTGTTCCAAGGCGTGATGGATGTccagct                >  2:556612/1‑246 (MQ=255)
                                                                                                                                                                                                                                    gTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAACCGCATCCGATGCCGCAGGCGCTGTTCACAACGTCGTTTGCGTTGTTGCACGTAGTGATGGATGTCACGCTg               >  1:337473/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gTAAGTTTCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATACGATGACGCAGGCGCTGTTCACAACGTCGGTTGCGTTGTTGCAAGGTGTGATGGATGTAACGCTg               >  2:398148/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gttccccGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCGCAGGCGCTGTTCACAACGTCGTTTGCGTTGTTGCAAGGAGTGATGGATGTCACGCTGCGTc           <  2:322066/247‑1 (MQ=255)
                                                                                                                                                                                                                                        gtctCCCGCCGAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCGCAGGCGCTGTTCACAACGTCGTTTGCGTTGTTGCAAGGAGTGATGGATGTCACGCTGCGTc           <  2:154429/248‑1 (MQ=255)
                                                                                                                                                                                                                                             ccgccgAGACGATGCTTTATCTTGGTGCGCTTCCGGATCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCGCAGGCGCTGTTCACAACGTCGTTTGCGTTGTTGCAAGGAGTGATGGATGTCACGCTGCGTCGCGAc      <  1:24672/251‑1 (MQ=255)
                                                                                                                                                                                                                                              cgccgAGACGATGCTTTATCTTGGTGCGCTACCGGCGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCGCAGGCGCTGTTCACAACGTCGTTTGAGTTGTTGCAAGGAGTGGTGGATGTCACGCTGCGGCGCGACg     >  2:371097/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gccgAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCa                                                                                                              >  1:342156/1‑146 (MQ=255)
                                                                                                                                                                                                                                                 cgAGACGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCGCAGGCGCTGTTCACAACGTCGTTTGCGTTGTTGCAAGGAGTGATGGATGTCACGCTGCGTCGCGACGGc   >  1:211489/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   agaCGATGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCCGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCGCAGGCGCTGTTCACAACGTCGTTTGCGTTGTTGCAAGGAGTGGTGGATGTCACGCTGCGTCGCGACGGcg  >  2:520322/1‑248 (MQ=255)
                                                                                                                                                                                                                                                   agaCGATGCGTTATCTTGGTGCGCGACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCGCAGGCGCTGTTCACAACGTCGTTTGCGTTGTTGCAAGGAGTGATGGATGTCACGCTGCGTCGCGACgg    <  1:520322/247‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
AGAAGATCAGCCAGACAACGAAATGCGGTGCATTGAGCACCTTTTACAGCGTCAGGTTGATGCCATTATTGTTTCGACGTCGTTGCCTCCTGAGCATCCTTTTTATCAACGCTGGGCTAACGACCCGTTCCCGATTGTCGCGCTGGACCGCGCCCTCGATCGTGAACACTTCACCAGCGTGGTTGGTGCCGATCAGGATGATGCCGAAATGCTGGCGGAAGAGTTACGTAAGTTTCCCGCCGAGACGGTGCTTTATCTTGGTGCGCTACCGGAGCTTTCTGTCAGCTTCCTGCGTGAACAAGGTTTCCGTACTGCCTGGAAAGATGATCCGCGCGAAGTGCATTTCCTGTATGCCAACAGCTATGAGCGGGAGGCGGCTGCCCAGTTATTCGAAAAATGGCTGGAAACGCATCCGATGCCGCAGGCGCTGTTCACAACGTCGTTTGCGTTGTTGCAAGGAGTGATGGATGTCACGCTGCGTCGCGACGGCA  >  NC_000913/88324‑88814

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