Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 401,924 A→G 5.8% V201A (GTC→GCC)  PP_0335 ← glutathione S‑transferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947401,9240AG 5.8% 200.2 / 2.2 69V201A (GTC→GCC) PP_0335glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base A (33/32);  new base G (2/2);  total (35/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 = 8.06e-01

CACATAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCA  >  NC_002947/401776‑402068
                                                                                                                                                    |                                                                                                                                                
ccaaTAATGGGGTGTCAGGCACCGGCTTCGGCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGa                                                                                                                                                  <  8:87956/146‑1 (MQ=255)
 acaTAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGAc                                                                                                                                                 >  3:72795/1‑149 (MQ=255)
   aTAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCa                                                                                                                                               >  2:323920/1‑149 (MQ=255)
    tAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCAt                                                                                                                                              >  6:239496/1‑149 (MQ=255)
     aaTGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATc                                                                                                                                             >  5:278836/1‑149 (MQ=255)
     aaTGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATc                                                                                                                                             >  8:6964/1‑149 (MQ=255)
      aTGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATc                                                                                                                                             >  8:68821/1‑148 (MQ=255)
          ggTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATgg                                                                                                                                        <  6:295125/149‑1 (MQ=255)
           gtgtCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGGCCATCGATggt                                                                                                                                       >  5:279863/1‑148 (MQ=255)
              gcccgCACCGGCTTCGCCGGTGTGCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGGCCATCGATGGTCGc                                                                                                                                    <  6:279863/145‑1 (MQ=255)
                aGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCtt                                                                                                                                  <  6:278836/149‑1 (MQ=255)
                aGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCtt                                                                                                                                  >  6:311084/1‑149 (MQ=255)
                 ggCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGTGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAATGGGCCCTCAGGCCACTGGCGCTGTAGCGGCCCTCGATGGTCGcga                                                                                                                                  >  6:51743/1‑146 (MQ=38)
                 ggCACCGGCTTCGCCGGTGTTCGCGGGCACGACCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGCCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCtt                                                                                                                                  >  6:105664/1‑148 (MQ=255)
                  gCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGTGCTGTAGCGACCATCGATGGCCGCTTGa                                                                                                                                <  6:133350/149‑1 (MQ=255)
                   cACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGTGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGGCCATCGATGGTCGCtacag                                                                                                                               <  5:51743/149‑5 (MQ=255)
                    aCCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGac                                                                                                                               >  3:342378/1‑148 (MQ=255)
                         cTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTc                                                                                                                         <  6:298200/149‑1 (MQ=255)
                               cGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATa                                                                                                                   <  4:72795/149‑1 (MQ=255)
                                        cgGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCAcgccgc                                                                                                          <  5:239496/149‑1 (MQ=255)
                                               gCCCGCTCCCAAAAGCCCACCACAGATCGCGAACACTGTGTAAAACCTTGTCGGAGCGGGTATACCCGCGAAGGGGCCCTCCGGCCACTGGCGCTGTAGCGACCATCGAt                                                                                                                                          <  5:252246/110‑1 (MQ=255)
                                               gCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCTTCGAt                                                                                                                                          >  6:252246/1‑110 (MQ=255)
                                               gCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGAt                                                                                                                                          <  5:252247/110‑1 (MQ=255)
                                               gCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGAt                                                                                                                                          >  6:252247/1‑110 (MQ=255)
                                                       ccAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGgcg                                                                                           >  8:157307/1‑149 (MQ=255)
                                                       ccAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGgc                                                                                            <  3:247175/148‑1 (MQ=255)
                                                       ccAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGgc                                                                                            <  5:311084/148‑1 (MQ=255)
                                                        cAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGgcgc                                                                                          >  6:359564/1‑149 (MQ=255)
                                                              cccACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACgc                                                                                    >  2:192970/1‑149 (MQ=255)
                                                                caccacAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACgcca                                                                                  <  7:68821/149‑1 (MQ=255)
                                                                caccacAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACgcca                                                                                  <  3:59405/149‑1 (MQ=255)
                                                                  ccacAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACgccagc                                                                                <  7:50425/149‑1 (MQ=255)
                                                                   cacaGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACgccagcc                                                                               <  1:323920/149‑1 (MQ=255)
                                                                    acaGATAGCGAACACTGTGGAAAACCTTGTAGGAGTGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACgccagcc                                                                               >  6:353137/1‑148 (MQ=255)
                                                                          aGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGTCGCTTGACCACCTCGGCATACCACGCCGCCAACCCa                                                                                                   <  4:4463/122‑1 (MQ=255)
                                                                          aGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGTCGCTTGACCACCTCGGCATACCACGCCGCCAACCCa                                                                                                   >  3:4463/1‑122 (MQ=255)
                                                                           gCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCg                                                                        <  6:213802/148‑1 (MQ=255)
                                                                              aaCACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCt                                                                    >  2:127396/1‑149 (MQ=255)
                                                                                 acTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGCGCTGAc                                                                 >  3:145857/1‑149 (MQ=255)
                                                                                     tgGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCa                                                              >  3:163637/1‑148 (MQ=255)
                                                                                            ccTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATc                                                          >  5:75521/1‑145 (MQ=255)
                                                                                                 tAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGc                                                  <  4:145857/148‑1 (MQ=255)
                                                                                                     aGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAgc                                             >  2:140684/1‑149 (MQ=255)
                                                                                                      gCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCa                                                        <  6:152848/137‑1 (MQ=255)
                                                                                                      gCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCa                                                        >  5:152848/1‑137 (MQ=255)
                                                                                                         ggTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGAt                                                     <  8:125378/137‑1 (MQ=255)
                                                                                                         ggTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGAt                                                     >  7:125378/1‑137 (MQ=255)
                                                                                                          gTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCg                                                                        <  2:184345/117‑1 (MQ=255)
                                                                                                          gTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCg                                                                        >  1:184345/1‑117 (MQ=255)
                                                                                                              aCCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGc                                     <  5:359564/148‑1 (MQ=255)
                                                                                                               cccGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACg                                                                >  4:211227/1‑120 (MQ=255)
                                                                                                               cccGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACg                                                                <  3:211227/120‑1 (MQ=255)
                                                                                                                ccGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACgc                                                                                    >  5:226653/1‑99 (MQ=255)
                                                                                                                ccGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACgc                                                                                    <  6:226653/99‑1 (MQ=255)
                                                                                                                  gcgAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGc                                >  3:11113/1‑149 (MQ=255)
                                                                                                                           cccTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATcgccg                       >  2:285048/1‑149 (MQ=255)
                                                                                                                                aGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAgg                  <  6:59453/149‑1 (MQ=255)
                                                                                                                                 ggCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGt                 >  5:59004/1‑149 (MQ=255)
                                                                                                                                  gCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTc                >  6:363601/1‑149 (MQ=255)
                                                                                                                                     aCTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGaa             <  6:75521/149‑1 (MQ=255)
                                                                                                                                         gcgcTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCg         <  4:163637/149‑1 (MQ=255)
                                                                                                                                           gcTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGga       <  7:157307/149‑1 (MQ=255)
                                                                                                                                           gcTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCAAAGGCCACGCCAATCGCCGCCAGGTCGAAGCCg         >  8:18888/1‑147 (MQ=255)
                                                                                                                                             tGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGgaga     <  4:342378/149‑1 (MQ=255)
                                                                                                                                               tAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGgaga     <  7:18888/147‑1 (MQ=255)
                                                                                                                                               tAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGgaga     <  1:127396/147‑1 (MQ=255)
                                                                                                                                               tAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGAt    >  8:20510/1‑148 (MQ=255)
                                                                                                                                                 gCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCa  <  1:140684/148‑1 (MQ=255)
                                                                                                                                                   gACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCgg        >  3:60860/1‑140 (MQ=255)
                                                                                                                                                   gACCATCGATGGCCGCTTGACCACCTCGGCATACCACGACGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCgg        <  4:60860/140‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CACATAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCA  >  NC_002947/401776‑402068

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: ‑