Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,344,242 G→A 10.1% E494E (GAG→GAA 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,2420GA10.1% 156.7 / 6.5 59E494E (GAG→GAApgi‑IIglucose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base A (3/3);  total (27/32)
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

CGAATTGCGCGCTAAAGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCC  >  NC_002947/5344098‑5344390
                                                                                                                                                |                                                                                                                                                    
cgaaTTGCGCGCTAAAGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAacacc                                                                                                                                                  >  6:286689/1‑148 (MQ=17)
          gcTAAAGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAg                                                                                                                                                >  7:253385/1‑141 (MQ=32)
          gcTAAAGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAg                                                                                                                                                <  8:253385/141‑1 (MQ=32)
              aaGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAACACCTTGTGCTCGTAc                                                                                                                                      <  8:319506/147‑1 (MQ=11)
              aaGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAACACCTTGTGCTCGTAc                                                                                                                                      <  5:286689/147‑1 (MQ=11)
              aaGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAACACCTTGTGCTCGTAc                                                                                                                                      <  8:323910/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)
                                               ccccGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAg                                                                                                                                                      >  6:133372/1‑98 (MQ=33)
                                                cccGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACaa                                                                                                                                                 >  4:65643/1‑102 (MQ=32)
                                                cccGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACaa                                                                                                                                                 <  3:65643/102‑1 (MQ=32)
                                                cccGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGt                                                                                                  >  7:278303/1‑149 (MQ=32)
                                                cccGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGt                                                                                                  >  1:189126/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                                                                                                                                        >  2:173118/1‑81 (MQ=33)
                                                                              gAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGt                                                                                                                                        <  1:173118/81‑1 (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                                                   <  8:89736/135‑1 (MQ=34)
                                                                                                             ccGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAgcgc                                                   >  7:89736/1‑135 (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)
                                                                                                                                                |                                                                                                                                                    
CGAATTGCGCGCTAAAGGCCTCAACGAAGCAGACATCGAAAAGCTGGCCCCGCACAAGGTGATCCCGGGTAACCGCCCGAGCAACACGTTGGTGGTCGAGCGCATCAGCCCGCGCCGGCTGGGTGCGCTGGTGGCCATGTACGAGCACAAGGTGTTCGTACAGAGCGTTATCTGGGGTATCAACGCCTTCGATCAGTGGGGCGTCGAGCTGGGCAAGGAGCTGGGCAAGGGCGTTTACCAGCGCTTGGTCGGTAGCCTGGAAGACAGCGCCGAAGACGGTTCCACCCAAGGCC  >  NC_002947/5344098‑5344390

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: