Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,589,705 A→T 7.6% L101Q (CTG→CAG)  yjjA ← putative DUF2501 family periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,589,7050AT 7.6% 276.7 / 10.0 85L101Q (CTG→CAG) yjjAputative DUF2501 family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (35/43);  new base T (3/4);  total (38/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GGGAAAAGGGCAGGAGAGAGCCGCGTCAGCCGGGTAACAGCTAAAACGCGGCGTTTCGCGGTTCAGGAAATGAAGTTCAGCCCCTGTTTTAACACCAGATCGCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCA  >  CP009273/4589450‑4589974
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                             
gggAAAAGGGCAGGAGAGAGCCGCGTCAGCCGGGTAACAGCTAAAACGCGGCGTTTCGCGGTTCAGGAAATGAAGTTCAGCCCCTGTTTTAACACCAGATCGCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCGAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGt                                                                                                                                                                                                                                                             >  1:505561/1‑274 (MQ=255)
          cAGGAGAGAGCCGCGTCAGCCGGGTAACAGCTAAAACGCGGCGTTTCGCGGTTCAGGAAATGAAGTTCAGCCCCTGTTTTAACACCAGATCGCAGGCTTTGGTTTTCACCTTTTCTGCNAGCGGAGTTGTTACGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGc                                                                                                                                                                                                                                                      <  1:478023/271‑1 (MQ=255)
                                      agaaaaaaAGGGGCGTTTGGAGGGGTAGGAAATGACGTTAAGAACTTGTTTTAACACCAGATCGCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTcgcg                                                                                                                                                                                                                       <  2:323415/270‑1 (MQ=255)
                                      aGCTAATACGCGGCGTTTCGCGGTTCAGGAAATGAAGTTCAGCCCCTGTTTTAACACCAGATCGCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCACCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTCTTGCTTcgcg                                                                                                                                                                                                                       <  1:560195/274‑1 (MQ=255)
                                      aGCAAAAACGGGCNGTTTCGCGGTTAAGGCAATGAAGTTCCGCCCCTGTTTTAACTCCAGTTCGCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTcgcg                                                                                                                                                                                                                       <  2:420127/274‑1 (MQ=255)
                                        ctataatgtgtcgTTTCGCGGTTCAGGAAATGAAGTGCAGCCCCTGTTTTAACACCAGATCGCAGGTTTTTGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTTTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCa                                                                                                                                                                                                                     <  1:188066/263‑1 (MQ=255)
                                        cTATAACGCGGTGTTTTGCGGTTCAGGAAATGAAGTTCAGCACCTGTTTTAACACCAGATCGCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCa                                                                                                                                                                                                                     <  1:246245/274‑1 (MQ=255)
                                            aaCGCGGCGTTTCGCGGTTCAGGAAATGAAGTTCAGCCCCTGTTTTAACACCAGATCGCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATTTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCtt                                                                                                                                                                                                                           >  1:292823/1‑264 (MQ=255)
                                                                       gAAGTTCAGCCCCTGTTTTAACACCAGATCGCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCtt                                                                                                                                                                                                                           >  2:27120/1‑237 (MQ=255)
                                                                           ttCAGCCCCTGTTTTAACACCAGATCGCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTcgcg                                                                                                                                                                                                                       <  1:27120/237‑1 (MQ=255)
                                                                                  ccTCTTTTAACCCCAGTATGCATGTTTTGTTTTTTCCCTTTTCTGCCGGTGGAGTTTGTCCGATGTTATCCAGAGTGAGGTGCGGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGGCGAGATAGTTGGTGTCTGCTTTGTGCTCTTCACGGTTCAGGCCCAGCTTTTACAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTGGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCt                                                                                                                                                                           <  2:505554/274‑1 (MQ=255)
                                                                                                aGATCGCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGtttt                                                                                                                                                                                                                                                          >  1:456093/1‑181 (MQ=255)
                                                                                                aGATCGCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTc                                                                                                                                                                                                                                                         >  1:352947/1‑182 (MQ=255)
                                                                                                 gATCGCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTTCATAATGCCTGCGGCCTTGTTCCTGTTTTCTTCGCTTaattcctttttttct                                                                                                                                                           >  2:391865/1‑262 (MQ=255)
                                                                                                    cGCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTTGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGc                                                                                                                                                                                                                                                      <  2:456086/181‑1 (MQ=255)
                                                                                                    cGCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCa                                                                                                                                                                                                                                                     <  2:352942/182‑1 (MQ=255)
                                                                                                     gCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAANCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTCAGGCCTGGTTTCCgc                                                                                                                                                          <  1:391871/272‑1 (MQ=255)
                                                                                                      cAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCgctgct                                                                                                                                                      <  1:245916/275‑1 (MQ=255)
                                                                                                              ggTTTTCACCGTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCttt                                                                                                                                                                                                                                >  1:563921/1‑193 (MQ=255)
                                                                                                                  ttCACCTTTTCTGCCATCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCt                                                                                                                                                                                                                            <  2:563913/193‑1 (MQ=255)
                                                                                                                                                     aTCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGa                                                                                                                                                                                                                                          >  2:166438/1‑144 (MQ=255)
                                                                                                                                                     aTCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGc                                                                                                                                                                                                >  2:33701/1‑186 (MQ=255)
                                                                                                                                                     aTCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAg                                                                                                                                                   >  2:555244/1‑231 (MQ=255)
                                                                                                                                                         aGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGcc                                                                                                                                                                                                                                      <  1:166439/144‑1 (MQ=255)
                                                                                                                                                         aGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTg                                                                                                                                                                                            <  1:33701/186‑1 (MQ=255)
                                                                                                                                                         aGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAag                                                                                                                                               <  1:555252/231‑1 (MQ=255)
                                                                                                                                                           aTTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTAc                                                                                                                                                                                                                                             >  2:562899/1‑135 (MQ=255)
                                                                                                                                                           aTTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGc                                                                                                                                                                                                                             >  1:112501/1‑151 (MQ=255)
                                                                                                                                                           aTTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCtt                                                                                                                                                                                                                           >  1:209858/1‑153 (MQ=255)
                                                                                                                                                               aGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGAc                                                                                                                                                                                                                                         <  1:562907/135‑1 (MQ=255)
                                                                                                                                                               aGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTcgcg                                                                                                                                                                                                                       <  2:209856/153‑1 (MQ=255)
                                                                                                                                                               aGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTcg                                                                                                                                                                                                                         <  2:112501/151‑1 (MQ=255)
                                                                                                                                                                gTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCAt                                                                                                                                                                                              >  1:313033/1‑177 (MQ=255)
                                                                                                                                                                   gCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTc                                                                                                               <  1:335092/253‑1 (MQ=255)
                                                                                                                                                                    cTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGc                                                                                                                                                                                                                             >  2:413957/1‑142 (MQ=255)
                                                                                                                                                                    cTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGtt                                                                                                                                                                                          <  2:313029/177‑1 (MQ=255)
                                                                                                                                                                    cTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTCGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCtt                                                                                                                                                                                                                           >  1:543558/1‑144 (MQ=255)
                                                                                                                                                                     tGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCtt                                                                                                                                                                                                                           >  1:349425/1‑143 (MQ=255)
                                                                                                                                                                        ccATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTcgcg                                                                                                                                                                                                                       <  2:543550/144‑1 (MQ=255)
                                                                                                                                                                        ccATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTcg                                                                                                                                                                                                                         <  1:413964/142‑1 (MQ=255)
                                                                                                                                                                         cATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTcgcg                                                                                                                                                                                                                       <  2:349420/143‑1 (MQ=255)
                                                                                                                                                                                               cTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCa                                                                                                                                                                                                                                                                               <  1:121050/65‑1 (MQ=255)
                                                                                                                                                                                                       cGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCACGTTCGCTGGCGGCGCTGGCGAGCGCATCttt                                                      >  1:178349/1‑274 (MQ=255)
                                                                                                                                                                                                        gTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGc                                                                                                                                                                                                                             >  1:511065/1‑106 (MQ=255)
                                                                                                                                                                                                            aGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTcg                                                                                                                                                                                                                         <  2:511058/106‑1 (MQ=255)
                                                                                                                                                                                                              aTAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCtt                                                                                                                                                                                                                           >  1:86122/1‑102 (MQ=255)
                                                                                                                                                                                                                  ttGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTcgct                                                                                                                                                                                                                       <  2:86122/102‑2 (MQ=255)
                                                                                                                                                                                                                    ggTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAg                                                                                                                                                   >  1:18178/1‑168 (MQ=255)
                                                                                                                                                                                                                        tcttctTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAag                                                                                                                                               <  2:18178/168‑1 (MQ=255)
                                                                                                                                                                                                                               tttGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCtt                                                                                                                                                                                                                           >  2:523604/1‑85 (MQ=255)
                                                                                                                                                                                                                                  gCTCTTCACTGTTCAGGCCCAGCTTTTCCTGCTCCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTGTGCTTCGCGCCGTANTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGGCGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTAGCCACCACGCAGCCTGCACCCCGGTGGaa                          <  1:394633/275‑1 (MQ=255)
                                                                                                                                                                                                                                   ctctTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTcgcg                                                                                                                                                                                                                       <  1:523612/85‑1 (MQ=255)
                                                                                                                                                                                                                                             ttCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCg                                                                                                                                                                                                                                                        >  2:125953/1‑42 (MQ=255)
                                                                                                                                                                                                                                             ttCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCgg                                                                                                                                                                                                                                                 >  2:278328/1‑49 (MQ=255)
                                                                                                                                                                                                                                             ttCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCa                                                                                                                                                                                                                     >  1:167045/1‑77 (MQ=255)
                                                                                                                                                                                                                                                aGGCCCATCTTTTCCAGCACCTGGTTCTTGATGTTTTCTGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTggc            <  1:141425/275‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCAt                                                                                                                                                                                                                                                    <  1:125954/42‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTAc                                                                                                                                                                                                                                             <  1:278330/49‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTAt                                                                                                                                                                                                                 <  2:167044/77‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCtt                                                                                                                                                                                                                           >  1:490220/1‑66 (MQ=255)
                                                                                                                                                                                                                                                  gCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGc                                                                                                                                                                                                >  1:75007/1‑93 (MQ=255)
                                                                                                                                                                                                                                                  gCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGGGGCGCTGGGGAGCGCAGCTTTCCCGGACGCAGCAGGCCACCCGGTGGGACNTAATGCGCTgctc           >  2:422625/1‑271 (MQ=255)
                                                                                                                                                                                                                                                  gCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTAACGACGCCAGCTTTTGCtt                                                                                                                                                                                                                           >  2:399664/1‑66 (MQ=255)
                                                                                                                                                                                                                                                     cccttttcCCTGCTCCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCCGTATTGCAGAATGCCTTCGGc                                                                                                                                                                                                >  2:320506/5‑90 (MQ=255)
                                                                                                                                                                                                                                                     cccttttcCCTGCTCCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGc                                                                                                                                                                                                >  2:215807/5‑90 (MQ=255)
                                                                                                                                                                                                                                                     cccttttcCCTGCTCCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGc                                                                                                                                                                                                >  1:16242/5‑90 (MQ=255)
                                                                                                                                                                                                                                                     cAGCTTTTCCAGCACCTGGTTCTTGATTTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCtt                                                                                                                                                                                                                           >  1:333771/1‑63 (MQ=255)
                                                                                                                                                                                                                                                      aGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTcgcg                                                                                                                                                                                                                       <  2:490213/66‑1 (MQ=255)
                                                                                                                                                                                                                                                      aGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTg                                                                                                                                                                                            <  2:75007/93‑1 (MQ=255)
                                                                                                                                                                                                                                                      aGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACAGACGCAAGCTTTTGCTTcgcg                                                                                                                                                                                                                       <  1:399670/66‑1 (MQ=255)
                                                                                                                                                                                                                                                       gCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCtt                                                                                                                                                                                                                           >  2:143004/1‑61 (MQ=255)
                                                                                                                                                                                                                                                       gCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGc                                                                                                                                                                                                   >  2:31088/1‑85 (MQ=255)
                                                                                                                                                                                                                                                       gCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGGCCAACCGCCTTCCTGTTTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGGCGCATGCACCCCGTGGGAAATTATTGAGCTCGCGGCAATg     >  1:246234/1‑275 (MQ=255)
                                                                                                                                                                                                                                                       gCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACCCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCCCCCCGGTGGAAATTAATGCGCTGGCGGCACTg     >  1:450994/1‑275 (MQ=255)
                                                                                                                                                                                                                                                       gCTTTTCCAGCACCTGGTTCTTGATGTATTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGttt                                                                                                                                                              >  1:336037/1‑122 (MQ=255)
                                                                                                                                                                                                                                                         ttttCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTcgcg                                                                                                                                                                                                                       <  2:333766/63‑1 (MQ=255)
                                                                                                                                                                                                                                                         tttcCCTGCTCCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTTCAGAATGCCTGCGGCATTg                                                                                                                                                                                            <  1:320511/90‑1 (MQ=255)
                                                                                                                                                                                                                                                         tttcCCTGCTCCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTg                                                                                                                                                                                            <  1:215809/90‑1 (MQ=255)
                                                                                                                                                                                                                                                         tttcCCTGCTCCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTg                                                                                                                                                                                            <  2:16242/90‑1 (MQ=255)
                                                                                                                                                                                                                                                          tttCCAGCTCCTTGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTCTCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGgca  <  1:302165/275‑1 (MQ=255)
                                                                                                                                                                                                                                                          tttCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGNTGCTAAGCAAGTAAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGgca  <  1:475638/275‑1 (MQ=255)
                                                                                                                                                                                                                                                           ttCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTcgcg                                                                                                                                                                                                                       <  1:143005/61‑1 (MQ=255)
                                                                                                                                                                                                                                                           ttCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCa                                                                                                                                                                                               <  1:31088/85‑1 (MQ=255)
                                                                                                                                                                                                                                                           ttCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCgc                                                                                                                                                          <  2:336032/122‑1 (MQ=255)
                                                                                                                                                                                                                                                           ttCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTGACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTTGTTGCCAAGTTCGCTggcg                                                                          <  1:465698/202‑1 (MQ=255)
                                                                                                                                                                                                                                                           ctCCAGGACCTGGTTCTGGATGTTTTCGGCATCGGTTACCGATGTCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTGATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGGTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGgca  <  2:178348/273‑1 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                             
GGGAAAAGGGCAGGAGAGAGCCGCGTCAGCCGGGTAACAGCTAAAACGCGGCGTTTCGCGGTTCAGGAAATGAAGTTCAGCCCCTGTTTTAACACCAGATCGCAGGCTTTGGTTTTCACCTTTTCTGCCAGCGGAGTTGTTCCGATGTTATCCAGATTGAGTTGCTGACCATCTTTGGTTTTCAGCAAACCCTGAATGCCGTCCAGATAGTTGGTGTCTTCTTTTTGCTCTTCACTGTTCAGGCCCAGCTTTTCCAGCACCTGGTTCTTGATGTTTTCGGCATCGGTTACCGACGCCAGCTTTTGCTTCGCGCAGTATTGCAGAATGCCTGCGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCA  >  CP009273/4589450‑4589974

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