Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,261,122 G→T 7.4% P176Q (CCG→CAG)  PP_4637 ← 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,261,1220GT 7.4% 343.2 / 10.0 108P176Q (CCG→CAG) PP_46375‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase family protein
Reads supporting (aligned to +/- strand):  ref base G (53/47);  new base T (3/5);  total (56/52)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 111 total aligned reads displayed.

GCAGCCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATGGGTTCA  >  NC_002947/5260973‑5261271
                                                                                                                                                     |                                                                                                                                                     
gcagcCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCg                                                                                                                                                       >  1:240179/1‑150 (MQ=255)
 cagcCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCgg                                                                                                                                                      >  1:1365622/1‑150 (MQ=255)
    ccccAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAg                                                                                                                                                   <  1:2359620/150‑1 (MQ=255)
      ccAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGt                                                                                                                                                 >  2:1010040/1‑150 (MQ=255)
         gggTCTTGCCCAAAACACACGGCGAGGCGTTGCGGCCGTCGACAAGTCCCAGCGACAGTACCTTGTGGGCCGCCAGGCGATCGAGAATGGTCCGGTACTGCTCGGGCGCCTGAACCAGATCCACGTGCAGGCCATCAACCGGCAGGTTgg                                                                                                                                              <  2:161375/150‑1 (MQ=255)
           gTCGTGCCCAAAACACAAGGCGAGGCGTTGAGACCGCCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCTGCCGGTTgatg                                                                                                                                            <  2:501997/150‑4 (MQ=255)
             cTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCGCCAACCTGCCGGTTgatggc                                                                                                                                          <  2:2193052/150‑6 (MQ=255)
              ttGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGcc                                                                                                                                         >  2:157525/1‑150 (MQ=255)
                  ccAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGc                                                                                                                                     >  1:2100410/1‑150 (MQ=255)
                   cAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCGCCAACCTGCCGGTTgatgg                                                                                                                                           <  2:315412/143‑5 (MQ=255)
                   cAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCGCCAACCTGCCGGTTgatgg                                                                                                                                           >  1:315412/1‑139 (MQ=255)
                   cAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGcctgcac                                                                                                                                    <  1:453397/150‑6 (MQ=255)
                   cAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGcc                                                                                                                                    >  1:578768/1‑150 (MQ=255)
                   cAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGcc                                                                                                                                    >  2:1450468/1‑150 (MQ=255)
                    aaaTCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCg                                                                                                                                   >  1:1388718/1‑150 (MQ=255)
                          caAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTg                                                                                                                             >  2:2194673/1‑150 (MQ=255)
                           aaGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGc                                                                                                                            <  1:1540897/150‑1 (MQ=255)
                           aaGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGc                                                                                                                            >  1:2170066/1‑150 (MQ=255)
                              gcgaggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCTGCCGGTTgatggc                                                                                                                                          >  1:2290360/1‑128 (MQ=255)
                              gcgaggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCTGCCGGTTgatggc                                                                                                                                          <  2:2290360/133‑6 (MQ=255)
                                gaggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTcc                                                                                                                       <  1:2194673/150‑1 (MQ=255)
                                 aggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCa                                                                                                                      <  2:1388718/150‑1 (MQ=255)
                                  ggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAg                                                                                                                     <  2:240179/150‑1 (MQ=255)
                                           aCCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGaa                                                                                                            <  1:1801103/150‑1 (MQ=255)
                                           aCCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCGTCCAGGCCACCGaa                                                                                                            <  2:740406/150‑1 (MQ=255)
                                             cGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCATATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGt                                                                                                          >  1:1708581/1‑150 (MQ=255)
                                             cGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGGTCGACGTGCAGGCAATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGt                                                                                                          >  2:1303501/1‑150 (MQ=255)
                                              gTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAgg                                                                                                                    <  1:1248934/139‑1 (MQ=255)
                                              gTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAgg                                                                                                                    >  2:1248934/1‑139 (MQ=255)
                                              gTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTa                                                                                                         >  2:182905/1‑150 (MQ=255)
                                              gTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTa                                                                                                         >  1:2241778/1‑150 (MQ=255)
                                                 gACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGt                                                                                                      >  2:2054468/1‑150 (MQ=255)
                                                   ccAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCATATCGACGTGCAGGCCATCAACCGGCAg                                                                                                                                                   <  1:1856229/103‑1 (MQ=255)
                                                   ccAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCATATCGACGTGCAGGCCATCAACCGGCAg                                                                                                                                                   >  2:1856229/1‑103 (MQ=255)
                                                   ccAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGCCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTgg                                                                                                    >  2:879202/1‑150 (MQ=255)
                                                     aGGCCCAGCGACAGTACCTTGTAGGACGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCg                                                                                                  >  1:2269970/1‑150 (MQ=255)
                                                         ccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCa                                                                                              >  1:932740/1‑150 (MQ=255)
                                                         ccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCa                                                                                              >  1:163566/1‑150 (MQ=255)
                                                         ccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCa                                                                                              >  2:861380/1‑150 (MQ=255)
                                                         ccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCa                                                                                              >  2:1132923/1‑150 (MQ=255)
                                                           aGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCAAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGc                                                                                            <  2:2301831/150‑1 (MQ=255)
                                                           aGCGACAGTACCTGGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGc                                                                                            <  2:730947/150‑1 (MQ=255)
                                                             cGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGAGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGctt                                                                                          <  1:25369/150‑1 (MQ=255)
                                                                   tACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGc                                                                                                                   <  2:738081/119‑1 (MQ=255)
                                                                   tACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGc                                                                                                                   >  1:738081/1‑119 (MQ=255)
                                                                      cTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCTGCCGGTTgatg                                                                                                                                            >  2:895200/1‑88 (MQ=255)
                                                                      cTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCTGCCGGTTgatg                                                                                                                                            <  1:895200/91‑4 (MQ=255)
                                                                      cTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGcttc                                                                                         <  1:999333/142‑1 (MQ=255)
                                                                       ttGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGt                                                                                <  2:942380/150‑1 (MQ=255)
                                                                          tAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTgcg                                                                             >  1:1109926/1‑150 (MQ=255)
                                                                           aGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTgcgc                                                                            <  2:2170066/150‑1 (MQ=255)
                                                                             gCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGt                                                                          >  2:243506/1‑150 (MQ=255)
                                                                                 gCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTgta                                                                      <  1:1098935/150‑3 (MQ=255)
                                                                                   aGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAgg                                                                    >  1:2229403/1‑150 (MQ=255)
                                                                                   aGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAgg                                                                    <  2:2098297/150‑1 (MQ=255)
                                                                                     gCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGCGCGCGTTGAAGGAt                                                                  <  2:393841/150‑1 (MQ=255)
                                                                                       gATCGAGAATGGTCTGGTACTGCGCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGAGTCGATCAGCTTCTTCAACGGTGCGCGTTGAAggttct                                                                >  1:500505/1‑146 (MQ=255)
                                                                                        aTCGAGAGTGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGAtgtt                                                               <  1:1823300/150‑1 (MQ=255)
                                                                                         tCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGtt                                                                                                                              >  1:1537368/1‑86 (MQ=255)
                                                                                         tCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGAtgttg                                                              <  2:445839/150‑1 (MQ=255)
                                                                                         tCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGAtgttg                                                              <  1:243506/150‑1 (MQ=255)
                                                                                          cGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGata                                                                                                                                >  2:1913504/1‑81 (MQ=255)
                                                                                          cGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGata                                                                                                                                <  1:1913504/83‑3 (MQ=255)
                                                                                          cGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGt                                                                                                          >  1:710288/1‑105 (MQ=255)
                                                                                                     tGGTACTGCTCGGGCGCCTGAACCATATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCa                                                  <  2:1708581/150‑1 (MQ=255)
                                                                                                       gTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATa                                                <  1:1450468/150‑1 (MQ=255)
                                                                                                        tACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAg                                               <  1:1010040/150‑1 (MQ=255)
                                                                                                         aCTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCa                                                                                                                      >  2:1174146/1‑78 (MQ=255)
                                                                                                         aCTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCa                                                                                                                      <  1:1174146/78‑1 (MQ=255)
                                                                                                          cTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTggc                                                                                                                                             >  1:701706/1‑54 (MQ=255)
                                                                                                          cTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTggc                                                                                                                                             <  2:701706/54‑1 (MQ=255)
                                                                                                          cTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGCGGCTATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGc                                             <  2:2432681/150‑1 (MQ=255)
                                                                                                           tGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTGGTATCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCg                                            <  1:157525/150‑1 (MQ=255)
                                                                                                            gCTCCGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGt                                           <  2:949471/150‑1 (MQ=255)
                                                                                                                 ggCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTc                                      <  2:2229403/150‑1 (MQ=255)
                                                                                                                   cgcCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCa                                    <  1:879202/150‑1 (MQ=255)
                                                                                                                         aaCCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGtt                                                                                                                              >  2:332820/1‑54 (MQ=255)
                                                                                                                         aaCCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGtt                                                                                                                              <  1:332820/54‑1 (MQ=255)
                                                                                                                          aCCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGACGGATGTTGTAGACGCGTTCATATGCGTTTTTTCATTCCTgg                             <  2:983301/150‑1 (MQ=255)
                                                                                                                           ccAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTggg                            >  1:1030251/1‑150 (MQ=255)
                                                                                                                           ccAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTGTGTCCAGTCCTggg                            >  1:1696622/1‑150 (MQ=255)
                                                                                                                             aGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAg                                                                                             >  1:2059982/1‑83 (MQ=255)
                                                                                                                              gATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCa                         >  1:1582415/1‑150 (MQ=255)
                                                                                                                               aTCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAg                        >  2:1551503/1‑150 (MQ=255)
                                                                                                                                    cGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCt                   <  1:1551503/150‑1 (MQ=255)
                                                                                                                                         aGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCg              >  1:1809410/1‑150 (MQ=255)
                                                                                                                                           gCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAg            >  1:1816101/1‑150 (MQ=255)
                                                                                                                                             cATCAACCGGCAGGTTGTCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGa          >  1:2136078/1‑150 (MQ=255)
                                                                                                                                             cATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGCCGCGTTAATAGGGGATCTGCCTGCGCTGGGGCAGCGCCTGGGCGAgag          >  2:1453869/1‑148 (MQ=255)
                                                                                                                                             cATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGa          >  1:380078/1‑150 (MQ=255)
                                                                                                                                             cATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGa          >  1:2032709/1‑150 (MQ=255)
                                                                                                                                              aTCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGAt         >  1:904813/1‑150 (MQ=255)
                                                                                                                                               tCAACCGGCAGGTTGTCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATg        <  2:2136078/150‑1 (MQ=255)
                                                                                                                                                cAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATgg       <  1:792139/150‑1 (MQ=255)
                                                                                                                                                cAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATgg       <  2:2211461/150‑1 (MQ=255)
                                                                                                                                                 aaCCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATggg      <  1:70503/150‑1 (MQ=255)
                                                                                                                                                    cGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGtt                                                                         >  1:1894664/1‑80 (MQ=255)
                                                                                                                                                    cGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGtt                                                                         <  2:1894664/80‑1 (MQ=255)
                                                                                                                                                    cGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATGGGTTc   >  1:1703984/1‑150 (MQ=255)
                                                                                                                                                     ggCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATGGGTTCa  >  1:1122922/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAGCCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATGGGTTCA  >  NC_002947/5260973‑5261271

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