Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,261,122 G→T 7.0% 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.0% 363.2 / 11.5 115P176Q (CCG→CAG) PP_46375‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase family protein
Reads supporting (aligned to +/- strand):  ref base G (54/53);  new base T (2/6);  total (56/59)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

Only 100 of 117 total aligned reads displayed.

GCAGCCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATGGGTTCA  >  NC_002947/5260973‑5261271
                                                                                                                                                     |                                                                                                                                                     
      ccAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGt                                                                                                                                                 <  1:775196/150‑1 (MQ=255)
      ccAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGt                                                                                                                                                 >  2:2045355/1‑150 (MQ=255)
       cAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGtt                                                                                                                                                >  2:1205539/1‑150 (MQ=255)
        aGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTg                                                                                                                                               <  1:2316282/150‑1 (MQ=255)
          ggTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCACCCGGCAGGTTtgc                                                                                                                                             <  2:33001/150‑1 (MQ=255)
          ggTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCTGCCGGTTgat                                                                                                                                             <  2:2110824/150‑3 (MQ=255)
           gTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCTGCCGGTTgatg                                                                                                                                            <  2:1627615/150‑4 (MQ=255)
             cTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCTGCCGGTTgatggc                                                                                                                                          <  2:2365080/150‑6 (MQ=255)
             cTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCTGCCGGTTgatggc                                                                                                                                          <  1:1222328/150‑6 (MQ=255)
             cTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACAGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTggcggc                                                                                                                                          <  2:260282/150‑1 (MQ=255)
              ttGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGcc                                                                                                                                         >  1:1916418/1‑150 (MQ=255)
              ttGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGcc                                                                                                                                         >  1:2227394/1‑150 (MQ=255)
              ttGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGcc                                                                                                                                         >  1:1453972/1‑150 (MQ=255)
                  ccAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGc                                                                                                                                     <  2:561095/150‑1 (MQ=255)
                  ccAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGc                                                                                                                                     >  1:1712766/1‑150 (MQ=255)
                  ccAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCGGAACCGGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGc                                                                                                                                     >  1:1468481/1‑150 (MQ=255)
                   cAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTACTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAg                                                                                                                                                   <  1:1908013/135‑1 (MQ=255)
                   cAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTACTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAg                                                                                                                                                   >  2:1908013/1‑135 (MQ=255)
                     aaTCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGa                                                                                                                                  <  2:1010557/150‑1 (MQ=255)
                       gcacaAGGCGCGGCGTGGCGACCGGCGACCAGGCGCAGCGACAGTACTTTGTAGGGCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGGGTCGGAACCCGATCGAGGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAgg                                                                                                                                <  2:696159/149‑1 (MQ=255)
                         acaAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCTACTTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGtt                                                                                                                              >  2:1475592/1‑150 (MQ=255)
                         acaAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCGCCAACCTGCCGGTTgatggc                                                                                                                                          >  1:2358867/1‑133 (MQ=255)
                         acaAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGtt                                                                                                                              >  1:167873/1‑150 (MQ=255)
                         acaAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGtt                                                                                                                              >  1:1752070/1‑150 (MQ=255)
                         acaAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGtt                                                                                                                              <  2:290739/150‑1 (MQ=255)
                         acaAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGtt                                                                                                                              >  1:2050547/1‑150 (MQ=255)
                          caAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTg                                                                                                                             >  1:1342990/1‑150 (MQ=255)
                                 aggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCa                                                                                                                      <  2:2171009/150‑1 (MQ=255)
                                  ggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGTCGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAg                                                                                                                     <  1:1985885/150‑1 (MQ=255)
                                  ggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAg                                                                                                                     <  2:1041953/150‑1 (MQ=255)
                                  ggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAg                                                                                                                     <  1:535149/150‑1 (MQ=255)
                                     gTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGGGCCTGAACCAGATTGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCTAGTcc                                                                                                                  >  1:18428/1‑150 (MQ=255)
                                     gTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGcc                                                                                                                  >  2:1046665/1‑150 (MQ=255)
                                       tGCGACCGTCGACCAGGCCTAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTTGCGGCCAGGCCTATGTTGCCTTCCAGGCCAc                                                                                                                <  2:500533/150‑1 (MQ=255)
                                           aCCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCCGGTTGGCGGCCAGGCAGAGGTTGCCTTCCAGGCCCCCGaa                                                                                                            <  2:343239/150‑1 (MQ=255)
                                            ccGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAg                                                                                                           >  1:1403016/1‑150 (MQ=255)
                                             cGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGt                                                                                                          <  2:1712766/150‑1 (MQ=255)
                                               tCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAgaag                                                                                                        >  1:610726/1‑147 (MQ=255)
                                                   ccAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTgg                                                                                                    <  1:871629/150‑1 (MQ=255)
                                                      ggCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGaa                                                                                                            >  2:515511/1‑139 (MQ=255)
                                                      ggCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGaa                                                                                                            <  1:515511/139‑1 (MQ=255)
                                                         ccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGc                                                                                                                                                     >  2:1638746/1‑95 (MQ=255)
                                                         ccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGc                                                                                                                                                     <  1:1638746/95‑1 (MQ=255)
                                                         ccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCa                                                                                              >  1:1095988/1‑150 (MQ=255)
                                                         ccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCa                                                                                              <  1:1746001/150‑1 (MQ=255)
                                                         ccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCa                                                                                              >  2:321444/1‑150 (MQ=255)
                                                              gACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCTGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGcttc                                                                                         <  2:90366/150‑1 (MQ=255)
                                                              gACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGAGTAGcttc                                                                                         >  1:1307852/1‑150 (MQ=255)
                                                              gACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCATGCCAGGCCACCGAAGTACGTGGCGATCAGcttc                                                                                         <  1:2371447/150‑1 (MQ=255)
                                                              gACAGTACCTTGTAGGCCGGCAGGAGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCACCCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGcttc                                                                                         <  2:412182/150‑1 (MQ=255)
                                                               aCAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGcttct                                                                                        >  1:1103328/1‑150 (MQ=255)
                                                                  gTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCCCCGAAGTACGTGGCGATCAGCTTCTTCa                                                                                     >  2:925714/1‑150 (MQ=255)
                                                                     ccTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACg                                                                                  <  2:2227394/150‑1 (MQ=255)
                                                                      cTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCTGCCGGTTgatg                                                                                                                                            >  2:1611536/1‑88 (MQ=255)
                                                                         gTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGAGCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCATGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTgc                                                                              <  2:112759/150‑1 (MQ=255)
                                                                           aGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTgcgc                                                                            >  1:2065869/1‑150 (MQ=255)
                                                                               cGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTg                                                                        >  1:1043209/1‑150 (MQ=255)
                                                                               cGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTg                                                                        <  2:189822/150‑1 (MQ=255)
                                                                                   aGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAgg                                                                    <  2:1403016/150‑1 (MQ=255)
                                                                                    ggCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGtt                                                                         <  2:2173237/144‑1 (MQ=255)
                                                                                    ggCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGtt                                                                         >  1:2173237/1‑144 (MQ=255)
                                                                                     gCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGAt                                                                  <  1:1581005/150‑1 (MQ=255)
                                                                                       gATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTCCGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGAtgt                                                                >  1:633701/1‑150 (MQ=255)
                                                                                       gATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGAtgt                                                                <  1:1071722/150‑1 (MQ=255)
                                                                                         tCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGcttct                                                                                        <  1:2147805/124‑1 (MQ=255)
                                                                                         tCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGAtgttg                                                              <  2:665671/150‑1 (MQ=255)
                                                                                         tCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGCCCAgg                                                                                                                                      >  1:1227606/1‑78 (MQ=255)
                                                                                         tCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGCCCAgg                                                                                                                                      <  2:1227606/78‑1 (MQ=255)
                                                                                          cGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGAtgttgt                                                             <  1:1156443/150‑1 (MQ=255)
                                                                                            agaATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAGTTACGTTGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAg                                                           >  2:1871941/1‑150 (MQ=255)
                                                                                             gaATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGaa                                                                      <  1:1107636/138‑1 (MQ=255)
                                                                                             gaATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGaa                                                                      >  2:1107636/1‑138 (MQ=255)
                                                                                                  gTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGt                                                     >  1:1421264/1‑150 (MQ=255)
                                                                                                     tGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCa                                                  >  2:57554/1‑150 (MQ=255)
                                                                                                     tGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTAGAGCTCCGATTATTCAAAGGCGCGCGTGGAAGGAGGTTGTAGACGCGTTCa                                                  >  2:2070067/1‑150 (MQ=255)
                                                                                                      ggTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCAt                                                 <  2:1752070/150‑1 (MQ=255)
                                                                                                        tACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAg                                               <  2:1060437/150‑1 (MQ=255)
                                                                                                           tGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGtt                                                                                                                              <  2:874131/68‑1 (MQ=255)
                                                                                                           tGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGtt                                                                                                                              >  1:874131/1‑68 (MQ=255)
                                                                                                           tGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCg                                            >  2:2240250/1‑150 (MQ=255)
                                                                                                            gCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGt                                           >  2:991990/1‑150 (MQ=255)
                                                                                                                   cgcCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCa                                    <  2:2065869/150‑1 (MQ=255)
                                                                                                                    gcCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAg                                   >  2:1053037/1‑150 (MQ=255)
                                                                                                                        gAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCt                               >  1:1323098/1‑150 (MQ=255)
                                                                                                                          aCCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTgg                             <  2:2216793/150‑1 (MQ=255)
                                                                                                                           ccAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCCGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATTGGCGTTTTTCCAGTCCTggg                            >  1:896421/1‑150 (MQ=255)
                                                                                                                                tCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAgc                       <  1:49143/150‑1 (MQ=255)
                                                                                                                                     gTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGaa                                                                                                            >  1:1569315/1‑60 (MQ=255)
                                                                                                                                     gTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGaa                                                                                                            <  2:1569315/60‑1 (MQ=255)
                                                                                                                                       gCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTggg                <  1:2070067/150‑1 (MQ=255)
                                                                                                                                        cAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTCGAAGGATGTTGTAGCCGCGTTCATAGGCGTTTTTCCAGTCCTGGGCCAGCGCCTgcgc               >  2:2259166/1‑150 (MQ=255)
                                                                                                                                         aGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCg              >  1:1563200/1‑150 (MQ=255)
                                                                                                                                          ggCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGa             <  2:1043209/150‑1 (MQ=255)
                                                                                                                                              aTCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGAt         >  1:2219891/1‑150 (MQ=255)
                                                                                                                                              aTCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGAt         <  2:2100821/150‑1 (MQ=255)
                                                                                                                                               tCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATg        >  1:2159689/1‑150 (MQ=255)
                                                                                                                                                 aaCCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATggg      >  2:1522926/1‑150 (MQ=255)
                                                                                                                                                   ccGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATGGGtt    <  1:2240250/150‑1 (MQ=255)
                                                                                                                                                    cGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGGAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATGGGTTc   >  2:1747330/1‑150 (MQ=255)
                                                                                                                                                     ggCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATGGGTTCa  >  2:2038448/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCAGCCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTCAACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATGGGTTCA  >  NC_002947/5260973‑5261271

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