Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,344,246 A→C 8.6% K496Q (AAG→CAG) ‡ pgi‑II → glucose 6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,344,2460AC 8.6% 157.5 / 7.9 58K496Q (AAG→CAG) ‡pgi‑IIglucose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base A (25/28);  new base C (2/3);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01

GCTAAAGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGAT  >  NC_002947/5344108‑5344394
                                                                                                                                          |                                                                                                                                                    
gcTAAAGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAg                                                                                                                                                    >  7:253385/1‑141 (MQ=32)
gcTAAAGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAg                                                                                                                                                    <  8:253385/141‑1 (MQ=32)
    aaGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAACACCTTGTGCTCGTAc                                                                                                                                          <  5:286689/147‑1 (MQ=11)
    aaGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAACACCTTGTGCTCGTAc                                                                                                                                          <  8:323910/147‑1 (MQ=11)
    aaGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAACACCTTGTGCTCGTAc                                                                                                                                          <  8:319506/147‑1 (MQ=11)
       gCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGc                                                                                                                                     <  3:304830/149‑1 (MQ=32)
              cGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCt                                                                                                                              >  7:158601/1‑149 (MQ=32)
                  gCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTgggg                                                                                                                          <  3:332601/149‑1 (MQ=32)
                    aGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGTTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGACAGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCTTTATCTGGGGTa                                                                                                                        <  6:104678/149‑1 (MQ=255)
                    aGACAGAGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTa                                                                                                                        <  8:151418/149‑1 (MQ=17)
                     gACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTAt                                                                                                                       <  8:158601/149‑1 (MQ=32)
                      aCATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATc                                                                                                                      <  2:238312/149‑1 (MQ=32)
                        aTCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTTTCTGGGGTATCaa                                                                                                                    >  3:255703/1‑149 (MQ=17)
                          cGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACg                                                                                                                  >  5:188542/1‑149 (MQ=32)
                           gAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGc                                                                                                                 >  4:218162/1‑149 (MQ=32)
                             aaaGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCt                                                                                                               <  2:313010/149‑1 (MQ=32)
                                 cTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGa                                                                                                           >  8:262497/1‑149 (MQ=32)
                                      cccGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACaa                                                                                                                                                     >  4:65643/1‑102 (MQ=32)
                                      cccGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACaa                                                                                                                                                     <  3:65643/102‑1 (MQ=32)
                                      cccGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGt                                                                                                      >  1:189126/1‑149 (MQ=32)
                                      cccGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGt                                                                                                      >  7:278303/1‑149 (MQ=32)
                                         gCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTggg                                                                                                   <  7:269069/149‑1 (MQ=32)
                                         gCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGATGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTggg                                                                                                   >  1:236158/1‑149 (MQ=17)
                                           acaAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGc                                                                                                 >  2:256885/1‑149 (MQ=32)
                                             aaGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTggg                                                                                                   >  7:256916/1‑145 (MQ=32)
                                             aaGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTggg                                                                                                   <  8:256916/145‑1 (MQ=32)
                                              aGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTc                                                                                              >  3:320263/1‑149 (MQ=32)
                                                gTGATCCCGGGTAACCGCCCGAGCAACACGATGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCga                                                                                            <  2:236158/149‑1 (MQ=17)
                                                   aTCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCgagct                                                                                         <  1:256885/149‑1 (MQ=32)
                                                      ccGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCgagctgg                                                                                       <  7:262497/148‑1 (MQ=32)
                                                      ccGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCAGCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCgagctgg                                                                                       <  4:235558/148‑1 (MQ=17)
                                                                gCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCgagctgggcaaggagctg                                                                            <  7:301281/149‑1 (MQ=34)
                                                                gCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCgagctgggcaaggagctg                                                                            >  5:344496/1‑149 (MQ=34)
                                                                    gAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGt                                                                                                                                            <  1:173118/81‑1 (MQ=33)
                                                                    gAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGt                                                                                                                                            >  2:173118/1‑81 (MQ=33)
                                                                     aGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCgagctgggcaaggagctgggc                                                                         >  2:307181/1‑147 (MQ=34)
                                                                     aGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCgagctgggcaaggagctgggc                                                                         >  2:57188/1‑147 (MQ=34)
                                                                          cacGTTGGTGGTCGAGCGCATCAACCCGCGCCGGCTGGGTGCGCTCGAACACAAATAAGAACACAAAGTTTTAGTACATATAATTATATAAAAAATAAACGCCTTCGATCATTTGGGCGTCGAGATGGGAAAGGAGCTGGGAAAGGGct                                                                  >  8:226708/1‑148 (MQ=2)
                                                                             gTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGttt                                                               <  6:344496/149‑1 (MQ=35)
                                                                                ggtggtCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGAACAAGGTGTTCATACAGAGCGTAATCTGGGGTATCAACGCCTTCGATCAGTTGGGCGTCTAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTAc                                                             <  6:341756/148‑1 (MQ=11)
                                                                                        aGCGCATCAGCCCGAGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTg                                                    <  7:226708/149‑1 (MQ=21)
                                                                                          cgcATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGt                                                  <  1:57188/149‑1 (MQ=34)
                                                                                                aGCCCGCGCCGGCTGGGTGCGCTGGTGTCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTa                                             <  6:64922/148‑1 (MQ=21)
                                                                                                aGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTa                                             <  3:212085/148‑1 (MQ=34)
                                                                                                 gCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGc                                           >  1:60305/1‑149 (MQ=34)
                                                                                                 gCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGc                                           >  5:212276/1‑149 (MQ=34)
                                                                                                   ccGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAgcgc                                                       >  7:89736/1‑135 (MQ=34)
                                                                                                   ccGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAgcgc                                                       <  8:89736/135‑1 (MQ=34)
                                                                                                     gcgcCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTg                                        >  5:40512/1‑148 (MQ=34)
                                                                                                                gTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGcc                            <  6:188542/149‑1 (MQ=34)
                                                                                                                 tGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGcc                            >  8:293970/1‑148 (MQ=34)
                                                                                                                 tGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCg                           <  8:278303/149‑1 (MQ=34)
                                                                                                                   cgcTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGaa                         <  4:320263/149‑1 (MQ=34)
                                                                                                                    gcTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCg                           >  5:22656/1‑146 (MQ=34)
                                                                                                                        gtggCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACgg                    >  2:200399/1‑149 (MQ=34)
                                                                                                                           gctctGTACGAACACCTTGTGCTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTc                 >  5:337403/5‑149 (MQ=11)
                                                                                                                           gCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCTTGGAAGACAGCGCCGAAGACGGtt                  <  8:182449/148‑1 (MQ=21)
                                                                                                                              aTGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCAc              <  6:22656/149‑1 (MQ=34)
                                                                                                                               tGTACGAACACCTTGTGCTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCAcc             >  6:108413/1‑149 (MQ=11)
                                                                                                                                      gCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGcc      >  1:38373/1‑149 (MQ=34)
                                                                                                                                       cacaAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGCCAGCTCCGAAGACGGTTCCACCCAAGGCCt     >  4:205616/1‑149 (MQ=21)
                                                                                                                                       cacaAGGTGTTCGTACAGAGCGATATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGcc      >  8:38561/1‑148 (MQ=21)
                                                                                                                                       cacaAGGTGTTCGTACAGAGCGATATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGcc      <  7:38561/148‑1 (MQ=21)
                                                                                                                                        acaAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTg    <  2:38373/149‑1 (MQ=34)
                                                                                                                                         caAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGAGCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTATACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCCCCCACGGCCTGa   >  8:5368/1‑149 (MQ=255)
                                                                                                                                          aaGGTGTTCTTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTATACCAGCGCTTTGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGAt  <  7:5368/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
GCTAAAGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCCTGAT  >  NC_002947/5344108‑5344394

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