Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,684,443 C→T K392K (AAG→AAA pdeK ← c‑di‑GMP phosphodiesterase PdeK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,684,4430CT97.9% 161.2 / ‑3.7 48K392K (AAG→AAApdeKc‑di‑GMP phosphodiesterase PdeK
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (30/17);  total (31/17)
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.00e-01

CCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGAT  >  NC_000913/3684210‑3684685
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
ccACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGc                                                                                                                                                                                                                                   <  1:14295/251‑1 (MQ=255)
      cACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAAATGAATCTGATTTTTGCCTTTATGGCGCGCGGTaa                                                                                                                                                                                                                              >  2:296892/1‑250 (MQ=255)
       aCTCAATGCGATAGATTAAGCCATCCGGCAGGTCCAAACTGCCATCCGGTTACTGGATACGCAGTAACACTTCCGCCCTAACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTCAAt                                                                                                                                                                                                                            >  2:145737/1‑251 (MQ=255)
               cgatcgatTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGa                                                                                                                                                                                                                    <  1:59772/251‑1 (MQ=255)
                 atcgatTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAACGCCGCAGaa                                                                                                                                                                                                                   >  1:201832/1‑250 (MQ=255)
                          gCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGAc                                                                                                                                                                                                         >  2:231526/1‑251 (MQ=255)
                            cATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTg                                                                                                                                                                                                       >  2:275365/1‑251 (MQ=255)
                             aaccGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGCCCAGTTTACCGCTGGTCATCTCGACCTTTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTGCCGTCATCCGCTGCTGTGCGGCTTCCATGTGCTGCGGATCAAAGGACTGAATCTGGTTTTTGCCTTTAGGGCAAGTGGTAAATGCCGTAGAGATTGCGCGACTGt                                                                                                                                                                                                      <  2:143563/249‑1 (MQ=255)
                                    aGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTg                                                                                                                                                                                               <  1:140880/251‑1 (MQ=255)
                                                 cATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTg                                                                                                                                                                                     >  1:9212/1‑248 (MQ=255)
                                                       gTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCg                                                                                                                                                                            >  1:102449/1‑251 (MQ=255)
                                                         tGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGt                                                                                                                                                                           <  1:19943/250‑1 (MQ=255)
                                                          gCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGt                                                                                                                                                                           >  1:225787/1‑249 (MQ=255)
                                                          gCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGGGAGATCGCCGt                                                                                                                                                                           >  2:106815/1‑249 (MQ=255)
                                                                                cGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTAc                                                                                                                                                   <  1:161240/251‑1 (MQ=255)
                                                                                  cACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTAc                                                                                                                                                   >  1:145518/1‑249 (MQ=255)
                                                                                  cACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTAc                                                                                                                                                   >  1:249/1‑249 (MQ=255)
                                                                                        ccAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAATTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGcgg                                                                                                                                            >  2:66396/1‑250 (MQ=255)
                                                                                         cAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGcggac                                                                                                                                          >  1:303475/1‑251 (MQ=255)
                                                                                             ttACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGcggacgga                                                                                                                                       <  1:100775/250‑1 (MQ=255)
                                                                                                     ggTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAAt                                                                                                                                <  2:52540/249‑1 (MQ=255)
                                                                                                       tCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGt                                                                                                                            >  1:130293/1‑251 (MQ=255)
                                                                                                             cGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATc                                                                                                                      >  2:75601/1‑251 (MQ=255)
                                                                                                             cGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATc                                                                                                                      <  1:163524/251‑1 (MQ=255)
                                                                                                                   gtgGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATCCGTTCAATCGGCAgg                                                                                                                >  2:201602/1‑251 (MQ=255)
                                                                                                                       cTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGcgct                                                                                                            >  1:105137/1‑251 (MQ=255)
                                                                                                                                 aTAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGAt                                                                                                  >  2:278035/1‑251 (MQ=255)
                                                                                                                                                        tGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAgg                                                                           >  2:141437/1‑251 (MQ=255)
                                                                                                                                                           aTTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTg                                                                         >  1:53455/1‑250 (MQ=255)
                                                                                                                                                            ttAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCTCTCGCTCATGATAGTGTGCCCTTGCTGCCCTAAGGGGAt                                                                       >  2:103515/1‑251 (MQ=255)
                                                                                                                                                            cgtggctcttccgatctttCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGAt                                                                       <  1:124224/234‑1 (MQ=255)
                                                                                                                                                             tAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGAtt                                                                      >  1:15905/1‑251 (MQ=255)
                                                                                                                                                                      cGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCAc                                                             >  1:174050/1‑251 (MQ=255)
                                                                                                                                                                         tctctTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGt                                                          >  1:17180/1‑251 (MQ=255)
                                                                                                                                                                          ctctTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGtt                                                         >  1:228371/1‑251 (MQ=255)
                                                                                                                                                                             ttCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGAt                                                                       >  2:124224/1‑234 (MQ=255)
                                                                                                                                                                                       cGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGtt                                            <  1:275365/251‑1 (MQ=255)
                                                                                                                                                                                        gCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTg                                           >  1:265962/1‑251 (MQ=255)
                                                                                                                                                                                                  ggcTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAg                                 <  1:90436/251‑1 (MQ=255)
                                                                                                                                                                                                             gctgcGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTc                          <  1:273822/247‑1 (MQ=255)
                                                                                                                                                                                                             gctgcGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTGCCTGTACACCGTTGGCAATGACAGCAAAGTCATa                       <  2:46922/250‑1 (MQ=255)
                                                                                                                                                                                                                gcGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAAcc                    <  2:265962/250‑1 (MQ=255)
                                                                                                                                                                                                                 cGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGc                  >  1:104433/1‑251 (MQ=255)
                                                                                                                                                                                                                    aTCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTa                >  2:134552/1‑250 (MQ=255)
                                                                                                                                                                                                                       aaaGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCt            <  1:35133/251‑1 (MQ=255)
                                                                                                                                                                                                                          gAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTgcg         <  2:53455/251‑1 (MQ=255)
                                                                                                                                                                                                                           aaCTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTgcgc        >  2:255601/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aTCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGAt  >  1:123918/1‑251 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                  
CCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGAT  >  NC_000913/3684210‑3684685

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