Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,050,224 G→T 6.1% T109T (ACC→ACA)  PP_1824 ← Smr domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,050,2240GT 6.1% 202.8 / 2.7 66T109T (ACC→ACA) PP_1824Smr domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base T (2/2);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

GAGCACTTGCGGGTGCTGGCGTAGCCAGGTGTTGACGTGGCTTTTGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCCACATCAGCTCGTCTTCGGCGCCTACGTCGATGACGAA  >  NC_002947/2050080‑2050373
                                                                                                                                                |                                                                                                                                                     
gAGCACTTGCGGGTGCTGGCGTAGCCAGGTGTTGACGTGGCTTTTGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTc                                                                                                                                                  >  1:1080199/1‑150 (MQ=255)
 aGCACTTGCGGGTGCTGGCGTAGCCAGGTGTTGACGGGGCTTTTGATCATGGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTcg                                                                                                                                                 >  1:669592/1‑150 (MQ=255)
    aCTTGCGGGTGCTGGCGTAGCCAGGTGTTGACGTGGCTTTTGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCgg                                                                                                                                              <  2:225805/150‑1 (MQ=255)
     cTTGCGGGTGCTGGCGTAGCCAGGTGTTGACGTGGCTTTTGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCggg                                                                                                                                             >  2:795354/1‑150 (MQ=255)
     cTTGCGGGTGCTGGCGTAGCCAGGTGTTGACGTGGCTTTTGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCAAAAGGTTTCGCggg                                                                                                                                             >  1:75780/1‑150 (MQ=255)
      ttGCGGGTGCTGGCGTAGCCAGGTGTTGACGTGGCTTTTGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGc                                                                                                                                            >  1:825553/1‑150 (MQ=255)
                tGGCGTAGCCAGGTGTTGACGTGGCTTTTGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCAcc                                                                                                                                  <  2:75780/150‑1 (MQ=255)
                tGGCGTAGCCAGGTGTTGACGTGGCTTTTGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCAcc                                                                                                                                  <  1:452344/150‑1 (MQ=255)
                tGGCGTAGCCAGGTGTTGACGTGGCTTTTGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCAcc                                                                                                                                  >  1:506379/1‑150 (MQ=255)
                    gTAGCCAGGTGTTGACGTGGCTTTTGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCa                                                                                                                              >  1:1164698/1‑150 (MQ=255)
                                   cGTGGCTTTTGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGa                                                                                                               <  2:564907/150‑1 (MQ=255)
                                     tGGCTTTTGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAgg                                                                                                             >  2:402036/1‑150 (MQ=255)
                                     tGGCTTTTGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAgg                                                                                                             >  2:36881/1‑150 (MQ=255)
                                     gggCTTTTGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAgg                                                                                                             <  2:24398/149‑1 (MQ=255)
                                         ttttGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAg                                                                                                                    <  2:71651/139‑1 (MQ=255)
                                         ttttGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAg                                                                                                                    >  1:71651/1‑139 (MQ=255)
                                               tcatcGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGa                                                                                                   >  1:1285111/1‑150 (MQ=255)
                                                catcGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGaa                                                                                                  >  1:836680/1‑150 (MQ=255)
                                                 atcGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAAc                                                                                                 >  2:1119703/1‑150 (MQ=255)
                                                 atcGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAAc                                                                                                 <  2:1164698/150‑1 (MQ=255)
                                                 atcGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAAc                                                                                                 <  2:825553/150‑1 (MQ=255)
                                                       cgcTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGATCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATc                                                                                           <  1:277026/150‑1 (MQ=255)
                                                          ttGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGc                                                                                        <  1:1119703/150‑1 (MQ=255)
                                                            gCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGccc                                                                                      <  1:365256/150‑1 (MQ=255)
                                                              cGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCca                                                                                    >  2:845995/1‑150 (MQ=255)
                                                               gTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCcag                                                                                   >  1:585980/1‑150 (MQ=255)
                                                                tCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCcagc                                                                                  >  1:929878/1‑150 (MQ=255)
                                                                  cAGGCGCGCGGCCGTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCcagctt                                                                                >  2:894549/1‑150 (MQ=255)
                                                                    ggCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCcagcttca                                                                              >  2:473727/1‑150 (MQ=255)
                                                                    ggCGCGCGGCCTTGCCGTACTGCACCCGCACACAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAACGTCTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCcagcttca                                                                              <  1:707272/150‑1 (MQ=255)
                                                                       gcgcgGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTACACCGTCATGCCGTGCAGGTCGAGGCTTCCCTCGAACGGAATTTGCCCcagcttcagct                                                                           >  1:589870/1‑150 (MQ=255)
                                                                          cgGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTgc                                                                        <  2:506379/150‑1 (MQ=255)
                                                                           ggCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTgcg                                                                       >  1:582039/1‑150 (MQ=255)
                                                                           ggCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTgcg                                                                       >  1:220835/1‑150 (MQ=255)
                                                                            gCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTgcgc                                                                      <  2:582039/150‑1 (MQ=255)
                                                                            gCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTgcgc                                                                      >  1:814118/1‑150 (MQ=255)
                                                                             ccTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCa                                                                     >  2:725953/1‑150 (MQ=255)
                                                                             ccTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCa                                                                     >  2:206074/1‑150 (MQ=255)
                                                                               ttGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTCTTGTGACTGCGCTGATGAAGTCCCACAAGGTTTGGCTGTCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGc                                                                   >  1:543167/1‑150 (MQ=255)
                                                                                   cGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGc                                                               <  1:845995/150‑1 (MQ=255)
                                                                                     tGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCcagc                                                                                  <  2:1027333/129‑1 (MQ=255)
                                                                                     tGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCcagc                                                                                  >  1:1027333/1‑129 (MQ=255)
                                                                                           aCCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGc                                                       >  2:822744/1‑150 (MQ=255)
                                                                                             ccGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAATGTTTCGCGGGCCTTTTCCACCGTCATGcc                                                                                                                          >  1:51493/1‑81 (MQ=255)
                                                                                             ccGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAATGTTTCGCGGGCCTTTTCCACCGTCATGcc                                                                                                                          >  1:51481/1‑81 (MQ=255)
                                                                                             ccGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAATGTTTCGCGGGCCTTTTCCACCGTCATGcc                                                                                                                          <  2:51481/81‑1 (MQ=255)
                                                                                             ccGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAATGTTTCGCGGGCCTTTTCCACCGTCATGcc                                                                                                                          <  2:51493/81‑1 (MQ=255)
                                                                                               gcacgcaGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTGGAGGCTGCCATCAACCGGAGTCTGACCGCGCGTTGGCTGGGGCAGCGGGCCTGCCGGCACGc                                                   >  2:140494/1‑150 (MQ=255)
                                                                                                        gcACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCgcggc                                          <  1:402036/150‑1 (MQ=255)
                                                                                                             tCCAGTTTTGTGGCTTCGGCGATGAAGTCCCCCAGGGGTTCGCTGACCATTGCCGCAGTAATGCCGTGCACGTCGACCACTACCACAAGAAAAATGTTCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTATCCGCCATCGCGGCGCcac                                     >  1:1255002/1‑150 (MQ=255)
                                                                                                             tCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTg                                                                 <  1:57591/122‑1 (MQ=255)
                                                                                                             tCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTg                                                                 >  2:57591/1‑122 (MQ=255)
                                                                                                                   tttGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGACCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGTACGCCATCGCGGCGCCACATCAGc                               <  2:1255002/150‑1 (MQ=255)
                                                                                                                     tgtgGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCCACATCAGCTc                             <  1:36881/150‑1 (MQ=255)
                                                                                                                        ggCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCCACATCAGCTCGTc                          <  2:585980/150‑1 (MQ=255)
                                                                                                                        ggCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCCACATCAGCTCGTc                          <  2:929878/150‑1 (MQ=255)
                                                                                                                          cTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCCACATCAGCTCGTCtt                        <  1:513745/150‑1 (MQ=255)
                                                                                                                           ttCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCCACATCAGCTCGTCTTc                       <  1:795354/150‑1 (MQ=255)
                                                                                                                                cGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCCACATCAGCTCGTCTTCGgcgc                  <  1:725953/150‑1 (MQ=255)
                                                                                                                                 gATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCCACATCAGCTCGTCTTCGGCGcc                 >  1:753171/1‑150 (MQ=255)
                                                                                                                                          ccACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCca                                      <  1:150754/120‑1 (MQ=255)
                                                                                                                                          ccACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCca                                      >  2:150754/1‑120 (MQ=255)
                                                                                                                                          ccACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCCACATCAGCTCGTCTTCGGCGCCTACGTCGAt        <  2:669592/150‑1 (MQ=255)
                                                                                                                                          ccACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCCACATCAGCTCGTCTTCGGCGCCTACGTCGAt        >  1:3657/1‑150 (MQ=255)
                                                                                                                                          ccACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCCACATCAGCTCGTCTTCGGCGCCTACGTCGAt        <  2:1080199/150‑1 (MQ=255)
                                                                                                                                          ccACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCCACATCAGCTCGTCTTCGGCGCCTACGTCGAt        <  1:473727/150‑1 (MQ=255)
                                                                                                                                                ggTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCCACATCAGCTCGTCTTCGGCGCCTACGTCGATGACGaa  <  2:3657/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GAGCACTTGCGGGTGCTGGCGTAGCCAGGTGTTGACGTGGCTTTTGATCATCGGGCGCTTGCCGTCCAGGCGCGCGGCCTTGCCGTGGGTCACCCGCACGCAGCGCACCTCCAGTTTTGTGGCTTCGGCGATGAAGTCCCACAAGGTTTCGCGGGCCTTTTCCACCGTCATGCCGTGCAGGTCGAGGCTGCCCTCGAACGGAATCTGCCCCAGCTTCAGCTTGCGCAGCTGGCCTTCCTGCACGCCATCGCGGCGCCACATCAGCTCGTCTTCGGCGCCTACGTCGATGACGAA  >  NC_002947/2050080‑2050373

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