Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,526,411 A→C 9.3% K219Q (AAG→CAG) ‡ PP_2215 → acetyl‑CoA acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,526,4110AC 9.3% 249.6 / 12.5 84K219Q (AAG→CAG) ‡PP_2215acetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (45/31);  new base C (5/3);  total (50/34)
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

TGCCCAGGCCAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTCCGAGGC  >  NC_002947/2526263‑2526559
                                                                                                                                                    |                                                                                                                                                    
tgCCCAGGCCAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGtcatca                                                                                                                                                      <  4:69358/149‑1 (MQ=255)
  cccAGGCCAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAg                                                                                                                                                    <  1:358785/149‑1 (MQ=255)
   ccAGGCCAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAg                                                                                                                                                    >  6:176731/1‑148 (MQ=255)
        ccAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGAc                                                                                                                                              <  5:106678/149‑1 (MQ=255)
          aaTGCCTTCAGCCGCGAGGCCCAGGACCAGTGCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAATCGCGTCATCAAGGATGACGa                                                                                                                                            <  6:113903/149‑1 (MQ=255)
                ttCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCTCGTCATCCTTGATGACgcgctttt                                                                                                                                      >  6:353170/1‑145 (MQ=255)
                ttCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAgcc                                                                                                                                      >  4:264512/1‑149 (MQ=255)
                 tCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAgccg                                                                                                                                     >  2:167832/1‑149 (MQ=255)
                 tCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAgccg                                                                                                                                     >  2:173247/1‑149 (MQ=255)
                  cAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAgccg                                                                                                                                     <  2:346552/148‑1 (MQ=255)
                              ccAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGa                                                                                                                                            <  5:14717/129‑1 (MQ=255)
                              ccAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGa                                                                                                                                            >  6:14717/1‑129 (MQ=255)
                              ccAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTc                                                                                                                        >  8:49621/1‑149 (MQ=255)
                              ccAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTc                                                                                                                        >  3:296035/1‑149 (MQ=255)
                               cAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCt                                                                                                                       <  1:173247/149‑1 (MQ=255)
                               cAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCt                                                                                                                       <  1:167832/149‑1 (MQ=255)
                               cAGGACCAGTTAGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGACGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCt                                                                                                                       >  3:236714/1‑149 (MQ=255)
                                aGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGa                                                                                                                                            <  6:152548/127‑1 (MQ=255)
                                aGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGa                                                                                                                                            >  5:152548/1‑127 (MQ=255)
                                  gACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATAGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGa                                                                                                                    <  5:11621/149‑1 (MQ=255)
                                    ccAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGAc                                                                                                                   >  5:71643/1‑148 (MQ=255)
                                     cAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACaa                                                                                                                 <  5:282561/149‑1 (MQ=255)
                                      aGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCTCGTCATCCTTGATGACgcgc                                                                                                                                          >  3:85815/1‑123 (MQ=255)
                                      aGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCTCGTCATCCTTGATGACgcgc                                                                                                                                          <  4:85815/123‑1 (MQ=255)
                                        ttCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGAt                                                                                                              >  6:177726/1‑149 (MQ=255)
                                        ttCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGAt                                                                                                              >  1:76365/1‑149 (MQ=255)
                                        ttCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGAt                                                                                                              >  6:310951/1‑149 (MQ=255)
                                            ccATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTgcg                                                                                                          >  3:56409/1‑149 (MQ=255)
                                             cATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTgcg                                                                                                          >  7:95852/1‑148 (MQ=255)
                                               tCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAg                                                                                                       <  6:71643/149‑1 (MQ=255)
                                                      gCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCaaa                                                                                                 <  3:330655/148‑1 (MQ=255)
                                                       cTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAAcc                                                                                               >  5:144045/1‑149 (MQ=255)
                                                            ccGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAAc                                                                                                <  4:201298/143‑1 (MQ=255)
                                                            ccGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAAc                                                                                                >  3:201298/1‑143 (MQ=255)
                                                                    aGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTg                                                                                    >  5:81197/1‑147 (MQ=255)
                                                                       aaGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAggc                                                                               >  5:630/1‑149 (MQ=255)
                                                                           cATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCa                                                                           >  8:203910/1‑149 (MQ=255)
                                                                               agcagcGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGt                                                                       >  8:89120/1‑149 (MQ=255)
                                                                               agcagcGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGt                                                                       <  3:264512/149‑1 (MQ=255)
                                                                                gcagcGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGt                                                                       >  3:173940/1‑148 (MQ=255)
                                                                                  agcGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGt                                                                                     >  6:61218/1‑132 (MQ=255)
                                                                                  agcGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGt                                                                                     <  5:61218/132‑1 (MQ=255)
                                                                                      gCCGTTTAGCCGCCGAGATCGTGCCGGTGGAAGGCACCGAGGGCAAGGAAAAGGGCGTCATGCAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGGTTCGTGAAGGCGGCACCGTGAcggcg                                                                <  6:179467/149‑1 (MQ=255)
                                                                                       ccGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGAcggcgg                                                               <  4:56409/149‑1 (MQ=255)
                                                                                       ccGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGAcggcgg                                                               <  6:630/149‑1 (MQ=255)
                                                                                       ccGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGAcggcgg                                                               <  7:49621/149‑1 (MQ=255)
                                                                                        cggtttGCCGCCGAGATCGTGCCGGTGGGAGTCCCCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGAGAAGATTGGGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACggcggc                                                              >  6:295078/4‑149 (MQ=255)
                                                                                               ccgccgAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGc                                                        >  7:105155/1‑148 (MQ=255)
                                                                                               ccgccgAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGc                                                        <  2:76365/148‑1 (MQ=255)
                                                                                                cgccgAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGAc                                                                                                                                              >  6:183493/1‑61 (MQ=255)
                                                                                                cgccgAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGAc                                                                                                                                              <  5:183493/61‑1 (MQ=255)
                                                                                                cgccgAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCCCCGTGACGGCGGCCAACGc                                                        >  6:178552/1‑147 (MQ=255)
                                                                                                cgccgAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCa                                                      <  4:296035/149‑1 (MQ=255)
                                                                                                   cgAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATGAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGAcggcg                                                                >  6:80462/1‑136 (MQ=255)
                                                                                                   cgAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGAcggcg                                                                <  5:80462/136‑1 (MQ=255)
                                                                                                   agagaTCGTGCCGGTGGAAGTCACCTAGGGAAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCTCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAACCCGGCGTTTCGTGAAGGCGGCACCGTGCCGGCGGCCAAGGCCAGtt                                                   >  8:177568/2‑149 (MQ=255)
                                                                                                           tGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCACCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTcc                                           >  8:280897/1‑149 (MQ=255)
                                                                                                                    aaGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgct                                  >  3:168271/1‑149 (MQ=255)
                                                                                                                      gTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGcc                                                             <  2:323039/120‑1 (MQ=255)
                                                                                                                      gTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGcc                                                             >  1:323039/1‑120 (MQ=255)
                                                                                                                       tCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgctgcg                               >  1:67310/1‑149 (MQ=255)
                                                                                                                         aCCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgctgcggc                             >  7:85740/1‑149 (MQ=255)
                                                                                                                         aCCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgctgcggc                             >  8:144331/1‑149 (MQ=255)
                                                                                                                             aGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgctgcggcgctg                         <  7:203910/149‑1 (MQ=255)
                                                                                                                              gggCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTgg                        <  6:144045/149‑1 (MQ=255)
                                                                                                                              gggCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTgg                        <  4:168271/149‑1 (MQ=255)
                                                                                                                               ggCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTgg                        >  1:307266/1‑148 (MQ=255)
                                                                                                                                 cAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTg                                                                      <  8:58349/100‑1 (MQ=255)
                                                                                                                                 cAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTg                                                                      >  7:58349/1‑100 (MQ=255)
                                                                                                                                  aacgAAATGGGCGTCATCAAGGATGACGAGCAGGCGCCCAAGGCGCGGCTGGACAAGATTGCGCAGCTCAAACCGGGGTTTCCTGAAGGCGGCACCGTGACGGCGTCCAACGCCAGTTCGATTTCCGACGCCGCTGGGCCGCTGGTACt                    <  6:5654/146‑1 (MQ=255)
                                                                                                                                  aaGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGTGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGGGCTGGTACt                    >  4:229937/1‑149 (MQ=255)
                                                                                                                                  aaGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACt                    >  2:293646/1‑149 (MQ=255)
                                                                                                                                  aaGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACt                    >  5:141524/1‑149 (MQ=255)
                                                                                                                                  aaGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACt                    >  7:279729/1‑149 (MQ=255)
                                                                                                                                  aaGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACt                    >  3:317500/1‑149 (MQ=255)
                                                                                                                                  aaGGAAAAGCGCGTCAACAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCCATTTCAGACGGCGCTGCGCCGCTGGTACt                    >  5:5654/1‑149 (MQ=255)
                                                                                                                                       gctgctcgTCATCCTTGATGACGCGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATgc               >  4:348047/7‑149 (MQ=255)
                                                                                                                                       aaaGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATgc               <  7:89120/149‑1 (MQ=255)
                                                                                                                                         tgctcgTCATCCTTGATGACGAGCAGCCGCCCACGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCCACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATgcgc             >  6:247598/5‑149 (MQ=255)
                                                                                                                                          tcgcgTCATCAAGGATGACGAGCAGCCGCCCAAGCCGCGTCTGGACAAGACTGCGCAGCTCAACCCGGAGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACTGCGCTGCGGCGCTGGTACTGATGcgcc            <  6:13066/148‑1 (MQ=255)
                                                                                                                                           ctcgTCATCCTTGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTc                                                  <  3:87717/108‑1 (MQ=255)
                                                                                                                                           ctcgTCATCCTTGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTc                                                  >  4:87717/3‑110 (MQ=255)
                                                                                                                                            gcgTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGcgccg           >  8:288755/1‑148 (MQ=255)
                                                                                                                                              gTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGcgccgc          >  7:201795/1‑147 (MQ=255)
                                                                                                                                              gTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGcgccgc          <  8:201795/147‑1 (MQ=255)
                                                                                                                                                 atcaAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGTGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTcc       >  2:240636/1‑147 (MQ=255)
                                                                                                                                                  tcaAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTCCGAg    >  6:39765/1‑149 (MQ=255)
                                                                                                                                                   caAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTCCGAgg   >  4:3359/1‑149 (MQ=255)
                                                                                                                                                    aaGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTCCGAGGc  <  7:144331/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCCCAGGCCAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTCCGAGGC  >  NC_002947/2526263‑2526559

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