Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,151,329 T→G 31.2% A590A (GCA→GCC ppc ← phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,151,3290TG31.2% 78.3 / 3.4 62A590A (GCA→GCCppcphosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base T (29/7);  new base G (7/14);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 5.83e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTT  >  NC_000913/4151094‑4151565
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
ggTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTgcg                                                                                                                                                                                                                                              <  1:1188548/239‑1 (MQ=255)
         cGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGCGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGGTGTGACCACAGCCCCGCATGAGCCGGGGCGCCGCcgcg                                                                                                                                                                                                                                    >  1:93919/1‑238 (MQ=255)
            aGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCTTTTCAGGCTTCCTGGCGGGTGTGACAGCAGCGCCGCATGAGCAGGGGCGCCGCCGCGAccc                                                                                                                                                                                                                                >  2:1079390/1‑238 (MQ=255)
            aGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCCCGACCa                                                                                                                                                                                                                                >  1:103604/1‑239 (MQ=255)
              gagaTGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTTTGACAGCAGCGCCGCATGAGCAGGGGCGCCGCCGCGACCAAt                                                                                                                                                                                                                             >  1:1157237/1‑239 (MQ=255)
                   gactgacAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGGGCGCCGCCGCGACCAATGGaa                                                                                                                                                                                                                         >  1:385862/1‑238 (MQ=255)
                         taGTTCCGCCATAGGGCGACGCCGGCTTTTTTTCGGCGGCGGGGGGGGCGGCAGGTTGGCTCCCAGACCCGCCCCGGTAAAAGGCGACAGGCTGCGGGGGGTTTTTTCTGGCAGACCATATTTAAAGCGGATCATCCCGCCCTTTTCGTTTACGCGCGGCCCGCTTTTCAGGCTTCCGGGCGTTTGTGACAGCAGCGCCGCAGGAGCGGGGGCGCCGCCGCGCCAAAGGGAACcgccgc                                                                                                                                                                                                                  <  2:914727/238‑1 (MQ=255)
                                       tGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGCCCGTGGaa                                                                                                                                                                                                        >  1:583073/1‑235 (MQ=255)
                                           aCGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACCGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTTGTACAGCAGCCCCGCCTCAGCAGGCGCCGCGCCCGCACCAATGGGACCGCCCCGCCAGCTGGACAACGtgc                                                                                                                                                                                                >  2:870513/1‑237 (MQ=255)
                                           aCGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCAAGAGCCGGTGCGCCGCCGCGGCCCATGGGACCGCCGCGCCCGGGGGACAACg                                                                                                                                                                                                   >  1:521087/1‑236 (MQ=255)
                                                  ctctctTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCCACCAATGGGACCGCCCCGCCCGTGGAACAACGTTAGCTcgac                                                                                                                                                                                         >  1:564336/1‑236 (MQ=255)
                                                  ctctctTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGGGGCAGCAACCCCGCATGAGCAGGGGCGCCGCCGCGACCAAAGGCACCGCCGCGACCGTGGAGCAACGGCAGCTCAAt                                                                                                                                                                                         >  2:744715/1‑239 (MQ=255)
                                                                 ggcggTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGGGGAACAACGTCACGTCAATACCCGCCTTTTCGCa                                                                                                                                                                          >  1:722216/1‑239 (MQ=255)
                                                                     ggggcggcGGGTGGGCTTCCAGACCCGCCCGGGTAGAAAGCGACGGGTTGCTGGCGTTGTTTTCTGGCAGCCTATTTTTAAGGGGTATCCTCCCGCCTTGTTGGGTTCCGCGCGGGCCGCCTTTCGGGCTTCCTGGGGGTGGGGCCGGCGGGCCCGCAGGAGCGGGCGCGCCGCCGCGACAATTGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGt                                                                                                                                                                       <  2:103604/236‑1 (MQ=255)
                                                                            agggTTGGCTCACAGAATCGCCCGGTTAAAAAGCGACGGGCTGCTGACGGTGATTTCGGGCAGACAATATTTAAGGGGGATCTTCTCGCCCTGTTGGGTTCCGCGCGGGCCCCTTTCCAGCTTTCGTGGGGGTTGTGACAGCGGCGCGGCATGGGCGGGGGCGCCGCCGCGCCAAGTGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGAtt                                                                                                                                                                <  2:1237554/236‑1 (MQ=255)
                                                                                gtGGCCTCCAGACTCCCCCGGGTATAAGGCGACAGGTTGCGGGGGGTGTTTCCTGCCAGACTATATTTAAAGCGACTCCTCCCCCTTTGTGCGGTCCCGCGCGGGCCGCTTTTCGGGCTTCCGGGGGGTTGTGCCAGCAGCGCCGCAGGAGCGGGGGCGCCGCCGCGACCAAGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATg                                                                                                                                                            <  2:1157237/237‑1 (MQ=255)
                                                                                         aGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCtgtg                                                                                                                                                   >  1:366494/1‑238 (MQ=255)
                                                                                             tCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGa                                                                                                                                              >  2:950754/1‑239 (MQ=255)
                                                                                                               aCAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATg                                                                                                                                                            >  1:720331/1‑207 (MQ=255)
                                                                                                               aCAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGaa                                                                                                                               >  1:25981/1‑236 (MQ=255)
                                                                                                               aCAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAggt                                                                                                                            >  1:754438/1‑237 (MQ=255)
                                                                                                                cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCCATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCt                                                                                                                            >  1:147317/1‑238 (MQ=255)
                                                                                                                cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGAGCAGCCGCGACCAATGGAACCGCCGCGACCGGGGAACAACGTAAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATACATCCTGTGCCTGAGTTTGGCCCCAGGAAGCt                                                                                                                            >  2:1085197/1‑238 (MQ=255)
                                                                                                                cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGCTGTGACAGCCGCGCCGCATGAGCAGGTGCGCCGCCCCGACCAATGGAACCGCCGCGACCGTGGAAAAACGTCAGCTCACTCCCCGCTTTTTCGCAGGGTTTGGTTAATGCATCCTGTGCCTGACATTCCGCCCAGGAAGCt                                                                                                                            >  2:7899/1‑238 (MQ=255)
                                                                                                                cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTGAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACCGCAGCGCCGCATGAGCAGGGGAGCCGCAGAGACCAATGGTAACGCCGCGACCGTGGGACAACGTCAGCACAATACCCCATTTTTCTCTGGATTTTATTTATGCGTCCTGGGCCTGATATAGTGCCCAGAAAGCt                                                                                                                            >  2:248115/1‑238 (MQ=255)
                                                                                                                                           aCCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCCGCCCAATACCCGCCTTTTCGCCGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGGAGCTGCCATCACTCCCTCATCTTTTGCttgg                                                                                                 >  2:1150789/1‑235 (MQ=255)
                                                                                                                                                      aGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAACGCCGCCTGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGGGGAACAACGTCAGCTCAATCCCCGCTTTTTCGCAGGCTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCACGAAGCTGCCATCACTCCCGCATCTTTTTCTGAGTCGGAATAGcgc                                                                                     >  2:1159379/1‑237 (MQ=255)
                                                                                                                                                          gATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGCTATTGCGCCCAGGAAGCTGCCATCACTCCCGCCTCTTTTGCTGAGTCGGAATAGCCAat                                                                                   >  1:271007/1‑237 (MQ=255)
                                                                                                                                                          gATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAatca                                                                                 >  1:316002/1‑239 (MQ=255)
                                                                                                                                                           atcatcTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGc                                                                                                                                                                                              >  2:860045/1‑129 (MQ=255)
                                                                                                                                                           atcatcTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGCGGCACAACGTCAGGTCAATACCCGCTTTTTCGCAGGGTTTTATTAATGTATCCTGGGCCTGATATTGCGCCCAGGAACCTGCTAACACTCCCGCATCTTTTCCTGAGGCGGAATAGCCAataa                                                                                 >  2:997699/1‑236 (MQ=255)
                                                                                                                                                            cccccccccccTGTCGGTTACCGCGCGGCCCCCCTTCCGGCCTCCCGGGCGGTGGTGACAGCAGCGCCGCGTGAGCGGGCGCCCCCCCGCGCCAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAatcat                                                                                <  2:25981/230‑1 (MQ=255)
                                                                                                                                                                tccGCCCTGTCGGTGACGGGCCGGCCGGCCGTCCGGGCTTCTGGGGGGTTGTGACGCAGGGCCAGAGGGGGAAGGGGCGCGGCCGCGCCCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCa                                                                           <  2:719938/237‑1 (MQ=255)
                                                                                                                                                                 cgcgccTGTTCCGGTGCCGCGCGGGCGGCCTTGCAGCCTTCCGGGGGGTTGTGACAGCAGCGCCGCAGGGGCAGGTGCGCCGCCGCGACCAATGGAGCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCa                                                                           <  1:121550/235‑1 (MQ=255)
                                                                                                                                                                     ccTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGGCGCATGAGCGGGGGCGCCGCCGCGGCCAAAGGAACCGGCGTGACCGTGGAACAACGTGAGCTCAATAACCGCTTTTTCGCAGGTTTTGGTTAAGGCCTGCTGTGCCTGATATTGCGCCCAGGAAGCTGCCCTCACTCCCGCATCTTTTGCTGATTCGGAATAGCCCATCATCAACATCTGt                                                                      >  2:89008/1‑239 (MQ=255)
                                                                                                                                                                         ttCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCCTCTTTTGCTGGGTCGGAATAGCCAATCATCACCATCTGTTTg                                                                   >  1:484350/1‑238 (MQ=255)
                                                                                                                                                                            gcgttaccgcgcgGGCCGCTTTACGGGCTCCCGGGCGGTTGTGACAGCAGCGCCACATGAGCGGGTGCGCCGCCGCGACAAAGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGcc                                                                 <  1:930179/232‑1 (MQ=255)
                                                                                                                                                                               ttACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGa                                                             >  1:843362/1‑239 (MQ=255)
                                                                                                                                                                                    ggccgggCCGCCTTCCAGGCTCCGGGGCGGTGTTGACACCCGCCGCGCGTGAGCGGGGGCCCCCCCGCACCCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCa                                                         <  2:883544/234‑1 (MQ=255)
                                                                                                                                                                                     cgcAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCCATGGAACCGCCGCGACCGTGGAACAACGTCCGCTCAATACCCGCTTTTTCGCCGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCCCTCCCGCCTCTTTTGTTGAGTCGGAATAACCAATCATCCCCATCTGTTTGCCCTGAATccg                                                        >  2:706456/1‑237 (MQ=255)
                                                                                                                                                                                     cgcAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTTCTGAGTCGGAATAGCCAATCATCCCCCTCTGTTTGCCCTGAATCAg                                                        >  2:282279/1‑239 (MQ=255)
                                                                                                                                                                                        aGGCCGCCTTTCGGGCTTCCTGGCGGTTGTGACGGCAGCGCCGCAGGAGCGGGGGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGc                                                      <  1:1068068/238‑1 (MQ=255)
                                                                                                                                                                                          gccgccTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCAc                                                   >  2:741318/1‑239 (MQ=255)
                                                                                                                                                                                          gccgccTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCAc                                                   >  1:883544/1‑239 (MQ=255)
                                                                                                                                                                                           gtccccTTTCCTGCGTCCCGGTGGTGTGTGACGGGAGCGCCGCATGGGCGGCTCCGCCCCCGCGCACATAGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCa                                                    <  2:564336/235‑1 (MQ=255)
                                                                                                                                                                                               tcttcagggCTTCCTGGGGGTGGGCACAGCGGCGCCGGAGGAGCGGGGGCGCCGCCGCGACAAAGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa                                               <  2:271007/232‑1 (MQ=255)
                                                                                                                                                                                               ccTTTCAGGCTTTCGGGCGGTTGTGACGACAGCGGCGTGGGGGCAGGTGCGCCCCCGCAACCAGTGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa                                               <  1:283325/239‑1 (MQ=255)
                                                                                                                                                                                               ccTTTCAGGCTTCCTGGCGGTTGTGACAGCCGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTTCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGAt                                                >  1:969262/1‑238 (MQ=255)
                                                                                                                                                                                                cgtccAGGCTTCCGGGGGGTTGCAGCAGCCGCCCCGCAGGAGCAGGGGCCCCCCCGCGCCAAATGGACCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa                                               <  2:366494/234‑1 (MQ=255)
                                                                                                                                                                                                    cAGCCTTCCTGGTGGTTGTGCCAGGAGCGCCGCAGGAGCAGGGCCCCCCCCGCGCCCAAGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGt                                          <  2:316002/239‑1 (MQ=255)
                                                                                                                                                                                                                    tgtggACAGCAGCCCCGCATGAGCAGGGCCCCCGCCGCAACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGagcagc                            <  1:747851/234‑1 (MQ=255)
                                                                                                                                                                                                                              ggcgcCACAGGAGCGGGGGCGCCGCGGCGACAAAGGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGAc                 <  2:266304/238‑1 (MQ=255)
                                                                                                                                                                                                                              cgcgcGGCAGGAGCAGGGCGCCCCACGGGGCCATATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATCCCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCTCCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGa                  <  1:749273/238‑1 (MQ=255)
                                                                                                                                                                                                                              acggccgcGTGAGCGGGGCGCCCGCCGCGACCAATGCAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGAc                 <  1:672314/236‑1 (MQ=255)
                                                                                                                                                                                                                                cgccgcATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACa                >  2:454951/1‑238 (MQ=255)
                                                                                                                                                                                                                                  ccgcATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCATTATTAAGCAGCTGGGTCATGACATCg             >  2:1138433/1‑239 (MQ=255)
                                                                                                                                                                                                                                   cgcATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGt            >  1:56143/1‑239 (MQ=255)
                                                                                                                                                                                                                                    ggcATGGGCAGGAGCGCCGCCGCGACCAAGGGAGCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCGTGACATCGt            <  1:1159379/238‑1 (MQ=255)
                                                                                                                                                                                                                                    ggcAGGAGCAGGGGCGCCGCCGCGCCAAGTGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGt            <  1:282279/238‑1 (MQ=255)
                                                                                                                                                                                                                                          ggCAGGGGCGCCGCCCCAACCAATGGACCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCg       <  1:947579/237‑1 (MQ=255)
                                                                                                                                                                                                                                           ggaGGTCCCCCCCCGCGCCAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCg       <  1:141469/235‑1 (MQ=255)
                                                                                                                                                                                                                                               gTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCgttgtt  <  1:741318/238‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
GGTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTACGCGCAGGCCGCCTTTCAGGCTTCCTGGCGGTTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTT  >  NC_000913/4151094‑4151565

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