Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,348,162 T→G 5.6% E173A (GAG→GCG)  PP_3820 ← group II intron‑encoding maturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,348,1620TG 5.6% 218.3 / 2.9 81E173A (GAG→GCG) PP_3820group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base T (40/35);  new base G (3/2);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAATAAGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCAGTGCTTGCTGG  >  NC_002947/4348017‑4348263
                                                                                                                                                 |                                                                                                      
gcgAATAAGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGAt                                                                                                    <  1:1086182/150‑1 (MQ=32)
gcgAATAAGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGACCGAGCCAGCTCGAt                                                                                                    >  2:1877866/1‑150 (MQ=12)
   aaTAAGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGc                                                                                                 <  1:2371953/150‑1 (MQ=32)
    aTAAGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCt                                                                                                <  2:1396848/150‑1 (MQ=32)
    aTAAGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATACCCGTATTTCTCCAGATTGTGGGCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCt                                                                                                <  1:1485600/150‑1 (MQ=11)
     tAAGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCtt                                                                                               >  2:977455/1‑150 (MQ=32)
        gCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGAt                                                                                            >  2:2057764/1‑150 (MQ=32)
          cTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGc                                                                                          <  1:1099587/150‑1 (MQ=32)
           tGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCg                                                                                         <  2:1556893/150‑1 (MQ=32)
             aGTACGCACTTGTCGTTGACCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTa                                                                                       <  1:301650/150‑1 (MQ=17)
                    aCTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGCCCCGATCGAAAAACTTCTCAAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGc                                                                                <  1:1303537/150‑1 (MQ=17)
                    aCTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGc                                                                                >  2:1445140/1‑150 (MQ=32)
                     cTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGa                                                                                                     >  2:1791509/1‑128 (MQ=33)
                     cTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGa                                                                                                     <  1:1791509/128‑1 (MQ=33)
                           tttGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGAt                                                                                                    <  1:1039720/123‑1 (MQ=33)
                           tttGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGAt                                                                                                    >  2:1039720/1‑123 (MQ=33)
                           tttGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACg                                                                         >  2:1190709/1‑150 (MQ=32)
                           tttGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCGAGATGGCtcgtg                                                                                            >  1:722591/1‑127 (MQ=12)
                           tttGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCGAGATGGCtcgtg                                                                                            <  2:722591/131‑5 (MQ=255)
                           tttGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCACGCTGTCACTTGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACg                                                                         >  2:1638490/1‑150 (MQ=17)
                               aTCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGCCACATGAGCACGAGCCATCTCGATg                                                                                                   <  1:1800963/120‑1 (MQ=18)
                               aTCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGCCACATGAGCACGAGCCATCTCGATg                                                                                                   >  2:1800963/1‑120 (MQ=18)
                                   gACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCAt                                                                 >  1:1141798/1‑150 (MQ=32)
                                   gACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCAt                                                                 >  1:1557195/1‑150 (MQ=32)
                                   gACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCAt                                                                 >  2:751418/1‑150 (MQ=32)
                                    aCGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATa                                                                >  2:1670778/1‑150 (MQ=32)
                                          gATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCCGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTa                                                          >  2:2124869/1‑150 (MQ=17)
                                              cAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTcc                                                     <  2:1375134/150‑1 (MQ=32)
                                               aGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCg                                                    >  1:366977/1‑150 (MQ=32)
                                               aGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCg                                                    >  1:367339/1‑150 (MQ=32)
                                               aGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACAAGCCATCGAGATGGCtcgtgc                                                                                           <  2:616739/112‑6 (MQ=11)
                                               aGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACAAGCCATCGAGATGGCtcgtgc                                                                                           >  1:616739/1‑107 (MQ=11)
                                                ggCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCga                                                   <  1:1445140/150‑1 (MQ=32)
                                                  ccANGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCgaga                                                 >  2:1511722/1‑150 (MQ=17)
                                                   cATNAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGaa                                                >  2:1516898/1‑150 (MQ=17)
                                                   cATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGaa                                                >  1:1229474/1‑150 (MQ=32)
                                                    aTGAGGATGTCGTGATTGACCCGATCGAAAAACTACTGCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAAc                                               >  1:787466/1‑150 (MQ=17)
                                                        ggATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTg                                           <  1:452717/150‑1 (MQ=32)
                                                         gATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTgg                                          <  2:94146/150‑1 (MQ=32)
                                                          aTGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTggg                                         <  2:1557195/150‑1 (MQ=32)
                                                           tGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCt                                                             <  2:1474398/130‑1 (MQ=33)
                                                           tGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCt                                                             >  1:1474398/1‑130 (MQ=33)
                                                             tCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATg                                                                                                   >  2:1102737/1‑90 (MQ=35)
                                                             tCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATg                                                                                                   <  1:1102737/90‑1 (MQ=35)
                                                              cGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACCACTTCCACCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGa                                     >  1:1001797/1‑150 (MQ=11)
                                                               gTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGaa                                    <  1:977455/150‑1 (MQ=32)
                                                                 gATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGa                                  <  2:1229474/150‑1 (MQ=32)
                                                                  aTTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGAt                                 >  1:251296/1‑150 (MQ=32)
                                                                   ttGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATa                                >  1:7803/1‑150 (MQ=32)
                                                                       cccGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGgtg                            >  2:1289997/1‑150 (MQ=32)
                                                                       cccGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGgtg                            >  2:30912/1‑150 (MQ=32)
                                                                        ccGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGgtgt                           <  1:1670778/150‑1 (MQ=32)
                                                                        ccGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGgtgt                           >  2:194108/1‑150 (MQ=32)
                                                                         cGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGtgtg                          <  1:1511722/150‑1 (MQ=32)
                                                                           aTCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAg                        >  1:1670449/1‑150 (MQ=32)
                                                                              gAAAAACTTCTCAAGATCGAGTTGCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACAAGCAATCGCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGAGAGGTGGGAgcgg                     >  2:1149341/1‑148 (MQ=11)
                                                                               aaaaaCTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCga                                                   >  2:830715/1‑119 (MQ=33)
                                                                               aaaaaCTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCga                                                   <  1:830715/119‑1 (MQ=33)
                                                                               aaaaaCTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAgctgc                    <  1:2011860/150‑1 (MQ=32)
                                                                               aaaaaCTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAgctgc                    <  1:751418/150‑1 (MQ=32)
                                                                               aaaaaCTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGAGGGGTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAgctgc                    >  1:1993253/1‑150 (MQ=17)
                                                                                   aCTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGtgtgg                         >  1:1922326/1‑140 (MQ=17)
                                                                                   aCTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGa                         <  2:1922326/141‑1 (MQ=32)
                                                                                      tctcCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCa             <  2:366977/150‑1 (MQ=32)
                                                                                      tctcCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCa             <  2:367339/150‑1 (MQ=32)
                                                                                        tcCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCAGt           <  1:30912/150‑1 (MQ=32)
                                                                                         ccAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTAc                                                                                      >  1:541874/1‑75 (MQ=34)
                                                                                         ccAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTAc                                                                                      <  2:541874/75‑1 (MQ=35)
                                                                                         ccAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCAGtg          >  2:180503/1‑150 (MQ=32)
                                                                                         ccAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCAGtg          >  2:1992888/1‑150 (MQ=32)
                                                                                          cAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCAGtgc         >  2:987662/1‑150 (MQ=32)
                                                                                           aGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCAGtgct        <  1:1190709/150‑1 (MQ=32)
                                                                                            gATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGt                                                        <  2:708053/102‑1 (MQ=32)
                                                                                            gATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGt                                                        >  1:708053/1‑102 (MQ=32)
                                                                                            gATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCAGtgctt       >  2:1736494/1‑150 (MQ=32)
                                                                                             aTCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGt                                            <  2:72185/112‑1 (MQ=32)
                                                                                             aTCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGt                                            >  1:72185/1‑112 (MQ=33)
                                                                                             aTCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTCAGTGACATTAACAGCGGGTCGAAGAGAGGTGTGAGCTGCTGCTGCAGtgcttg      >  2:63238/1‑150 (MQ=11)
                                                                                              tCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCAGtgcttgc     >  2:2281144/1‑150 (MQ=32)
                                                                                                gAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCAGTGCTTGCTg   <  1:987662/150‑1 (MQ=32)
                                                                                                 aGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCAGTGCTTGCTgg  <  1:2281144/150‑1 (MQ=31)
                                                                                                                                                 |                                                                                                      
GCGAATAAGCCTGAGTACGCACTTGTCTTTGATCCGACGTTCGATGCAGGCCATGAGGATGTCGTGATTGACCCGATCGAAAAACTTCTCCAGATCGAGTTCCACGCACCAGCGGTGTCCCGCTGTCACATGAGCACGAGCCATCTCGATGGCTTGATGCGTACTTCTGCCCGGACGAAAACCATAGCTGTAGTCCGAGAACAGTGGGTCGAAGATAGGTGTGAGCTGCTGCTGCAGTGCTTGCTGG  >  NC_002947/4348017‑4348263

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