Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_012971 754,386 T→A G7G (GGA→GGT ECD_RS03725 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_012971754,3860TA98.0% 171.1 / NA 49G7G (GGA→GGTECD_RS03725hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (28/20);  minor base C (1/0);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGAGGTGGGAATATCCTGCATTCCCGAACCCATCGACGAACTGTTTCAAGGCTTCTTGGACGTCGCTGGCGTGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGCTTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTACGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGCTCTCCATTCCATTCTCCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACCCGCCGTAGCGAGTTCAGATAAAATAAATCCCCGCGAGTGCGAGGATTGTTATGTAATATTGGGTTTAATCATCTATATGTTT  >  NC_012971/754096‑754512
                                                                                                                                                                                                                                                                                                  |                                                                                                                              
aCCGGAGGAGGGAATATCCTGCATTCGCGCACCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACt                                                                                                                       <  1:60919/300‑1 (MQ=255)
    gaggagGGAATATCCTGCATTCGCGCACCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGa                                                                                                                   >  1:478337/1‑300 (MQ=255)
      gtagtgAATATCTTGCATTCGCACACCCATCGACGAACGGGTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAg                                                                                                                 <  2:2968/295‑1 (MQ=255)
      gtagggAATATCCTGCATTCGCGCACCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAg                                                                                                                 <  2:500009/297‑1 (MQ=255)
      gtagggAATATCCTGCATTCGCGCACCCATCGACGAACTGTTTCAAGGCGTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAg                                                                                                                 <  2:77879/297‑1 (MQ=255)
      ggagggAATATCCTGCATTCGCGCACCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCACTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAg                                                                                                                 <  2:500060/297‑1 (MQ=255)
       gagggAATATCCTGCATTCGCGCACCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGTAAg                                                                                                                 <  2:2679/297‑1 (MQ=255)
         gggAATATCCTGCATTCGCGCACCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTAGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTATGTTGAATGGTTCTCCATTCCATTCACCTGTGACTAGGAAGTGc                                                                                                              >  2:304106/1‑300 (MQ=255)
         gggAATATCCTGCATTCGCGCACCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTACGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGc                                                                                                              >  2:304740/1‑300 (MQ=255)
          ggAATATCCTGCATTCGCGCACCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGNAGCAGTCGTTGATTTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCa                                                                                                             >  1:573547/1‑300 (MQ=255)
          ggAATATCCTGCATTCGCGCACCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGAAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTAAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAGTGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCa                                                                                                             >  2:111682/1‑300 (MQ=255)
           gAATATCCTGCATTCGCGCACCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCAt                                                                                                            <  2:518064/300‑1 (MQ=255)
                tCCTGCATTCGCGCACCCATCGACGAACTGTTTCAATGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTatc                                                                                                       <  1:122458/300‑1 (MQ=255)
                tCCTGCATTCGCGCACCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGANTATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATTTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTatc                                                                                                       <  2:573547/300‑1 (MQ=255)
                   ttcATTCGCGCACCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTatcatc                                                                                                    <  2:449851/298‑1 (MQ=255)
                   tGCATTCGCGCACCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTatcatc                                                                                                    <  1:724492/300‑1 (MQ=255)
                      attcgcgcaCCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTcc                                                                                                 <  2:44109/292‑1 (MQ=255)
                       ttcgcgcaCCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCAtctaca                                                                                                >  1:77319/8‑297 (MQ=255)
                         cgcgcaCCCATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTACATa                                                                                              >  1:691715/6‑300 (MQ=255)
                              aCCCATCGACGAACTGGTTCAAGGCTTCTTGGGCGCCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATaaaac                                                                                          <  1:284618/299‑1 (MQ=255)
                               cccATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTATATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGACGTGCATTTATCATCTACATaaaacaa                                                                                        >  2:372910/1‑300 (MQ=255)
                               cccATCGACGAACTGTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAGGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATaaaacaa                                                                                        >  1:531525/1‑300 (MQ=255)
                                             gTTTCAAGGCTTCTTGGGCGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTACATAAAACAAAACTCGCCGTAGc                                                                           >  1:713400/1‑299 (MQ=255)
                                                            gggcGTCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGataaaaga                                                           <  1:454506/297‑3 (MQ=255)
                                                                 tCGCTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTAAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATcc                                                      >  2:447278/1‑300 (MQ=255)
                                                                    cTGGCGAGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATcccc                                                    >  1:706122/1‑299 (MQ=255)
                                                                              ttCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTACATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGGAATCCCCGCGAgggc                                            >  1:279297/1‑294 (MQ=255)
                                                                              ttCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTACATAAAACAAAACTCGCAGTAGCGAGTTCAGATAAAAGAAATCCACGCGAGTGCGa                                          >  2:747790/1‑299 (MQ=255)
                                                                                   cTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATcc                                                      >  2:38017/1‑282 (MQ=255)
                                                                                   cTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGAtt                                     >  1:388150/1‑299 (MQ=255)
                                                                                         aaGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCATTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGATTGTTAttc                              <  2:648564/300‑3 (MQ=255)
                                                                                         aaGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCACCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGATTGTTAttc                              <  2:387250/300‑3 (MQ=255)
                                                                                           gtgtCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTNATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGATTGTTAttcat                            >  1:153585/1‑296 (MQ=255)
                                                                                           gtgtCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTACGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTTAGATAAAAGAAATCCCCGAGAGTGTGAGGATTGTTAttcat                            >  2:481020/1‑296 (MQ=255)
                                                                                             gtCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTAAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTACATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGGTTGTTAttcatt                           >  2:143261/1‑294 (MQ=255)
                                                                                             gtCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCAATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATTTACATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCACGCGAGTGCGAGGATTGTTAttcatt                           >  2:143246/1‑294 (MQ=255)
                                                                                                      aTCGCAAAGTCTCCGCAATTACACGCCAGACAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTACGTCTGCCTGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCCCCATTCCATTCCCCTGTGACTCGGAAGTGCACTTATCATCTACATAAAACAAAACTCGCCGTAGCGCGTTCAGATAAAAGAAATCCACGCGAGTGCGAGGATTGTTAtttattgccggtctt                  >  1:248777/1‑285 (MQ=255)
                                                                                                      aTCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGATTGTTAttcattgccgatattc                 >  1:249364/1‑285 (MQ=255)
                                                                                                      aTCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTACATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGATTGTTAttcattgccgctattc                 >  1:113050/1‑285 (MQ=255)
                                                                                                      aTCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTACATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGATTGTTAttcattgccgatattc                 >  1:113015/1‑285 (MQ=255)
                                                                                                      aTCGCAAAGTCTACGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGGTTGTTAttcattgccgatatt                  >  1:249565/1‑285 (MQ=255)
                                                                                                           aaaGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGATTGTTAttcattgccgatattcacct             <  2:555583/299‑20 (MQ=255)
                                                                                                             aGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGATTGTTAttcattgccgatattcaccttta          >  2:713412/1‑278 (MQ=255)
                                                                                                              gTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGATTGTTAttcattgccgatattcacctttat         >  1:658323/1‑277 (MQ=255)
                                                                                                                ctcCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGATTGTTAttcattgccgatattcacctttatc        >  2:468170/1‑275 (MQ=255)
                                                                                                                 tacGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATNTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGATTGTTAttcattgccgatattcacctttatcg       <  2:263751/297‑26 (MQ=255)
                                                                                                                     cAATTACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGATTGTTAttcattgccgatattcacctttatcgcgaac  >  1:480065/1‑270 (MQ=255)
                                                                                                                         tACACGCAAGAAAAAACCGCCATCAGGCGGATTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTATGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTGACTCGTAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGATTGTTAttcattgccgatattcacctttatcg       <  1:623321/291‑26 (MQ=255)
                                                                                                                                  gAAAAAACCGCCATCAGGCTGGTTGGTGTTCTTTCAGTTCTTCAATTCGAAGATTGGTTACTTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATATTCTCCGCTTCGATAACTCTGTTGAATGGTTCTCCATTCCATTCACCTGTTACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACTCGCCGTAGCGAGTTCAGATAAAAGAAATCCCCGCGAGTGCGAGGATTGTTAttcattgccgatattcacctttatcg       <  2:373281/282‑26 (MQ=255)
                                                                                                                                                                                                                                                                                                  |                                                                                                                              
ACCGGAGGTGGGAATATCCTGCATTCCCGAACCCATCGACGAACTGTTTCAAGGCTTCTTGGACGTCGCTGGCGTGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACACGCAAGAAAAAACCGCCATCAGGCGGCTTGGTGTTCTTTCAGTTCTTCAATTCGAATATTGGTTACGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATGTTCTCCGCTTCGATAACTCTGTTGAATGGCTCTCCATTCCATTCTCCTGTGACTCGGAAGTGCATTTATCATCTCCATAAAACAAAACCCGCCGTAGCGAGTTCAGATAAAATAAATCCCCGCGAGTGCGAGGATTGTTATGTAATATTGGGTTTAATCATCTATATGTTT  >  NC_012971/754096‑754512

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