Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 47,099 G→A Q431Q (CAG→CAA yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091347,0990GA96.4% 160.5 / ‑4.3 57Q431Q (CAG→CAAyaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (32/22);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGTCGTATTGGCACCATTGTTTCGACCTGGGCACTACCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAGACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTCCCATGCGAATAAACGGATGCTTGAACAGGCAAGGACGCTGGAAGGCGTCGAA  >  NC_000913/46853‑47338
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
cGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAAccc                                                                                                                                                                                                                                             >  2:115025/1‑250 (MQ=255)
 gggCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCCCTGGCGCAAAcc                                                                                                                                                                                                                                              >  1:670789/1‑249 (MQ=255)
     cGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGCTTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGGTGTCACTGGCGCAAACCAGCact                                                                                                                                                                                                                                        >  2:268390/1‑249 (MQ=255)
                  gTGGCGCTATCATAATGATCTTGTCCCGACAGTTATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGCCCTGGGCACTGCTGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGCCGTTCcgcg                                                                                                                                                                                                                           <  2:243067/251‑1 (MQ=255)
                   tGGCTCTATCCTAATGAACTCTTCCCGACAGAGATCCGAGCCTCTGCCGTGGGCGGGATTATGTCCTTAAGCCGGATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGAATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTTCCCCGGAGACTCGAAGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCgg                                                                                                                                                                                                                          >  1:120723/1‑251 (MQ=255)
                       tctATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAatat                                                                                                                                                                                                                                      >  1:104138/1‑235 (MQ=255)
                       tctATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAatat                                                                                                                                                                                                                                      <  2:104138/235‑1 (MQ=255)
                          aTCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCagag                                                                                                                                                                                                                    <  1:22015/250‑1 (MQ=255)
                          aTCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAgaga                                                                                                                                                                                                                   <  1:445783/251‑1 (MQ=255)
                               aataatatcttcccaccagagttCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGGATTGGCACCATTGTGTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATggg                                                                                                                                                                                                              <  2:520398/229‑1 (MQ=255)
                               aaTGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTTGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATgg                                                                                                                                                                                                               <  1:614875/250‑1 (MQ=255)
                                 tGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATg                                                                                                                                                                                                                <  1:633316/247‑1 (MQ=255)
                                 tGAACTCTTCCCAAAAGATATCCGCGCCTCTGCCGTGGGCGCGGTTATGTCCTTCAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGTTTTCCGTAACGTTTGACACGGAGAATCGAGGGATGTAACTGGCGCAAACCAGCACTATGACGACGCGCGGGCAGGGGATgggg                                                                                                                                                                                                             >  2:633316/1‑250 (MQ=255)
                                           ccGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCCGCCATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTc                                                                                                                                                                                                   >  2:506845/1‑250 (MQ=255)
                                                    atCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTtctct                                                                                                                                                                                         >  1:487794/1‑251 (MQ=255)
                                                                   gTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAat                                                                                                                                                                          >  1:239336/1‑251 (MQ=255)
                                                                   gTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAat                                                                                                                                                                          >  2:483332/1‑251 (MQ=255)
                                                                     gggCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAata                                                                                                                                                                         <  1:569073/250‑1 (MQ=255)
                                                                      ggCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTCTCAATAATTACGGAATCAGTAACACGATGCTAATGTGGGCGGGGATCTCGCAGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGTGTCTCGAGGGATGTCACTGGCGAAACCCAGAATTATGAACATCCGAGGGCTGGGATTGGGGGAATTTTTTCAGGTTTCTCCCTTTTCTTAAtacataat                                                                                                                                                                       >  2:619842/1‑244 (MQ=255)
                                                                                      tAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCcg                                                                                                                                                        <  1:483332/250‑1 (MQ=255)
                                                                                       aagccttATTGGCCCCTTTGTTTCGACCTGGGCACTGACGTTCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCcgc                                                                                                                                                       <  2:209545/244‑1 (MQ=255)
                                                                                             tATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACAGTATCAGTAACACGATGCTAATGGGGGAGGGTATCTCGCTGTTTGACTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACACGAGGGATGTCACTTGCGCAAACCAGCAATATGACGTTCCGCGGGCGGAGAAGTGGGTAAATTGTTCAGATTTCTCTCATTTCTGAATCAATATTATTGACAATAAGCCGCGTGAat                                                                                                                                                 >  2:197842/1‑250 (MQ=255)
                                                                                               ttGGCACCACTGTGTCGACCTGGGCTCTGCCGATCTTTATCAATAATGACTGTATCAGGAACACGATGCTAATGGGGGTGGGTATCTCGCGGTTTGTCTTGTTGATTTCCGTAGCGGTTGCCCCGGAGTCTCGAGGGATTTCACTGGCGCATACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAAGTGTTCATATTGCTCTCTTTTCTGAATCAATATTATTGACTACAAGCCGCTTGAatag                                                                                                                                               <  2:287504/250‑2 (MQ=255)
                                                                                                  gCACCATTGTTTCGACCTGGTCACTGCCGATCTTTATCAATAATTACTGTACCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGa                                                                                                                                           <  1:368067/251‑1 (MQ=255)
                                                                                                    cccATTGCTTCGACCTGGTCACTGCCGATCTTTATCAATAAGTACGGTAACAGTAACACGATGCGAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGTGCAGAGAATGGTGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACt                                                                                                                                         <  2:245293/250‑1 (MQ=255)
                                                                                                    aCCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATAAGACt                                                                                                                                         >  1:145634/1‑251 (MQ=255)
                                                                                                       aTTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTGCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTaca                                                                                                                                      <  1:89972/251‑1 (MQ=255)
                                                                                                          gTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACtt                                                                                                                                   >  1:661040/1‑251 (MQ=255)
                                                                                                          gTTTCGACCTGAGCACTGACGATCTTTATCACTAATTACGGTATCAGTACCACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGATGATTTCCGTAACGTTTGCCCCGGAGACTCGAGGGATGTCATTGGAGCAAACAAGCAATATGACGATCCGCGGGCAGAGAATGGGGTCAATTGTTCAGATTTCTCTATTTTCTGAATCAATTTTTTTGACTATAGGCCGCGTCACTATATTACtcaactt                                                                                                                                   >  2:276602/1‑245 (MQ=255)
                                                                                                              cGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTtgtg                                                                                                                               <  1:506845/251‑1 (MQ=255)
                                                                                                               gACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTAAAAGCCGCGTGAATATATGACTACACTTTGTgg                                                                                                                              >  2:554708/1‑251 (MQ=255)
                                                                                                                        actgccGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCCTGTTGATTTCCGTAGCGTTTGCCCCGGCGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATACCGATCCGCGGGCAGAGGCTGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAAACTTAATGCCTATAAGCCGCGAGAATATATTACTACACTTTGTGGGAAAACaa                                                                                                                      >  2:196233/5‑250 (MQ=255)
                                                                                                                               aTCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGAAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATTACTACACTTTGTGGGAAAACAAAGGCgaa                                                                                                               >  2:130377/1‑248 (MQ=255)
                                                                                                                                tCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAAt                                                                                                             >  1:478128/1‑251 (MQ=255)
                                                                                                                                    tATCAATAATGACGGTATCAGTAACACGATGCTAATGGGGGCGGACAAATCGCAGTTTGGATTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCAATGGCGCAAAACAGCCATATGACGGTCCGCGGGCAGAGAATAGGGTAAATTGTTCAGCTTTCTCTTTTTTCTGAATCAATATTATTGACTTTAAGACGCGGGGATATATTACAACACTTTTTGGGGAAACAAAgaagtacacagg                                                                                                         >  2:491737/1‑240 (MQ=255)
                                                                                                                                       cAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCgg                                                                                                      <  1:115025/251‑1 (MQ=255)
                                                                                                                                         ataatTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGc                                                                                                    >  1:537021/1‑251 (MQ=255)
                                                                                                                                         ataatTACGGTATCAGTAACACGATGCTAATGCGGGCGGGTATCTCGGTGTTTGGCTTGTTAATTTCCGTACCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCCGCAATATGACGATCCGCGGGGAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTAATGACTATAAGCTGCGTGAATATATGAATACACTTTGTGGGAAATAGACGGCGTAATCAcgcaggg                                                                                                    >  2:329716/1‑247 (MQ=255)
                                                                                                                                          taatTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGAAGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCt                                                                                                   >  1:631039/1‑251 (MQ=255)
                                                                                                                                           aatTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCt                                                                                                   <  1:34552/250‑1 (MQ=255)
                                                                                                                                               actgtagcAGTAACACGATGCTAATGGGGGCGGGGATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGTGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGTCTACGCTTTGTGGGAAAACAAAGGCGTAATCACGCGGGTTACCt                                                                                               <  1:570555/243‑1 (MQ=255)
                                                                                                                                                 ggTATCAGTCACACGATGCTAATGTGGGCGGGAATCTCGCGCTTTGGATTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCCAGGGATGTCACTGGCGCCAACCAGCACTATGACGATCCCCGGGCAGAGAATTGGGCAAATTTTTCCAATTTCTCTCTTTTCTGCATCCATATTATTGACTATAAGCCGCGCGAAAATATCACTACACTTTGTGGGAAAAAACAGGCGTAATCCCACGGGCTACCTATg                                                                                            >  2:641837/1‑251 (MQ=255)
                                                                                                                                                              cGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGGAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAAttt                                                                               >  2:103982/1‑251 (MQ=255)
                                                                                                                                                                    tAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAGACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATgcgc                                                                         >  1:350012/1‑251 (MQ=255)
                                                                                                                                                                     aaTGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGGGGGATGTCAATGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGAAATCACGCGGGCTACCTATGATTCTTATAATTTATgcgc                                                                         >  1:448900/1‑250 (MQ=255)
                                                                                                                                                                        gggggCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGAGCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTAGGCGCATcc                                                                     <  2:304951/251‑1 (MQ=255)
                                                                                                                                                                            gCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAGACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTAt                                                                 <  2:350012/251‑1 (MQ=255)
                                                                                                                                                                               ggTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGAGTCTTATAATTTATGCGCATCCGTATCCg                                                              <  2:661040/251‑1 (MQ=255)
                                                                                                                                                                                                     tgATTTCCGTAGCGTTTGCCCCGGAGACTCGCGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCCCGGGCAGAGCATCGGGCAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGCATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCCTCATTCCCATGCGAATAAAc                                        >  1:613018/1‑251 (MQ=255)
                                                                                                                                                                                                            cGTAGCGGTTGCCCCGGAGAATCGAGGGATGTCACTCGCGCAAACCAGCAATATGACGATCCGCGGGCCGAGAATGGGGTAAAATGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGAGAAAACACAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCAACCGTATCAGCCTCATTCCCATGCGAATAAACGGATGCt                                 >  2:642828/1‑251 (MQ=255)
                                                                                                                                                                                                                      gCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGAGCATCCGTATCCGCATCATTCCCATGCGAATAAACGGATGCTTGAACAGGc                        >  1:369026/1‑250 (MQ=255)
                                                                                                                                                                                                                             agaCTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTCCCATGCGAATAAACGGATGCTTg                               <  1:181162/236‑1 (MQ=255)
                                                                                                                                                                                                                             agaCTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTCCCATGCGAATAAACGGATGCTTg                               >  2:181162/1‑236 (MQ=255)
                                                                                                                                                                                                                               aCTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTAAACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGATACCTATGATTCTTATACTTTATGCGCATCCGTATCCGCATCATTCCCATGCGAATAAGCGGATGCTTGAACAGGCAAGTACGCTg              >  2:451605/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTCCCATGCGAATAAACGGATGCTTGAACAGGCAAGGACGCTGGAAgg         <  1:554708/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTCCCATGCGAATAAACGGATGCTTGAACAGGCAAGGACGCTGGAAGGCGTCGaa  >  1:425912/1‑250 (MQ=255)
                                                                                                                                                                                                                                            cACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTCCCATGCGAATAAACGGATGCTTGAACAGGCAAGGACGCTGGAAGGCGTCGaa  <  1:103982/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               
CGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGTCGTATTGGCACCATTGTTTCGACCTGGGCACTACCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAGACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAAAACAAAGGCGTAATCACGCGGGCTACCTATGATTCTTATAATTTATGCGCATCCGTATCCGCATCATTCCCATGCGAATAAACGGATGCTTGAACAGGCAAGGACGCTGGAAGGCGTCGAA  >  NC_000913/46853‑47338

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