Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,345,404 G→A A52V (GCC→GTC)  yjdF ← DUF2238 domain‑containing inner membrane protein YjdF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,345,4040GA98.7% 243.3 / NA 75A52V (GCC→GTC) yjdFDUF2238 domain‑containing inner membrane protein YjdF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (34/40);  minor base C (0/1);  total (34/41)
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

GAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGC  >  NC_000913/4345157‑4345649
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
gAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACTGATAACGTCTGACgg                                                                                                                                                                                                                                                      <  1:170904/251‑1 (MQ=255)
  aCGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTg                                                                                                                                                                                                                                                    <  1:50254/251‑1 (MQ=255)
               aCGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATATAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGc                                                                                                                                                                                                                                       <  1:92671/250‑1 (MQ=255)
                 gTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGc                                                                                                                                                                                                                                    <  1:349608/251‑1 (MQ=255)
                    cgcgCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGAAGCTGCACg                                                                                                                                                                                                                                 >  1:18875/1‑251 (MQ=255)
                           gTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCCGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACaata                                                                                                                                                                                                                           >  2:265861/1‑250 (MQ=255)
                                 cccGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCg                                                                                                                                                                                                                    >  2:358008/1‑251 (MQ=255)
                                   cgcgAACGAGAATTTCTCGTGCCACCAGTGCAGGCGCCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGgt                                                                                                                                                                                                                  <  2:162161/251‑1 (MQ=255)
                                    gcgAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGgtg                                                                                                                                                                                                                 <  2:564254/251‑1 (MQ=255)
                                      gAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCGTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTc                                                                                                                                                                                                               >  2:25350/1‑251 (MQ=255)
                                       aaCGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTc                                                                                                                                                                                                               >  1:321850/1‑250 (MQ=255)
                                         agagaATTTATCGTGCCACCAGTGCAGGCACCAGCAACTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACGTTCGCGTAGGTGTATTGTCCGGCGACCATCAGGATGATGGAGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACt                                                                                                                                                                                                            <  2:164496/250‑1 (MQ=255)
                                            gaATTTCTCGTGCACCCAGGGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTcc                                                                                                                                                                                                         <  2:202038/251‑1 (MQ=255)
                                               tttCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATc                                                                                                                                                                                                      >  1:164067/1‑251 (MQ=255)
                                                    cgggCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCTCTTCCATCAGCCa                                                                                                                                                                                                 <  2:312484/248‑1 (MQ=255)
                                                    cGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGcc                                                                                                                                                                                                  <  1:451332/250‑1 (MQ=255)
                                                    cGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCa                                                                                                                                                                                                 <  1:314943/251‑1 (MQ=255)
                                                     gTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACAATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGcc                                                                                                                                                                                                  >  2:64352/1‑249 (MQ=255)
                                                            caGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGtt                                                                                                                                                                                         <  1:99926/251‑1 (MQ=255)
                                                              gTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTg                                                                                                                                                                                       >  1:508392/1‑251 (MQ=255)
                                                                  aGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGtaat                                                                                                                                                                                   >  1:517892/1‑248 (MQ=255)
                                                                  aGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCAt                                                                                                                                                                                  <  2:367064/251‑1 (MQ=255)
                                                                      cccAGCTCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTTGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGc                                                                                                                                                                              <  2:88730/250‑1 (MQ=255)
                                                                        cAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCg                                                                                                                                                                            >  1:448598/1‑251 (MQ=255)
                                                                           ccccTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGaaa                                                                                                                                                                         <  1:403602/251‑1 (MQ=255)
                                                                              cTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGc                                                                                                                                                                      >  2:199150/1‑251 (MQ=255)
                                                                               tGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGcc                                                                                                                                                                     >  1:286663/1‑251 (MQ=255)
                                                                                   aaaaGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGt                                                                                                                                                                  >  1:554774/1‑250 (MQ=255)
                                                                                        tGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCCTTCCTGTACCTCACAACCAACGGGGACTTTCGCGTGGGTGAATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGGTAACGGATAACGTCTGACGGCGGCAAGCAGTGGCGGAACGACAATAATAACCGGTGTCACTTCCATCAGCCAGGGGGGTTGGTAATGGGCCGCATTGCCGTAATAAAt                                                                                                                                                            >  1:545208/1‑251 (MQ=255)
                                                                                        tGCCCCAGCTTGTAATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGAAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCAGGTATAAAt                                                                                                                                                            >  2:145888/1‑251 (MQ=255)
                                                                                         gCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGCAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATc                                                                                                                                                           >  2:488495/1‑251 (MQ=255)
                                                                                         gCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATc                                                                                                                                                           >  1:372571/1‑251 (MQ=255)
                                                                                            ccAGCTTGTCCTTCGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGCGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATCACGTCTGACGGTGGCCAGCAGTAGCTGCCCGTCAATAATCACCGGTGTCACTTCCATCAGCCAGGTGCGTTTGTCATGGGCAGACATGCCGGTACAAATCAgg                                                                                                                                                        >  2:324193/1‑251 (MQ=255)
                                                                                                    tCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGt                                                                                                                                                >  2:612696/1‑251 (MQ=255)
                                                                                                      ctACGGATTACGCCTCAACACCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGCGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAATATGATCTGATAGAGGAGCGGCGTTAACGGATAACGTGAGACGGTTGCAAGTAGTAGCTGCTCGACAATAATCACCTTTGTCACTTCCGTCAGCCAGGTGAGTGTGTCATGTGCCTAAATGCCGGTCTAAACAAGGATTAActgt                                                                                                                                               >  2:331422/2‑247 (MQ=255)
                                                                                                        aCGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGttt                                                                                                                                                                                        >  1:250031/1‑208 (MQ=255)
                                                                                                        aCGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGttt                                                                                                                                                                                        <  2:250031/208‑1 (MQ=255)
                                                                                                                cggcacAACCCCAGCAACTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCAACCATCAGGCTGATGGCGTGAAGGTAAATGAGCGTCTAGAGAAGCGGCGTTAACGGATGACGTCTGACGGCGGCAAGCAGCAGCTGCACGACAATAATCACCGGTGTCACCTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGACGGTATAAATTAGTATTAACGTTAGCGTCAGTgc                                                                                                                                    >  2:318856/6‑251 (MQ=255)
                                                                                                                cGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTgc                                                                                                                                    <  1:488495/251‑1 (MQ=255)
                                                                                                                      aaCCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGgt                                                                                                                              <  1:363017/251‑1 (MQ=255)
                                                                                                                       aCCCCAGCCATTCCTGCACCACAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAGATGAGCGTATAGAGCAGCTGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGGAGCTGGACGACAATAAGCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGAGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGgtg                                                                                                                             <  1:91413/251‑1 (MQ=255)
                                                                                                                        ccacAGCCATTCCGGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGTCGACCATAAGGATGATGGGGTGACGGAAAATGAGCGTAAAGAGCCGAGGCGTTAACGGATAACGTCTGACGGAGGCAAGCAGTAGCTGCACGACAATAATCACAGGAGTCACATCCATCAGCCAGGAGAGTTTGTCATGGGCAGAAATGCAGGTATAAATTAGGATTAACGTTAGAGTCAGTGCGCTGntgt                                                                                                                            >  2:604459/4‑251 (MQ=255)
                                                                                                                             gCCATTCCTGCACCTCAAAACCAACGGGGACGTTCGCGTAGGTGTATTGTCAGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACGTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTACAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGa                                                                                                                       <  1:324193/251‑1 (MQ=255)
                                                                                                                              ccATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGaa                                                                                                                      <  1:80606/251‑1 (MQ=255)
                                                                                                                                      gCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGc                                                                                                              <  1:366746/251‑1 (MQ=255)
                                                                                                                                       cACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCt                                                                                                             <  2:192509/251‑1 (MQ=255)
                                                                                                                                            cAAAACCAACGGGGACTTTCGAGTAGGTGTATTGTCCGCAGACACTCAGGATGCTGGCGTGAAGGAAAATGAGCGTATAGCGCAGCGGCGTTAACGGCTCACGCCTGACGGAGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTGCTTCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTCTAAATCAGGATTAACGTTAGCGTCAGTGCGATGGTGTTAAGAATTAACGGCTTgggt                                                                                                        >  2:316777/1‑248 (MQ=255)
                                                                                                                                                 ccAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCg                                                                                                   >  1:65857/1‑251 (MQ=255)
                                                                                                                                                     cGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTc                                                                                               <  1:638252/251‑1 (MQ=255)
                                                                                                                                                      ggggACTTTCGCGTAGGTGTATTGTGCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGTCAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCTGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGATTGAGTTTGCGGGTCa                                                                                              <  2:5761/251‑1 (MQ=255)
                                                                                                                                                          aCTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACCATAATCACCGGTGTCACTTCCATCAGc                                                                                                                                                                                                   >  1:418201/1‑147 (MQ=255)
                                                                                                                                                          aCTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGt                                                                                          <  1:109175/251‑1 (MQ=255)
                                                                                                                                                          aCTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGt                                                                                          <  1:390205/251‑1 (MQ=255)
                                                                                                                                                          aCTTGCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGc                                                                                                                                                                                                   <  2:418201/147‑1 (MQ=255)
                                                                                                                                                               cgcgTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGtt                                                                                     <  2:32776/251‑1 (MQ=255)
                                                                                                                                                                gcgTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTc                                                                                    <  1:202096/251‑1 (MQ=255)
                                                                                                                                                                    aGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCaa                                                                                <  2:164067/251‑1 (MQ=255)
                                                                                                                                                                     ggTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGAAGCTGCACGACAATACTCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTAAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGGGTGCGGGTCATGGTGAGTTAGCaa                                                                                >  1:616524/1‑250 (MQ=255)
                                                                                                                                                                     ggTGTAGTGTCCGACGACCACCAGGATGAGGGCGTGAAGGAAAATTAGCGTATATAGCGGCGGCGGTAACGGATAACGTCTGCCGGTGACAAGCAGTAGCTGCACGACTACAAACACCGTTGTCACTTCCATGAGCCAGTTGAGTGTGGCAGTGGGCGAAAGGTCTGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTTGTAAGAATTAACTGCTTTCGTGTTCGTGTCATGGttaggtggcaaa                                                                               <  2:408325/251‑12 (MQ=255)
                                                                                                                                                                        gtATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATc                                                                             >  1:451455/1‑250 (MQ=255)
                                                                                                                                                                           ttngccgccgACCTGCAGTATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAAAGGATGACGTCTGACTGTGTGGAAGTGTAGCTGCACGACAATAACCACGGGTGTCACTTCCATCGTCCAGGTGAGTTTGTTGTGGGCCGAAATGCCGGTATAAAGCACGATCTACGTCAGCGTCAGTGCGCTCGTGTTAAGTATTAACGGCTTGAGCGCGCGGGTTATGGTgtttttgcagatcagg                                                                          <  2:594181/246‑16 (MQ=255)
                                                                                                                                                                             gtcagccgACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTCTAGAGCAGCGGCGTTAACGTATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCCCCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCATGATTAACGTTAGCGTTAGTGCGCTGGTGTTAATAATTAACGGATTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGa                                                                        <  2:6631/246‑1 (MQ=255)
                                                                                                                                                                                 gccgACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATa                                                                   <  2:372571/251‑1 (MQ=255)
                                                                                                                                                                                 gccgACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATa                                                                   <  2:567511/251‑1 (MQ=255)
                                                                                                                                                                                   cgACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGACAATACt                                                                 >  1:559948/1‑251 (MQ=255)
                                                                                                                                                                                           aGGATGATGGCGGGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGCGAGTTTGTCATGGGCGGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGATTGTGCGGGTCAAGGGGAGGTGGCAAATCAGGAAAAAGACTGAcatttc                                                         >  2:53373/1‑246 (MQ=255)
                                                                                                                                                                                                    gCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAg                                                <  1:64352/251‑1 (MQ=255)
                                                                                                                                                                                                         aaGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGa                                           <  1:265861/251‑1 (MQ=255)
                                                                                                                                                                                                             aaaaTGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAAtt                                       >  1:570160/1‑251 (MQ=255)
                                                                                                                                                                                                                             gcagcGGCGTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAgg                       <  1:358008/251‑1 (MQ=255)
                                                                                                                                                                                                                                gcggcgTTAACGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTgg                    >  2:542122/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cGGATAACGTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGAt          <  1:542122/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTCTGACGGTGGCAAGCAGTAGCTGCAGGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGa                                           >  1:253700/1‑210 (MQ=255)
                                                                                                                                                                                                                                                  gTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGa                                           <  2:253700/210‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTCTGACGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGgc  >  1:462667/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
GAACGCCACCATTTTACGTCCGCGCACGTACATCCCGCGAACGAGAATTTCTCGTGCCACCAGTGCAGGCACCAGCCCCTGGAAAAAGTGCCCCAGCTTGTCATACGGATTACGACTCAACCCCAGCCATTCCTGCACCTCAAAACCAACGGGGACTTTCGCGTAGGTGTATTGTCCGCCGACCATCAGGATGATGGCGTGAAGGAAAATGAGCGTATAGAGCAGCGGCGTTAACGGATAACGTCTGGCGGTGGCAAGCAGTAGCTGCACGACAATAATCACCGGTGTCACTTCCATCAGCCAGGTGAGTTTGTCATGGGCCGAAATGCCGGTATAAATCAGGATTAACGTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGC  >  NC_000913/4345157‑4345649

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