Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,196,328 C→A 7.7% intergenic (+126/+243) PP_5618 → / ← PP_5619 Cro/CI family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,196,3280CA 7.7% 230.5 / 3.9 74intergenic (+126/+243)PP_5618/PP_5619Cro/CI family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/34);  new base A (2/4);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGAA  >  NC_002947/4196179‑4196475
                                                                                                                                                     |                                                                                                                                                   
tCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGc                                                                                                                                                     >  1:1528737/1‑150 (MQ=255)
      gTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGa                                                                                                                                               <  1:316586/150‑1 (MQ=255)
          gcAGGGCCTGGGGATGGGAAGCCACTGCTCGCTTCTGAGCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGgt                                                                                                                                           <  2:600612/149‑1 (MQ=255)
                ccTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGAc                                                                                                                                     <  2:1712829/150‑1 (MQ=255)
                    gTGACGGTAATCCACTGCTCGCTTCTGCTCTACCCAAGGCGCGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGtct                                                                                                                                 <  1:1556297/150‑1 (MQ=255)
                    gTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGtct                                                                                                                                 >  1:1925188/1‑150 (MQ=255)
                    gTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGtct                                                                                                                                 >  1:950380/1‑150 (MQ=255)
                    gTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGtct                                                                                                                                 <  2:1528737/150‑1 (MQ=255)
                            aaGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGACGTCGACACCGTCGACGGCTCggttcac                                                                                                                         <  2:1435885/150‑6 (MQ=255)
                              gCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTAGGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCg                                                                                                                       <  1:1234334/150‑1 (MQ=255)
                               ccACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGCCGACGGTGTCGACGTCTCGGcgccacgg                                                                                                                      >  2:855188/1‑143 (MQ=255)
                                cACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCggg                                                                                                                     >  2:917161/1‑150 (MQ=255)
                                  cTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCa                                                                                                                   <  1:438383/150‑1 (MQ=255)
                                        gCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCNCGCANGAGNTCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTAGGCCCTCCGGGCAGCCATc                                                                                                             >  1:1495842/1‑150 (MQ=255)
                                             tGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGACGTCGACACCGTCGACGGCTCggttca                                                                                                                          >  1:954741/1‑128 (MQ=38)
                                             tGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGACGTCGACACCGTCGACGGCTCggttca                                                                                                                          <  2:954741/132‑5 (MQ=38)
                                              gCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGCGCAGCCATCAGTCAc                                                                                                       <  2:818860/150‑1 (MQ=255)
                                               ctctACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACg                                                                                                      >  1:1926383/1‑150 (MQ=255)
                                                tctACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGCCACGc                                                                                                     >  2:1706974/1‑150 (MQ=255)
                                                  tACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCtt                                                                                                   >  2:1407768/1‑150 (MQ=255)
                                                       gaggagGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCg                                                                                                                       >  2:1118899/1‑125 (MQ=255)
                                                       gaggagGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCg                                                                                                                       <  1:1118899/125‑1 (MQ=255)
                                                       gaggagGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTgcgcg                                                                                              >  1:698083/1‑150 (MQ=255)
                                                          gaggGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGccc                                                                                           <  1:855188/150‑1 (MQ=255)
                                                             ggTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAg                                                                                                                  >  2:1096925/1‑124 (MQ=255)
                                                             ggTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAg                                                                                                                  <  1:1096925/124‑1 (MQ=255)
                                                             ggTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCAt                                                                                                              >  1:1218874/1‑128 (MQ=255)
                                                             ggTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCAt                                                                                                              <  2:1218874/128‑1 (MQ=255)
                                                             ggTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCAc                                                                                        <  2:1495842/150‑1 (MQ=255)
                                                               tAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACtt                                                                                      <  2:950380/150‑1 (MQ=255)
                                                                aaTGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGtct                                                                                                                                 <  1:1879294/106‑1 (MQ=255)
                                                                aaTGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGtct                                                                                                                                 >  2:1879294/1‑106 (MQ=255)
                                                                                cGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTc                                                                     <  1:1187646/150‑1 (MQ=255)
                                                                                cGCATGAGATCGCCAAGCATATCTTGCTAAGGACTGCGTGCCCTATGGGACCGGGTGTCGGGATTTCATTCAGGCACAGTGTCGACGTCGCGCCCCTCCGGGCCGCCATCACTCACGCGTGCGCTCCCCACTTGCCTGCAGCAGAAATTc                                                                     <  1:1407768/150‑1 (MQ=255)
                                                                                 gCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCACTCCGGGCCGCAATCAGTAATACTTGCTCGCCCCACTTGCCTGAAGCAGAAATTCt                                                                    >  2:1911144/1‑150 (MQ=255)
                                                                                       gattgCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACt                                                                                       <  2:1236693/121‑1 (MQ=255)
                                                                                       gaTTGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACt                                                                                       >  1:1236693/1‑125 (MQ=255)
                                                                                          cGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGct                                                           >  1:881961/1‑150 (MQ=255)
                                                                                             cTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGc                                                                                                                 >  2:736515/1‑93 (MQ=255)
                                                                                             cTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGc                                                                                                                 <  1:736515/93‑1 (MQ=255)
                                                                                                         gCTAAGGACTGCCTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACCGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGTACCATAAATTCTGAACATGCTCTGTTCGACCGCTCa                                            <  1:1911144/150‑1 (MQ=39)
                                                                                                           tAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATg                                          >  2:1169072/1‑150 (MQ=255)
                                                                                                           tAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATg                                          >  1:253108/1‑150 (MQ=255)
                                                                                                            aaGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCt                                                                    >  2:1358518/1‑123 (MQ=255)
                                                                                                            aaGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCt                                                                    <  1:1358518/123‑1 (MQ=255)
                                                                                                            aaGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATgg                                         >  2:360111/1‑150 (MQ=255)
                                                                                                            aaGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATgg                                         <  2:1385800/150‑1 (MQ=255)
                                                                                                               gACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGAtt                                      >  1:46090/1‑150 (MQ=255)
                                                                                                                    cGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCt                                                                    <  2:406700/115‑1 (MQ=255)
                                                                                                                    cGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCt                                                                    >  1:406700/1‑115 (MQ=255)
                                                                                                                    cGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATtgctgc                                 <  1:1169072/150‑1 (MQ=255)
                                                                                                                      tGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGt                               >  2:798517/1‑150 (MQ=255)
                                                                                                                       gCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGtt                              >  1:192239/1‑150 (MQ=255)
                                                                                                                         ccTATGGCCACGGCCGCCGACAAGCGAGCCGGCGACGCTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTg                            >  2:1197633/1‑150 (MQ=255)
                                                                                                                          catagggCACGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTgg                           <  1:1232578/146‑1 (MQ=255)
                                                                                                                          catagggCACGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTgg                           >  2:1232578/5‑150 (MQ=255)
                                                                                                                          cTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTgg                           <  1:798517/150‑1 (MQ=255)
                                                                                                                          cTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTgg                           <  2:46090/150‑1 (MQ=255)
                                                                                                                          cTATGGGCCCGGCTGTCGTGAACCGAGAAGTCGCCGGAGTCGGGCGCGCGTGCCTCCGGGCAGCCATCAGTCAGGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTgg                           <  1:917161/150‑1 (MQ=255)
                                                                                                                           tATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCATCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTggg                          <  2:1926383/150‑1 (MQ=255)
                                                                                                                           tATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCt                                                                    >  1:524841/1‑108 (MQ=255)
                                                                                                                           tATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCt                                                                    <  2:524841/108‑1 (MQ=255)
                                                                                                                              gggCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGACGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTAt                       <  2:253108/150‑1 (MQ=255)
                                                                                                                               ggCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATg                      >  1:921735/1‑150 (MQ=255)
                                                                                                                                gCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGa                     >  1:1488606/1‑150 (MQ=255)
                                                                                                                                  ccGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACt                   <  2:1925188/150‑1 (MQ=255)
                                                                                                                                       tGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGaa              >  2:1203543/1‑150 (MQ=255)
                                                                                                                                       tGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGaa              <  2:698083/150‑1 (MQ=255)
                                                                                                                                          cGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACga           <  1:1113062/150‑1 (MQ=255)
                                                                                                                                          cGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACga           <  1:1203543/150‑1 (MQ=255)
                                                                                                                                           gTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACgag          >  2:206696/1‑150 (MQ=255)
                                                                                                                                              aaCCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACa       <  2:1884818/150‑1 (MQ=255)
                                                                                                                                                 cGAGACGTCGACACCGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATtgctgc                                 <  2:164685/121‑1 (MQ=255)
                                                                                                                                                 cGAGACGTCGACACCGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATtgctgc                                 >  1:164685/1‑121 (MQ=255)
                                                                                                                                                  gagaCGTCGACACCGTCGACGGCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGa                                                   >  2:964654/1‑102 (MQ=38)
                                                                                                                                                   aGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCgaa  >  2:189267/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCTACCGTTGCCAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCTTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGAA  >  NC_002947/4196179‑4196475

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