Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,321,488 G→A A70V (GCC→GTC)  phnG ← carbon‑phosphorus lyase core complex subunit PhnG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,321,4880GA95.5% 187.7 / NA 66A70V (GCC→GTC) phnGcarbon‑phosphorus lyase core complex subunit PhnG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (29/34);  minor base C (1/2);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGGCGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGC  >  NC_000913/4321244‑4321729
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
tCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCggg                                                                                                                                                                                                                                              >  1:59432/1‑251 (MQ=255)
   cgcttccgatcggcgAACCATCGTAAAGAAGGCGACCCGGCTGGCGTTCACTTCGGCCGGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGGTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTGATCACGCCCCTGCGCCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCa                                                                                                                                                                                                                                           <  1:17187/239‑1 (MQ=255)
         gTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTgg                                                                                                                                                                                                                                     <  1:143995/251‑1 (MQ=255)
         gTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTg                                                                                                                                                                                                                                      <  1:223871/250‑1 (MQ=255)
         gTCTCCGCGAACCATCGGAAAGAAGTCGACCCGGCGGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTGATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTgg                                                                                                                                                                                                                                     <  1:118490/251‑1 (MQ=255)
               gcgAACCATAGTACAGCAGTCGAACCGGCTGCCGTTCACTTCGGCCTCGCGTGCGACAAAGCGTGCCATACGCTCAGCGTCCAGCGGGGCACTAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGACAGACGCACGGCGACGCGGGTCa                                                                                                                                                                                                                                           >  2:17187/1‑239 (MQ=255)
                cgAACCCTCGTAAAGAAGTCGACCCGGCGGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGGCAGTCAGGCGCACGGCGACGCGGGTCATCGTGGCGTCGcc                                                                                                                                                                                                                              <  2:186770/251‑1 (MQ=255)
                 gAACCATCGTATAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTTCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCg                                                                                                                                                                                                                             <  1:105410/251‑1 (MQ=255)
                    ccATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCa                                                                                                                                                                                                                          <  1:159553/251‑1 (MQ=255)
                                        cGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTgcc                                                                                                                                                                                                      <  2:78105/251‑1 (MQ=255)
                                            tGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGccc                                                                                                                                                                                                  >  1:66304/1‑251 (MQ=255)
                                             ggCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATTCGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAATTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCTCCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGGCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCa                                                                                                                                                                                                 <  1:235855/251‑1 (MQ=255)
                                                  tCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGAGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCCCGGCGACGCGGGTCAGCGTGGCGTCGCCGGCCAACAACCGTTCGCCGGTGCCGCCCACCcgc                                                                                                                                                                                            >  2:132541/1‑251 (MQ=255)
                                                   cACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGGTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCcgc                                                                                                                                                                                            <  1:172245/250‑1 (MQ=255)
                                                       gcGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGGCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGGCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTATCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGTGGGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGcc                                                                                                                                                                                         <  1:126327/248‑1 (MQ=255)
                                                        cGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGa                                                                                                                                                                                      <  2:28319/251‑1 (MQ=255)
                                                         ggCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGaa                                                                                                                                                                                     <  1:128640/251‑1 (MQ=255)
                                                            cTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCt                                                                                                                                                                                  >  2:43408/1‑251 (MQ=255)
                                                                     ggCAATGCGTGCCATACGGCCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGc                                                                                                                                                                         <  1:126897/251‑1 (MQ=255)
                                                                        aaggCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGAGTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAg                                                                                                                                                                      <  2:62382/248‑1 (MQ=255)
                                                                        aaTGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGcc                                                                                                                                                                        >  2:135214/1‑249 (MQ=255)
                                                                             gTGCCATATGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCGATCAGCGCGTAGCGTTCGTCTTGCTGTTTATCACGCCCCGGCTCCCAGCTGGAGCCGAGCGTGCCGTCAGTCAGGCGCCCGGCGACGCGGGGCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCTATCCGCGCCTGAATCTGTACCAGGCCAGTTTCa                                                                                                                                                                 <  2:62594/251‑1 (MQ=255)
                                                                                       gTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCgg                                                                                                                                                        >  1:238316/1‑250 (MQ=255)
                                                                                       gTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGACGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCgg                                                                                                                                                        >  1:182674/1‑250 (MQ=255)
                                                                                        tCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATCCGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCTCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCGAAACAACGTTCGCAGGGGCCGCCCATCCGCTCCTGACTCTGTACCAGGCCCGTTTAAgcgcggcggcc                                                                                                                                                      >  2:151336/1‑243 (MQ=255)
                                                                                             gTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGCTCAccc                                                                                                                                                 >  2:90837/1‑250 (MQ=255)
                                                                                                 aGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATa                                                                                                                                             <  1:244100/251‑1 (MQ=255)
                                                                                                    ggggCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTc                                                                                                                                          <  2:209288/251‑1 (MQ=255)
                                                                                                            aagggttgtggataggtgatGAAAGTGGCGAGAATGCGTCAGCAGCGCGGCAAGCAGCGCGCAGCTTTCGGCGTGCGGTTTATCACGCCCCTGCACCCATCTTTAGCCGAGCGGGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCTCCCATCCGCGCCTGAGTCTTTGCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGa                                                                                                                                  <  1:160019/232‑1 (MQ=255)
                                                                                                                   tCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTa                                                                                                                            <  1:209395/250‑1 (MQ=255)
                                                                                                                       aTAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCCGTCAGGCGCACGGCGCCGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGCATCTGTACCAGGCCAGGTTCAGCGGCGCGGCTCACCTCATCGACGGCGGTGGCGTTTAGcgcg                                                                                                                       >  2:61721/1‑251 (MQ=255)
                                                                                                                           gTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCTGTGACGTTTAGCGCGTTCa                                                                                                                   <  2:262246/251‑1 (MQ=255)
                                                                                                                           gTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCa                                                                                                                   <  1:55414/251‑1 (MQ=255)
                                                                                                                                      ggCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGAGCTGTTTATCACGCCCCTGCCCCCCGCTGTAGCAGAGCGTGCCGTCAGCCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGCAGCCGCCCATCCGCGCCTGCATCTGTACCAGGCCAGATTACGCGGCGCGGCCCACCTCATAGTCGGCGGTGACGTTTAGCGAGTTCAGGCGTGCTGcc                                                                                                        >  1:124538/1‑251 (MQ=255)
                                                                                                                                        cGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGcc                                                                                                        >  2:208241/1‑249 (MQ=255)
                                                                                                                                         gagaTTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGGTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGTTGTAGTCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCGGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCGGTACCAGGCCAGTTTCAGCGGCGCGGATCACCGCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGt                                                                                                     <  1:266361/251‑1 (MQ=255)
                                                                                                                                          agaTTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCCGGCTTGCTGCCCGtt                                                                                                    >  2:38912/1‑251 (MQ=255)
                                                                                                                                             ttGCTGCATAAGCGCGTCAATCAGCGCGCAGCGTTCAGCGTGCTGTTTATCACGCCACTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATACTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGACAGTTc                                                                                                   >  2:188838/1‑249 (MQ=255)
                                                                                                                                              tgctgcATCAGCGAGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGGGACGTTTAGCGCGTTCAGGCGTGCTGCCCGTTCAGc                                                                                                >  2:267196/1‑251 (MQ=255)
                                                                                                                                                     tCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTgg                                                                                         >  1:151797/1‑251 (MQ=255)
                                                                                                                                                           cgTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTCTCACGCCCCTGCACCCAGATGTAGCCGCGCGTGCCGTCCGTCCGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCCTCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGGTCACCTCCTAGTCGGCGGCGACGTTTAGGGCATTCCGGCGGGCTGCCAGCTAAGCCGGTTGGCTGTgc                                                                                   >  2:123528/1‑251 (MQ=255)
                                                                                                                                                            gTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTgcg                                                                                  >  1:133951/1‑251 (MQ=255)
                                                                                                                                                                         cAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCCCGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAACAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGCATCTGTACCAGGCCAGTTTCAGCGGCGCGGCTCACCTCATAGTCGGCGGTGACGTTTCGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTGTGCTGTGCGCCAGCACGGCCAt                                                                     >  2:133354/1‑251 (MQ=255)
                                                                                                                                                                             gTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATcc                                                                   >  1:40972/1‑249 (MQ=255)
                                                                                                                                                                                cGGCGTGCTGGTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGc                                                              <  2:213290/251‑1 (MQ=255)
                                                                                                                                                                                  gCGTGCTGTTTATCACGCCCCAGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGGGCGCCAGCACGGACAGCCAgctctg                                                            >  2:53217/1‑246 (MQ=255)
                                                                                                                                                                                   cGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACCTCCAGTGCTgg                                                           >  1:170651/1‑251 (MQ=255)
                                                                                                                                                                                   cGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTgg                                                           >  1:102293/1‑251 (MQ=255)
                                                                                                                                                                                           gtATCACGCCCCGGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTcg                                                    <  2:264435/249‑1 (MQ=255)
                                                                                                                                                                                             aTCACGCCCCTGCACCCAGCTGTCGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGGGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGGCGCgg                                                 >  2:21521/1‑251 (MQ=255)
                                                                                                                                                                                              ccACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGt                                                <  1:162182/250‑1 (MQ=255)
                                                                                                                                                                                                 cGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATc                                             <  2:31731/251‑1 (MQ=255)
                                                                                                                                                                                                  gCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGCCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTGCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTTGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCt                                            <  1:175826/251‑1 (MQ=255)
                                                                                                                                                                                                               gCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGt                              <  1:90837/251‑1 (MQ=255)
                                                                                                                                                                                                                  gTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTg                             <  1:62603/249‑1 (MQ=255)
                                                                                                                                                                                                                  gTAGCCGAGCGGGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCt                           <  2:59432/251‑1 (MQ=255)
                                                                                                                                                                                                                   gagcggtgCGTGCCGTCAGGCAGGCGCACGGCGACGCGGGTCAGCGGGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTTCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTGGCGGTATCTGCGTGCAGTCAGTGCTc                          <  2:42072/244‑1 (MQ=255)
                                                                                                                                                                                                                       cGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTGTCTGCGTGCATTCAGTGATCCAt                       >  1:132416/1‑250 (MQ=255)
                                                                                                                                                                                                                            gTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGGGCTGGCGGGGCCCGGTATCTGGGTGGATTCCGGGCTcccttgtga                  >  2:8513/1‑243 (MQ=255)
                                                                                                                                                                                                                                 gtcggtccggCGCACGGCGCCGGGGGCCACCGTGGCATCGCCTGCTAAAAACCGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGTCCATTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGTTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCaa            <  2:221193/243‑1 (MQ=255)
                                                                                                                                                                                                                                 gtcagtcaGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCCTTCAGTGCCCCATAGTGAATTCa             >  1:139255/1‑250 (MQ=255)
                                                                                                                                                                                                                                 ggcaggcagGCGCACGGCGACGCGGGTCAGCGTGGCGGCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCaa            <  2:198222/245‑1 (MQ=255)
                                                                                                                                                                                                                                   cagtcagGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGGCGCGGGATCTGCGTGCGTTCTAGGCCCCATAGTGGATTCAAt           >  2:29207/1‑250 (MQ=255)
                                                                                                                                                                                                                                       cagGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTGTCTGCGGGCATTCCGTGCTCCCTAGAGAATTTAAACATGt      >  1:50434/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gcgcACGGCGCCGCTGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGGTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCgg   <  2:1383/251‑1 (MQ=255)
                                                                                                                                                                                                                                           cgcACGGCGACGCGGGTCAGCGGGGCGGCGCCGGAAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGGACCAGGCCAGTTTCAGCGGCGCGGATCACCTCAGAGTCGGCGGTGACGTTTAGCGCGGTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGgc  <  1:106677/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
TCATGCGTTGTCTCCGCGAACCATCGTAAAGAAGTCGACCCGGCTGGCGTTCACTTCGGCCTGGCGTGCGGCAATGCGTGCCATACGGTCAGCGTCCAGCGGGGCAATAAGGGTTTCTGATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGGCGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGC  >  NC_000913/4321244‑4321729

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