Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 126,115 A→G 5.9% intergenic (+25/+25) PP_0120 → / ← thrB zinc ABC transporter/homoserine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947126,1150AG 5.9% 242.8 / 4.6 84intergenic (+25/+25)PP_0120/thrBzinc ABC transporter/homoserine kinase
Reads supporting (aligned to +/- strand):  ref base A (42/37);  new base G (2/3);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

GTTGCCGGTACGCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGGCGCGA  >  NC_002947/125967‑126262
                                                                                                                                                    |                                                                                                                                                   
gTTGCCGGTACGCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCa                                                                                                                                                     <  3:245089/149‑1 (MQ=255)
    ccGGTACGCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGc                                                                                                                                                  <  7:208153/148‑1 (MQ=255)
     cGGTACGCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGccc                                                                                                                                                >  2:252231/1‑149 (MQ=255)
      ggTACGCCTTGCCGAGCTGGATGCCTTGGGGGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCcg                                                                                                                                               <  3:123200/149‑1 (MQ=255)
       gTACGCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCcgc                                                                                                                                              <  1:266097/149‑1 (MQ=255)
       gTACGCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCcgc                                                                                                                                              >  6:302029/1‑149 (MQ=255)
          cGCCTTGCCGAGCTGGATGCCTTGGGGGGCACCGATCCTGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCGTGCCCGCgaa                                                                                                                                           <  2:354457/149‑1 (MQ=255)
           gCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCcgcg                                                                                                                                             >  2:157533/1‑146 (MQ=255)
           gCCTAGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCcgcg                                                                                                                                             >  4:288950/1‑146 (MQ=255)
            ccTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaaga                                                                                                                                         <  1:157533/149‑1 (MQ=255)
            ccTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCCGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaaga                                                                                                                                         <  7:56110/149‑1 (MQ=255)
             cTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCGTGCCCGCgaaggg                                                                                                                                        <  5:152820/149‑3 (MQ=255)
             cTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaag                                                                                                                                          >  5:233873/1‑147 (MQ=255)
                gCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGc                                                                                                                                      <  8:289096/148‑1 (MQ=255)
                gCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGCAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGc                                                                                                                                      <  3:43699/148‑1 (MQ=255)
                 ccGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAg                                                                                                                                    <  8:91846/149‑1 (MQ=255)
                 ccGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACCCCCCCgaagaac                                                                                                                                       >  3:3786/1‑145 (MQ=255)
                   gAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGAc                                                                                                                                  >  2:23045/1‑149 (MQ=255)
                      cTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGAcgcg                                                                                                                               >  8:307888/1‑149 (MQ=255)
                         gATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCgg                                                                                                                              <  1:130149/147‑1 (MQ=255)
                         gATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGt                                                                                                                             <  1:252231/148‑1 (MQ=255)
                          aTGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTc                                                                                                                            >  4:226108/1‑148 (MQ=255)
                                   gTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaa                                                                                                                                           <  5:333754/124‑1 (MQ=255)
                                   gTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaa                                                                                                                                           >  6:333754/1‑124 (MQ=255)
                                    tGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGaa                                                                                                                 <  7:184149/149‑1 (MQ=255)
                                       cACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGACGCAGCAGACGCGGTCTACAGCGCGAACgg                                                                                                              >  4:46487/1‑149 (MQ=255)
                                        aCCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGc                                                                                                                                      <  8:219423/124‑1 (MQ=255)
                                        aCCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGc                                                                                                                                      >  7:219423/1‑124 (MQ=255)
                                        aCCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGc                                                                                                             <  3:421487/149‑1 (MQ=255)
                                        aCCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACACCGCGAACgg                                                                                                              >  6:268125/1‑148 (MQ=255)
                                          cGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTc                                                                                                                            >  1:394165/1‑132 (MQ=255)
                                          cGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTc                                                                                                                            <  2:394165/132‑1 (MQ=255)
                                          cGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAg                                                                                                           <  5:302029/149‑1 (MQ=255)
                                               cGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGc                                                                                                       <  8:265380/148‑1 (MQ=255)
                                                 gTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAAt                                                                                                    >  8:301896/1‑149 (MQ=255)
                                                    gATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGc                                                                                                  >  5:401572/1‑148 (MQ=255)
                                                             gggTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACAGCAGGTGGCATTGCACctgccgt                                                                                        <  7:290869/149‑4 (MQ=255)
                                                                   gAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCCGACGCGGTCTACAGCGCGAACGGCAGGTGCCATGCCACCTGCTGGCGCTgc                                                                                  >  2:418302/1‑149 (MQ=255)
                                                                     gcgcTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTgcgc                                                                                >  6:362397/1‑149 (MQ=255)
                                                                     gcgcTTGCTGGAGAAGCTGGGCGGTGACCAGGCAGGCTGCCTGGAGCAGTTGTAAGCCAGTGCCGGCCACATAGCGGGCACGCCCGCCAACAAGCAGACGCGGTCTACAGCGCGAACGGCACGTg                                                                                                        >  3:29347/1‑125 (MQ=255)
                                                                            cTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGc                                                                         <  2:71700/149‑1 (MQ=255)
                                                                              ggAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCAc                                                                       >  6:122541/1‑149 (MQ=255)
                                                                              ggAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGAGGGCACGCCCGCGAAGAAGCAGCCGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTCCTCGCGCTGCGCCAGACGCa                                                                        >  5:364094/1‑148 (MQ=255)
                                                                               gagaAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCa                                                                        <  7:307888/147‑1 (MQ=255)
                                                                                agaAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACtt                                                                     <  3:375961/149‑1 (MQ=255)
                                                                                agaAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACtt                                                                     >  8:90113/1‑149 (MQ=255)
                                                                                agaAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACtt                                                                     <  3:288950/149‑1 (MQ=255)
                                                                                     cTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAAc                                                                >  5:404477/1‑149 (MQ=255)
                                                                                       gggCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAAc                                                                >  1:105156/1‑147 (MQ=255)
                                                                                       gggCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAAc                                                                <  2:105156/147‑1 (MQ=255)
                                                                                           ggTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTgcgc                                                                                <  6:205130/127‑1 (MQ=255)
                                                                                           ggTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTgcgc                                                                                >  5:205130/1‑127 (MQ=255)
                                                                                           ggTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAActcgctc                                                          <  7:199784/149‑1 (MQ=255)
                                                                                                ccTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATc                                                     >  5:348396/1‑149 (MQ=255)
                                                                                                 cTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAActcgctc                                                          <  7:296225/143‑1 (MQ=255)
                                                                                                 cTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAActcgctc                                                          >  8:296225/1‑143 (MQ=255)
                                                                                                    gCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTgg                                                 <  6:401572/149‑1 (MQ=255)
                                                                                                      aGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGat                                               >  3:377891/1‑149 (MQ=255)
                                                                                                                   cAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCAAACTCTCTCGGAACGTGGCTAATCAcatcaatc                                  >  4:365772/1‑149 (MQ=255)
                                                                                                                     gTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCCTCACATCCATccc                                >  4:282159/1‑149 (MQ=255)
                                                                                                                        gTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCCTCACATCCATCCCGGc                             >  6:6353/1‑149 (MQ=255)
                                                                                                                           aGCGTGTGCCGTCCCCTTCGCGGGCACGCTCGCGAAGATGCAGACGGGGTCTACAGCGGGAAGGGCAGGTTCAATGCGACGTGCTGTCGATGCGCCAGACTCTCTTCGATCTCGCTCGGATCGTGGATCATCACTTCCATCCCGGGGaa                          <  6:364094/149‑1 (MQ=255)
                                                                                                                            gCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGAtcatca                                          <  8:89839/132‑1 (MQ=255)
                                                                                                                            gCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGAtcatca                                          >  7:89839/1‑132 (MQ=255)
                                                                                                                            gCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAAc                         <  6:407019/149‑1 (MQ=255)
                                                                                                                            gCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAAc                         <  6:9640/149‑1 (MQ=255)
                                                                                                                             ccTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACg                        >  4:78684/1‑149 (MQ=255)
                                                                                                                              cTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGa                       <  1:23045/149‑1 (MQ=255)
                                                                                                                               tgtgCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCAcat                                       <  2:286705/132‑1 (MQ=255)
                                                                                                                               tgtgCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCAcat                                       >  1:286705/1‑132 (MQ=255)
                                                                                                                               tgtgCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGa                       >  4:203345/1‑148 (MQ=255)
                                                                                                                                 tgCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTc                    <  5:122541/149‑1 (MQ=255)
                                                                                                                                    cGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGCAGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATAGTGGATCATACCATCCATCCCGGCGAACGACTAGGc                 >  6:291839/1‑149 (MQ=255)
                                                                                                                                    cGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGc                 >  4:411402/1‑149 (MQ=255)
                                                                                                                                      gCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGAt                                                      <  5:422506/110‑1 (MQ=255)
                                                                                                                                      gCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGAt                                                      >  6:422506/1‑110 (MQ=255)
                                                                                                                                          cTTCGCGGGCGTGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGt                                                   <  1:43692/109‑1 (MQ=255)
                                                                                                                                          cTTCGCGGGCGTGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGt                                                   >  2:43692/1‑109 (MQ=255)
                                                                                                                                           ttCGCGGGCGTGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAAt           >  6:387380/1‑148 (MQ=255)
                                                                                                                                            tCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAAt           >  4:340949/1‑147 (MQ=255)
                                                                                                                                              gcgGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAgg       <  7:301896/149‑1 (MQ=255)
                                                                                                                                                gggCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGgcg     >  4:150034/1‑149 (MQ=255)
                                                                                                                                                gggCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGgcg     <  3:340949/149‑1 (MQ=255)
                                                                                                                                                  gCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCACTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGgcgc    >  6:10446/1‑148 (MQ=255)
                                                                                                                                                   cACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGGCGCGa  >  2:19822/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTTGCCGGTACGCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGGCGCGA  >  NC_002947/125967‑126262

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