Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 401,927 A→G 6.7% 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,9270AG 6.7% 223.1 / 4.6 74V201A (GTC→GCC) PP_0335glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base A (31/38);  new base G (3/2);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

ACACATAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCA  >  NC_002947/401778‑402071
                                                                                                                                                     |                                                                                                                                                
ACACATAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCCCTGGCGCTGTAGCGa                                                                                                                                                  <  2:269938‑M1/150‑1 (MQ=255)
  acaTAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGAcc                                                                                                                                                >  1:558791/1‑150 (MQ=255)
  acaTAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGAcc                                                                                                                                                >  2:1400395/1‑150 (MQ=255)
  acaTAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGACGCTGTAGCGAcc                                                                                                                                                >  2:740016/1‑150 (MQ=255)
     tAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATc                                                                                                                                             >  2:1096651/1‑150 (MQ=255)
       aTGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGa                                                                                                                                           >  2:1146058/1‑150 (MQ=255)
        tGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGAt                                                                                                                                          >  1:47025/1‑150 (MQ=255)
        tGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGAt                                                                                                                                          <  2:50767/150‑1 (MQ=255)
              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                                                                                                         >  2:183117‑M2/7‑150 (MQ=255)
                                         cgggcaCGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACgccgcc                                                                                                         <  1:1196842‑M2/144‑1 (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                                                                                                  >  2:1366539/1‑106 (MQ=255)
                                                                                            aCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAg                                                                                                  <  1:1366539/106‑1 (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                                       <  2:898083/143‑1 (MQ=255)
                                                                                                                  cgcgAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAg                                       >  1:898083/1‑143 (MQ=255)
                                                                                                                  cgcgAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGc                                <  1:859254/150‑1 (MQ=255)
                                                                                                                   gcgAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGcc                               >  1:146279/1‑150 (MQ=255)
                                                                                                                     gAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCaa                             >  1:316816/1‑150 (MQ=255)
                                                                                                                       aGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATc                           <  2:265296/150‑1 (MQ=255)
                                                                                                                       aGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATc                           <  1:61391/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)
                                                                                                                                                     |                                                                                                                                                
ACACATAATGGGGTGTCAGGCACCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCA  >  NC_002947/401778‑402071

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: