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→GCC) ppcphosphoenolpyruvate 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

GGTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTT  >  NC_000913/4151094‑4151565
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
ggTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTgcg                                                                                                                                                                                                                                              <  1:1188548/239‑1 (MQ=255)
         cGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGCGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGG‑TGTGACCACAGCCCC‑GCATGAGCCGGGGCG‑CCGC‑cgcg                                                                                                                                                                                                                                    >  1:93919/1‑238 (MQ=255)
            aGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCTTTTCAGGCTTCCTGGCGGG‑TGTGACAGCAGCGCC‑GCATGAGCAGGGGCG‑CCGC‑CGCGAccc                                                                                                                                                                                                                                >  2:1079390/1‑238 (MQ=255)
            aGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CCCGACCa                                                                                                                                                                                                                                >  1:103604/1‑239 (MQ=255)
              gagaTGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TTTGACAGCAGCGCC‑GCATGAGCAGGGGCG‑CCGC‑CGCGACCA‑At                                                                                                                                                                                                                             >  1:1157237/1‑239 (MQ=255)
                   gactgacAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGGGCG‑CCGC‑CGCGACCA‑ATGGaa                                                                                                                                                                                                                         >  1:385862/1‑238 (MQ=255)
                         taGTTCCGCCATAGGGCGACGCCGGCTTTTTTTCGGCGGCGGGGGGGGCGGCAGGTTGGCTCCCAGACCCGCCCCGGTAAAAGGCGACAGGCTGCGGGGGGTTTTTTCTGGCAGACCATATTTAAAGCGGATCATCCCGCCCTTTT‑CGTTTA‑CGCGCGGCCC‑GCTTTTCAGGCTTCCGGGCGTT‑TGTGACAGCAGCGCC‑GCAGGAGCGGGGGCG‑CCGC‑CGCGCCAA‑AGGGAACcgccgc                                                                                                                                                                                                                  <  2:914727/238‑1 (MQ=255)
                                       tGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGCCCGTGGaa                                                                                                                                                                                                        >  1:583073/1‑235 (MQ=255)
                                           aCGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACCGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TTGTACAGCAGCCCC‑GCCTCAGCAGGCGCC‑GCGC‑CCGCACCA‑ATGGGACCGCCCCGCCAGCTGGACAACGtgc                                                                                                                                                                                                >  2:870513/1‑237 (MQ=255)
                                           aCGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCAAGAGCCGGTGCG‑CCGC‑CGCGGCCC‑ATGGGACCGCCGCGCCCGGGGGACAACg                                                                                                                                                                                                   >  1:521087/1‑236 (MQ=255)
                                                  ctctctTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCCACCA‑ATGGGACCGCCCCGCCCGTGGAACAACGTTAGCTcgac                                                                                                                                                                                         >  1:564336/1‑236 (MQ=255)
                                                  ctctctTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGGGGCAGCAACCCC‑GCATGAGCAGGGGCG‑CCGC‑CGCGACCA‑AAGGCACCGCCGCGACCGTGGAGCAACGGCAGCTCAAt                                                                                                                                                                                         >  2:744715/1‑239 (MQ=255)
                                                                 ggcggTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGGGGAACAACGTCACGTCAATACCCGCCTTTTCGCa                                                                                                                                                                          >  1:722216/1‑239 (MQ=255)
                                                                     ggggcggcGGGTGGGCTTCCAGACCCGCCCGGGTAGAAAGCGACGGGTTGCTGGCGTTGTTTTCTGGCAGCCTATTTTTAAGGGGTATCCTCCCGCCTTGTT‑GGGTTC‑CGCGCGGGCC‑GCCTTTCGGGCTTCCTGGGGGT‑GGGGCCGGCGGGCCC‑GCAGGAGCGGGCGCG‑CCGC‑CGCGACAA‑TTGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGt                                                                                                                                                                       <  2:103604/236‑1 (MQ=255)
                                                                            agggTTGGCTCACAGAATCGCCCGGTTAAAAAGCGACGGGCTGCTGACGGTGATTTCGGGCAGACAATATTTAAGGGGGATCTTCTCGCCCTGTT‑GGGTTC‑CGCGCGGGCC‑CCTTTCCAGCTTTCGTGGGGGT‑TGTGACAGCGGCGCG‑GCATGGGCGGGGGCG‑CCGC‑CGCGCCAA‑GTGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGAtt                                                                                                                                                                <  2:1237554/236‑1 (MQ=255)
                                                                                gtGGCCTCCAGACTCCCCCGGGTATAAGGCGACAGGTTGCGGGGGGTGTTTCCTGCCAGACTATATTTAAAGCGACTCCTCCCCCTTTGTG‑CGGTCC‑CGCGCGGGCC‑GCTTTTCGGGCTTCCGGGGGGT‑TGTGCCAGCAGCGCC‑GCAGGAGCGGGGGCG‑CCGC‑CGCGACCA‑AGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATg                                                                                                                                                            <  2:1157237/237‑1 (MQ=255)
                                                                                         aGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCtgtg                                                                                                                                                   >  1:366494/1‑238 (MQ=255)
                                                                                             tCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGa                                                                                                                                              >  2:950754/1‑239 (MQ=255)
                                                                                                               aCAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATg                                                                                                                                                            >  1:720331/1‑207 (MQ=255)
                                                                                                               aCAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGaa                                                                                                                               >  1:25981/1‑236 (MQ=255)
                                                                                                               aCAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAggt                                                                                                                            >  1:754438/1‑237 (MQ=255)
                                                                                                                cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCC‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCt                                                                                                                            >  1:147317/1‑238 (MQ=255)
                                                                                                                cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGAG‑CAGC‑CGCGACCA‑ATGGAACCGCCGCGACCGGGGAACAACGTAAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATACATCCTGTGCCTGAGTTTGGCCCCAGGAAGCt                                                                                                                            >  2:1085197/1‑238 (MQ=255)
                                                                                                                cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGC‑TGTGACAGCCGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CCCGACCA‑ATGGAACCGCCGCGACCGTGGAAAAACGTCAGCTCACTCCCCGCTTTTTCGCAGGGTTTGGTTAATGCATCCTGTGCCTGACATTCCGCCCAGGAAGCt                                                                                                                            >  2:7899/1‑238 (MQ=255)
                                                                                                                cAGGCTGCTGACGGTGATTTCTGGCAGACCATATTGAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACCGCAGCGCC‑GCATGAGCAGGGGAG‑CCGC‑AGAGACCA‑ATGGTAACGCCGCGACCGTGGGACAACGTCAGCACAATACCCCATTTTTCTCTGGATTTTATTTATGCGTCCTGGGCCTGATATAGTGCCCAGAAAGCt                                                                                                                            >  2:248115/1‑238 (MQ=255)
                                                                                                                                           aCCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCCGCCCAATACCCGCCTTTTCGCCGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGGAGCTGCCATCACTCCCTCATCTTTTGCttgg                                                                                                 >  2:1150789/1‑235 (MQ=255)
                                                                                                                                                      aGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAACGCC‑GCCTGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGGGGAACAACGTCAGCTCAATCCCCGCTTTTTCGCAGGCTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCACGAAGCTGCCATCACTCCCGCATCTTTTTCTGAGTCGGAATAGcgc                                                                                     >  2:1159379/1‑237 (MQ=255)
                                                                                                                                                          gATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGCTATTGCGCCCAGGAAGCTGCCATCACTCCCGCCTCTTTTGCTGAGTCGGAATAGCCAat                                                                                   >  1:271007/1‑237 (MQ=255)
                                                                                                                                                          gATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAatca                                                                                 >  1:316002/1‑239 (MQ=255)
                                                                                                                                                           atcatcTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGc                                                                                                                                                                                              >  2:860045/1‑129 (MQ=255)
                                                                                                                                                           atcatcTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGCGGCACAACGTCAGGTCAATACCCGCTTTTTCGCAGGGTTTTATTAATGTATCCTGGGCCTGATATTGCGCCCAGGAACCTGCTAACACTCCCGCATCTTTTCCTGAGGCGGAATAGCCAataa                                                                                 >  2:997699/1‑236 (MQ=255)
                                                                                                                                                            cccccccccccTG‑T‑CGGTTACCGCGCGGCCC‑CCCTTCCGGCCTCCCGGGCGGT‑GGTGACAGCAGCGCC‑GCGTGAGCGGGCGCC‑CCCC‑CGCGCCAA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAatcat                                                                                <  2:25981/230‑1 (MQ=255)
                                                                                                                                                                tccGCCCTG‑T‑CGGTGA‑CGGGCCGGCCGGCCGTCCGGGCTTCTGGGGGGT‑TGTGAC‑GCAGGGCCAGAGGGGGAAGGGGCG‑CGGC‑CGCGCCCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCa                                                                           <  2:719938/237‑1 (MQ=255)
                                                                                                                                                                 cgcgccTGTTCCGGTGC‑CGCGCGGGCG‑GCCTTGCAGCCTTCCGGGGGGT‑TGTGACAGCAGCGCC‑GCAGGGGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAGCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCa                                                                           <  1:121550/235‑1 (MQ=255)
                                                                                                                                                                     ccTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGGC‑GCATGAGCGGGGGCG‑CCGC‑CGCGGCCA‑AAGGAACCGGCGTGACCGTGGAACAACGTGAGCTCAATAACCGCTTTTTCGCAGGTTTTGGTTAAGGCCTGCTGTGCCTGATATTGCGCCCAGGAAGCTGCCCTCACTCCCGCATCTTTTGCTGATTCGGAATAGCCCATCATCAACATCTGt                                                                      >  2:89008/1‑239 (MQ=255)
                                                                                                                                                                         tt‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCCTCTTTTGCTGGGTCGGAATAGCCAATCATCACCATCTGTTTg                                                                   >  1:484350/1‑238 (MQ=255)
                                                                                                                                                                            gcgttaccgcgcgGGCC‑GCTTTACGGGCTCCCGGGCGGT‑TGTGACAGCAGCGCC‑ACATGAGCGGGTGCG‑CCGC‑CGCGACAA‑AGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGcc                                                                 <  1:930179/232‑1 (MQ=255)
                                                                                                                                                                               ttA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGa                                                             >  1:843362/1‑239 (MQ=255)
                                                                                                                                                                                    ggccgggCC‑GCCTTCCAGGCTCCGGGGCGGT‑GTTGACACCCGCCGC‑GCGTGAGCGGGGGCC‑CCCC‑CGCACCCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCa                                                         <  2:883544/234‑1 (MQ=255)
                                                                                                                                                                                     cgcAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCC‑ATGGAACCGCCGCGACCGTGGAACAACGTCCGCTCAATACCCGCTTTTTCGCCGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCCCTCCCGCCTCTTTTGTTGAGTCGGAATAACCAATCATCCCCATCTGTTTGCCCTGAATccg                                                        >  2:706456/1‑237 (MQ=255)
                                                                                                                                                                                     cgcAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTTCTGAGTCGGAATAGCCAATCATCCCCCTCTGTTTGCCCTGAATCAg                                                        >  2:282279/1‑239 (MQ=255)
                                                                                                                                                                                        aGGCC‑GCCTTTCGGGCTTCCTGGCGGT‑TGTGACGGCAGCGCC‑GCAGGAGCGGGGGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGc                                                      <  1:1068068/238‑1 (MQ=255)
                                                                                                                                                                                          gcc‑gccTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCAc                                                   >  2:741318/1‑239 (MQ=255)
                                                                                                                                                                                          gcc‑gccTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCAc                                                   >  1:883544/1‑239 (MQ=255)
                                                                                                                                                                                           gtccccTTTCCTGCGTCCCGGTGGTGTGTGACGGGAGCGCC‑GCATGGGCGGCTCCG‑CCCC‑CGCGCACA‑TAGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCa                                                    <  2:564336/235‑1 (MQ=255)
                                                                                                                                                                                               tcttcagggCTTCCTGGGGGT‑GGGCACAGCGGCGCC‑GGAGGAGCGGGGGCG‑CCGC‑CGCGACAA‑AGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa                                               <  2:271007/232‑1 (MQ=255)
                                                                                                                                                                                               ccTTTCAGGCTTTCGGGCGGT‑TGTGACGACAGCGGC‑GTGGGGGCAGGTGCG‑CCCC‑CGCAACCA‑GTGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa                                               <  1:283325/239‑1 (MQ=255)
                                                                                                                                                                                               ccTTTCAGGCTTCCTGGCGGT‑TGTGACAGCCGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTTCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGAt                                                >  1:969262/1‑238 (MQ=255)
                                                                                                                                                                                                cgtccAGGCTTCCGGGGGGT‑TGCAGCAGCCGCCCC‑GCAGGAGCAGGGGCC‑CCCC‑CGCGCCAA‑ATGGACCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATa                                               <  2:366494/234‑1 (MQ=255)
                                                                                                                                                                                                    cAGCCTTCCTGGTGGT‑TGTGCCAGGAGCGCC‑GCAGGAGCAGGGCCC‑CCCC‑CGCGCCCA‑AGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGt                                          <  2:316002/239‑1 (MQ=255)
                                                                                                                                                                                                                    tgtggACAGCAGCCCC‑GCATGAGCAGGGCCC‑CCGC‑CGCAACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGagcagc                            <  1:747851/234‑1 (MQ=255)
                                                                                                                                                                                                                              ggcgcC‑ACAGGAGCGGGGGCG‑CCGC‑GGCGACAA‑AGGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGAc                 <  2:266304/238‑1 (MQ=255)
                                                                                                                                                                                                                              cgcgcG‑GCAGGAGCAGG‑GCG‑CCCCACGGGGCCATATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATCCCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCTCCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGa                  <  1:749273/238‑1 (MQ=255)
                                                                                                                                                                                                                              acggcc‑gcGTGAGC‑GGGGCGCCCGC‑CGCGACCA‑ATGCAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGAc                 <  1:672314/236‑1 (MQ=255)
                                                                                                                                                                                                                                cgcc‑gcATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACa                >  2:454951/1‑238 (MQ=255)
                                                                                                                                                                                                                                  cc‑gcATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCATTATTAAGCAGCTGGGTCATGACATCg             >  2:1138433/1‑239 (MQ=255)
                                                                                                                                                                                                                                   c‑gcATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGt            >  1:56143/1‑239 (MQ=255)
                                                                                                                                                                                                                                    ggcATGGGCAGGAGCG‑CCGC‑CGCGACCA‑AGGGAGCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCGTGACATCGt            <  1:1159379/238‑1 (MQ=255)
                                                                                                                                                                                                                                    ggcAGGAGCAGGGGCG‑CCGC‑CGCGCCAA‑GTGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGt            <  1:282279/238‑1 (MQ=255)
                                                                                                                                                                                                                                          ggCAGGGGCG‑CCGC‑CCCAACCA‑ATGGACCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCg       <  1:947579/237‑1 (MQ=255)
                                                                                                                                                                                                                                           ggaGGTCCC‑CCCC‑CGCGCCAA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCg       <  1:141469/235‑1 (MQ=255)
                                                                                                                                                                                                                                               gTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCgttgtt  <  1:741318/238‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
GGTAGACATCGCAGGAGATGACTGACAGTTCATCCATAATGCGACGCCAGCTCTCTTTCGGCTCCGGCGGTGGCAGCAGGTTGGCTTCCAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTT‑CGGTTA‑CGCGCAGGCC‑GCCTTTCAGGCTTCCTGGCGGT‑TGTGACAGCAGCGCC‑GCATGAGCAGGTGCG‑CCGC‑CGCGACCA‑ATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAATAGCCAATCATCACCATCTGTTTGCCCTGAATCAGGCCACGATACCAGTCAATATTGAGCAGCTGGGTCATGACATCGTTGGCGTTGTT  >  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: ‑