Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 401,938 C→T 6.9% R197R (CGG→CGA 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,9380CT 6.9% 215.3 / 4.7 73R197R (CGG→CGAPP_0335glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base C (27/40);  new base T (3/2);  total (30/43)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGT  >  NC_002947/401792‑402083
                                                                                                                                                  |                                                                                                                                                 
tgcctgacaccccactatgtgtgttcgcgggcaAGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGACGc                                                                                                                                                <  1:291840‑M2/117‑1 (MQ=255)
gtCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGCACCACATGAAGGAGCCGGCCTCCCCCCCAAGTCTCCCTCACGCCACTGGCGCTGTAGCGACCATCGATGGCCGc                                                                                                                                                <  1:1096651/150‑1 (MQ=25)
 tCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCt                                                                                                                                               >  2:507815/1‑150 (MQ=255)
  cAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCtt                                                                                                                                              <  1:660804/150‑1 (MQ=255)
     gCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGACCACCGTGGAAAACCTTCTATGATCGGGTTTACCCGCGAAGGTGCCCTCAGTGCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGac                                                                                                                                           <  2:161872/150‑1 (MQ=38)
     ccACCGGCTTCCCCGGTGTTCGCGGGCACGCCCGCTCCCAACAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGac                                                                                                                                           <  2:95391/149‑1 (MQ=255)
             ttCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACCGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTCCCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCgg                                                                                                                                   <  2:1025181/150‑1 (MQ=255)
             ttCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGCACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGAGCCTCAGTGCGCTGATGCTTTAGCGACCATCGATGGCCGCTTGACCACCTCgg                                                                                                                                   <  1:740016/150‑1 (MQ=255)
             ttCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCgg                                                                                                                                   >  1:112482/1‑150 (MQ=255)
                  cGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGTAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATAc                                                                                                                              <  1:1237466/150‑1 (MQ=255)
                       ggagcgggcgTGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGaccacc                                                                                                                                       >  2:544975‑M2/11‑136 (MQ=255)
                       ggagcgggcgTGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGaccacc                                                                                                                                       <  1:544975‑M2/126‑1 (MQ=255)
                          gcgggcaCGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCAcgccgc                                                                                                                      <  1:1197206‑M2/143‑1 (MQ=255)
                           cgggcaCGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACgccgcc                                                                                                                     <  1:1196842‑M2/144‑1 (MQ=255)
                           cgggcaCGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACgccgcc                                                                                                                     >  2:183117‑M2/7‑150 (MQ=255)
                                    ccGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGc                                                                                                            >  1:923292/1‑150 (MQ=255)
                                          ccAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGGCCATCGATGGTCGc                                                                                                                                                <  2:1174150/108‑1 (MQ=255)
                                          ccAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGGCCATCGATGGTCGc                                                                                                                                                >  1:1174150/1‑108 (MQ=255)
                                              aaGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCTGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCAc                                                                                                  <  1:1009599/150‑1 (MQ=255)
                                                gCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACgc                                                                                                >  1:250143/1‑150 (MQ=255)
                                                   caccacAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACgccag                                                                                             >  1:70826/1‑150 (MQ=255)
                                                     ccacAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGTGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTg                                                                                                          <  1:449297/135‑1 (MQ=255)
                                                     ccacAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTg                                                                                                          >  2:449297/1‑135 (MQ=255)
                                                     ccacAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACgccagcc                                                                                           <  1:1146058/150‑1 (MQ=255)
                                                       acaGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCaa                                                                                         >  1:1003232/1‑150 (MQ=255)
                                                        caGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCaaa                                                                                        <  2:1003232/150‑1 (MQ=255)
                                                             aGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCgg                                                                                   <  2:1347663/150‑1 (MQ=255)
                                                             aGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCgg                                                                                   <  2:1006231/150‑1 (MQ=255)
                                                               cGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGc                                                                                 <  2:70826/150‑1 (MQ=255)
                                                                gAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCt                                                                                >  2:926849/1‑150 (MQ=255)
                                                                   cacTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGAc                                                                             <  2:923292/150‑1 (MQ=255)
                                                                   cacTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGAc                                                                             <  2:558791/150‑1 (MQ=255)
                                                                    acTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACg                                                                            >  2:28093/1‑150 (MQ=255)
                                                                      tgtgGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGGCCATCGATGGTCGc                                                                                                                                                <  1:158292/80‑1 (MQ=255)
                                                                      tgtgGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGGCCATCGATGGTCGc                                                                                                                                                >  2:158292/1‑80 (MQ=37)
                                                                      tgtgGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCg                                                                                    >  1:899132/1‑140 (MQ=255)
                                                                      tgtgGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCg                                                                                    <  2:899132/140‑1 (MQ=255)
                                                                      tgtgGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCa                                                                          >  2:864784/1‑150 (MQ=255)
                                                                        tgGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCaaa                                                                        <  1:864784/150‑1 (MQ=255)
                                                                          gAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCAc                                                                                                  <  2:985563/122‑1 (MQ=255)
                                                                          gAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCAc                                                                                                  >  1:985563/1‑122 (MQ=255)
                                                                          gAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATc                                                                      >  1:1201280/1‑150 (MQ=255)
                                                                            aaaCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCa                                                                    <  1:1400395/150‑1 (MQ=255)
                                                                              aCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAg                                                                                                              <  1:1366539/106‑1 (MQ=255)
                                                                              aCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAg                                                                                                              >  2:1366539/1‑106 (MQ=255)
                                                                                cTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATa                                                                >  1:800528/1‑150 (MQ=255)
                                                                                cTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATa                                                                >  2:859254/1‑150 (MQ=255)
                                                                                   gTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGcc                                                             <  1:926849/150‑1 (MQ=255)
                                                                                         gCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCggg                                                                                  >  1:449123/1‑123 (MQ=255)
                                                                                         gCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCggg                                                                                  <  2:449123/123‑1 (MQ=255)
                                                                                           gggTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGcac                                                     >  1:265296/1‑150 (MQ=255)
                                                                                            ggTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGcaca                                                    >  1:577344/1‑150 (MQ=255)
                                                                                             gTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAg                                                   <  2:577344/150‑1 (MQ=255)
                                                                                              tttACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAgg                                                  <  1:770274/150‑1 (MQ=255)
                                                                                                   ccGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACg                                             >  1:206258/1‑150 (MQ=255)
                                                                                                    cgcgAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAg                                                   >  1:898083/1‑143 (MQ=255)
                                                                                                    cgcgAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAg                                                   <  2:898083/143‑1 (MQ=255)
                                                                                                    cgcgAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGc                                            <  1:859254/150‑1 (MQ=255)
                                                                                                     gcgAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGcc                                           >  1:146279/1‑150 (MQ=255)
                                                                                                       gAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCaa                                         >  1:316816/1‑150 (MQ=255)
                                                                                                         aGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATc                                       <  1:61391/150‑1 (MQ=255)
                                                                                                         aGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATc                                       <  2:265296/150‑1 (MQ=255)
                                                                                                            ggCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCGCAGGCCACGCCAATcgcc                                    <  2:47025/150‑1 (MQ=255)
                                                                                                                 tCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAg                               <  2:316816/150‑1 (MQ=255)
                                                                                                                         cTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACTCCGCGTACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGc                       <  2:970148/150‑1 (MQ=255)
                                                                                                                           ggCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCg                     <  1:28093/150‑1 (MQ=255)
                                                                                                                                  aaGCGGCCATCGATGGTCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCa              >  1:1335925/2‑150 (MQ=255)
                                                                                                                                       aCCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATa                                                                <  2:659011/95‑1 (MQ=255)
                                                                                                                                       aCCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATa                                                                >  1:659011/1‑95 (MQ=255)
                                                                                                                                           tCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGCTGATCTCCGGct     <  1:269938/150‑2 (MQ=255)
                                                                                                                                             gATGGCCGCTTGACCACCTCGGCATACCAcgccgc                                                                                                                      <  2:785584/35‑1 (MQ=255)
                                                                                                                                             gATGGCCGCTTGACCACCTCGGCATACCAcgccgc                                                                                                                      >  1:785584/1‑35 (MQ=255)
                                                                                                                                              aTGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGt  <  1:183117/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGT  >  NC_002947/401792‑402083

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