Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 514,434 C→T G142G (GGG→GGA ybbJ ← NfeD‑like family protein YbbJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913514,4340CT97.3% 222.2 / ‑4.5 73G142G (GGG→GGAybbJNfeD‑like family protein YbbJ
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (35/36);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTC  >  NC_000913/514188‑514671
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
cTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGGGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCtt                                                                                                                                                                                                                                            >  1:1592/1‑251 (MQ=255)
  ggTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTCTTCCTTCt                                                                                                                                                                                                                                          >  1:27068/1‑251 (MQ=255)
           gTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATg                                                                                                                                                                                                                                 <  1:456025/251‑1 (MQ=255)
             ttAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGa                                                                                                                                                                                                                                <  2:343711/250‑1 (MQ=255)
               aaCGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACtt                                                                                                                                                                                                                             >  1:399514/1‑251 (MQ=255)
               aaCGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACtt                                                                                                                                                                                                                             <  1:18619/251‑1 (MQ=255)
                  gACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCaa                                                                                                                                                                                                                          <  2:76325/251‑1 (MQ=255)
                      cGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGACGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTCTCGAAAATCTCTCCGGCTGCTGTCATAATCGGGCAGGGAGATTAAGACGAGACGGCCCGGCTATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATg                                                                                                                                                                                                                      >  1:199517/1‑251 (MQ=255)
                         aGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGt                                                                                                                                                                                                                   >  1:373369/1‑251 (MQ=255)
                           cGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTAc                                                                                                                                                                                                                 <  1:13783/251‑1 (MQ=255)
                           cGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTAc                                                                                                                                                                                                                 >  1:448452/1‑251 (MQ=255)
                               acaGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGGGACGGCTCGGGTATGTAGGGTTATTCCTTCTATCGCAAAGACTTCAACATGCGTACCTGc                                                                                                                                                                                                             >  2:354605/1‑251 (MQ=255)
                                caGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCa                                                                                                                                                                                                            <  1:186726/251‑1 (MQ=255)
                                  gcgcCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCAcc                                                                                                                                                                                                          <  1:284847/251‑1 (MQ=255)
                                       gACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGACAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCg                                                                                                                                                                                                         >  1:121789/1‑247 (MQ=255)
                                                      aCGCTGGGGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGa                                                                                                                                                                                      <  2:171699/251‑1 (MQ=255)
                                                      aCGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGa                                                                                                                                                                                      <  1:382678/251‑1 (MQ=255)
                                                        gcgggAGAAGGGGGCGGAGATCGCACGACACATTGAGGAGCTGCACACCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACa                                                                                                                                                                                    <  2:191744/248‑1 (MQ=255)
                                                         cTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGAc                                                                                                                                                                                     >  2:94831/1‑249 (MQ=255)
                                                           ggAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGAAGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGc                                                                                                                                                                                 >  1:32216/1‑251 (MQ=255)
                                                                aGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACGGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCGCGGATATGTAGGGTTATTCCTTCGATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGa                                                                                                                                                                            <  1:63779/251‑1 (MQ=255)
                                                                    ggCGGAGATCGAACGACCCATGTAGGAGCTGCATTCCATGCGCGACCAGCTGCTGGCCCTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGACCCTCGCTGGCGCTGACAGGCCATGAACt                                                                                                                                                                         <  2:187877/250‑1 (MQ=255)
                                                                    ggCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGGGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTg                                                                                                                                                                        >  1:120736/1‑251 (MQ=255)
                                                                             cGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGAc                                                                                                                                                                <  1:91498/250‑1 (MQ=255)
                                                                             cGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACCGTCACCGAcc                                                                                                                                                               >  2:324402/1‑250 (MQ=255)
                                                                                acggtacaTTGAGGAGCTGCATGCCATGCGAGACCCGCTGCTTGCTCTGGCGAATGCCTGCCCTGGCGATGCCAGCGCCGACTGCCCGATTATCGAAAATCTCGCCGGCTGTTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCGTCTATCGCAATGACTTCAACATGCGTCCCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACgcgc                                                                                                                                                            <  2:52294/246‑1 (MQ=255)
                                                                                 cGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGCCTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGGTATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCCCTGACCGGCCATGAACTGTCACCGAcgctcc                                                                                                                                                           >  2:293120/1‑248 (MQ=255)
                                                                                    cacaTTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCatat                                                                                                                                                        >  1:144354/1‑251 (MQ=255)
                                                                                     acaTTGAGGAGCTGCAAGCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATg                                                                                                                                                       <  1:471110/251‑1 (MQ=255)
                                                                                       atggaggagCTGCAATCCATGCGCGACCAGCTGCTGGCACGGGCGAAAGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCAGCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGTCAGGCCATGAACTGTCACCGACGCGCATATGAc                                                                                                                                                     <  2:234604/248‑1 (MQ=255)
                                                                                             gagtgtcAATCCAGGCGCGACCAGCTGCGGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCa                                                                                                                                                           <  2:178639/233‑1 (MQ=255)
                                                                                             gagCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCa                                                                                                                                                           >  1:178639/1‑239 (MQ=255)
                                                                                                cTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGt                                                                                                                                            >  1:128455/1‑251 (MQ=255)
                                                                                                  gCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGtt                                                                                                                                           <  2:144354/250‑1 (MQ=255)
                                                                                                    aaTCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGACCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACATGCACCGCGATCCTCGCTGGCGCTGACCGGCCATGAACTGTCACCGACGCGCCTATGACCGCGCCCGTTGAc                                                                                                                                        >  2:225501/1‑251 (MQ=255)
                                                                                                                  cAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATg                                                                                                                                                                                                                      <  1:113825/159‑1 (MQ=255)
                                                                                                                  cAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATg                                                                                                                                                                                                                      >  2:113825/1‑159 (MQ=255)
                                                                                                                        ctgGCACTGGCGAATGCCTGCACTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCCATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCCCCGCCAGGCCCTGAACCGTCCCCCACGCGCACATGCCCGCGCCCgcttaccaccgcacattccccaccc                                                                                                                     >  2:225402/1‑226 (MQ=255)
                                                                                                                          ggCACTGTCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACaa                                                                                                                    <  1:290995/249‑1 (MQ=255)
                                                                                                                                ggCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCgg                                                                                                            <  1:397540/251‑1 (MQ=255)
                                                                                                                                  cGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGcc                                                                                                          >  1:131056/1‑251 (MQ=255)
                                                                                                                                                   gATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCgag                                                                                                                                                                                                       >  2:365496/1‑141 (MQ=255)
                                                                                                                                                   gATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCgag                                                                                                                                                                                                       <  1:365496/141‑1 (MQ=255)
                                                                                                                                                    aTGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGTCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCAcgcg                                                                                        <  1:483327/251‑1 (MQ=255)
                                                                                                                                                        aaGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTAGCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTgg                                                                                    <  2:280846/250‑1 (MQ=255)
                                                                                                                                                         aGCGGCGACTGCCCGATTATCGAAAATCTCTCCGGATGTTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTAGCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCGCGCGGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGATGCCCGCGCTGGt                                                                                   <  2:399514/251‑1 (MQ=255)
                                                                                                                                                              cGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTaaa                                                                              <  1:464759/251‑1 (MQ=255)
                                                                                                                                                               gACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCCCAAAACGTCGGCCAATCAGCTGCTGCCCGCGCAGGTTTAAAt                                                                             >  2:41729/1‑251 (MQ=255)
                                                                                                                                                                aCTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCGGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATg                                                                            <  1:146692/251‑1 (MQ=255)
                                                                                                                                                                aCTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATg                                                                            <  1:238050/251‑1 (MQ=255)
                                                                                                                                                                     ccGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTAt                                                                       >  1:232223/1‑251 (MQ=255)
                                                                                                                                                                      cGATTAACGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATc                                                                      <  2:169791/251‑1 (MQ=255)
                                                                                                                                                                        aTTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCCATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGc                                                                    >  2:425627/1‑251 (MQ=255)
                                                                                                                                                                                    ctctcCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGtt                                                        >  2:24528/1‑251 (MQ=255)
                                                                                                                                                                                      ctcCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATTTGGAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACAGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGATTAAATGACTATCGCTGTGCTTTTGTTcg                                                      >  1:123406/1‑251 (MQ=255)
                                                                                                                                                                                      ctcCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGt                                                         <  1:293120/248‑1 (MQ=255)
                                                                                                                                                                                      ctcCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTcg                                                      >  1:123387/1‑251 (MQ=255)
                                                                                                                                                                                        ccGGCTGCTGTCCTCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATCTGACCGCGCCCGTTGACCAGCGGAGCTTCCAGCCCAACACGTCGGCCACTCAGCTGCTGCCCGCGCTGGTTTAACTGACtctagatcggaaaagcctcgtg                                                    >  2:340406/1‑230 (MQ=255)
                                                                                                                                                                                        ccGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTATGGTTATTCCTTCTAGCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTcgc                                                     <  2:12423/250‑1 (MQ=255)
                                                                                                                                                                                        ccGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGAAAAAATTCTAGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTAt                                                                       <  1:340406/232‑1 (MQ=255)
                                                                                                                                                                                        ccGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAAGCGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATTACCTCGCCCGTTGACCAGCGGATATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTc                                                       <  1:226287/248‑1 (MQ=255)
                                                                                                                                                                                        ccGGCTGCGGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTAGCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGGGCTTTTGTTcg                                                      <  2:448452/249‑1 (MQ=255)
                                                                                                                                                                                           gctgctGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCAc                                                 <  1:304257/251‑1 (MQ=255)
                                                                                                                                                                                                      catcGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGAc                                       >  1:371871/1‑250 (MQ=255)
                                                                                                                                                                                                         cGGGCAGGGTGATTAAGACGAGACGGCTCGGCTATGTAGGGTTATTCATTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAAcc                                   >  1:409959/1‑251 (MQ=255)
                                                                                                                                                                                                            gCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCCATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGAc                                       >  1:336093/1‑244 (MQ=255)
                                                                                                                                                                                                                ggTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTccc                           >  1:191399/1‑250 (MQ=255)
                                                                                                                                                                                                                    aTTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACccgcc                                            <  1:95685/231‑1 (MQ=255)
                                                                                                                                                                                                                    aTTAAGACGAGACGGCTCGGATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGGGCTTTTGTTCGCGCACccgcc                                            >  2:95685/1‑231 (MQ=255)
                                                                                                                                                                                                                                      ggATATGTAGGGTTATTCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGc                    >  2:403758/1‑236 (MQ=255)
                                                                                                                                                                                                                                      ggATATGGAGGTTTATCCCTTCTATCGCAATGACTCCAACAGGCGTCCCGGCCCCGAGACCCTCGCGGGCGCTGACAGGCCAGGAACTGTCACCGACGCGAATAGGACCGCGCCCGTTGCCCAGCGGAGATTCCAGCAAAAAACGGCGGCAAATCAGCGGCTGCCCGCGTTGGTTTAAATGACTATCGCGGTGTTTTTGTCCGCGCACCCGCCGCGACAACCATTTCCACCACAGc                    <  1:403758/236‑1 (MQ=255)
                                                                                                                                                                                                                                        atatGTCGGGTTATTCCTTATATCGCCATGACTTCTACATGCGTACCTGCGCCGAGATCCTCGCTGGCCCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTCGACCAGCGGAGCTTCCCGAACCAAACGTCGGCCCATCAGCTGCTGCCCGCGCTGGTTTCATTGACTACCGCAGTCCTTTTGTTCGTGCCCCCGCCGCGCAACCCATTTacacccaatcacggcggcaagcagcgt   >  2:51592/1‑229 (MQ=255)
                                                                                                                                                                                                                                          atgtgggtttATTCCTTCTACCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCCCGCTGGCGCTGACAGCCCATGACCTGCCTCCGACGCGCATATGCCCGCGCCCTTTGACCAGCGGAGATTCCAGCAAAAACCGTCGGCCAATCAGCTGCTCCCCGCGTTGGTTTAAATGCCTATCGCTGTGTTTTTGTGCGCTCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTc  <  1:393374/242‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
CTGGTGAATCTGTTTAACGACCCGCAGCGGCACAGCGCCGACGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTC  >  NC_000913/514188‑514671

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