Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 126,116 C→T 5.9% intergenic (+26/+24) 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,1160CT 5.9% 255.4 / 3.3 84intergenic (+26/+24)PP_0120/thrBzinc ABC transporter/homoserine kinase
Reads supporting (aligned to +/- strand):  ref base C (42/37);  new base T (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 = 9.63e-01

CCGGTACGCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGGCGCGA  >  NC_002947/125971‑126262
                                                                                                                                                 |                                                                                                                                                  
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                                                                                                                                      >  7:219423/1‑124 (MQ=255)
                                    aCCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGc                                                                                                                                      <  8:219423/124‑1 (MQ=255)
                                    aCCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGc                                                                                                             <  3:421487/149‑1 (MQ=255)
                                    aCCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACACCGCGAACgg                                                                                                              >  6:268125/1‑148 (MQ=255)
                                      cGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTc                                                                                                                            <  2:394165/132‑1 (MQ=255)
                                      cGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTc                                                                                                                            >  1:394165/1‑132 (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                                                                     >  8:90113/1‑149 (MQ=255)
                                                                            agaAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACtt                                                                     <  3:288950/149‑1 (MQ=255)
                                                                            agaAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACtt                                                                     <  3:375961/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                                                                                >  5:205130/1‑127 (MQ=255)
                                                                                       ggTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTgcgc                                                                                <  6:205130/127‑1 (MQ=255)
                                                                                       ggTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAActcgctc                                                          <  7:199784/149‑1 (MQ=255)
                                                                                            ccTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATc                                                     >  5:348396/1‑149 (MQ=255)
                                                                                             cTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAActcgctc                                                          >  8:296225/1‑143 (MQ=255)
                                                                                             cTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAActcgctc                                                          <  7:296225/143‑1 (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                                          >  7:89839/1‑132 (MQ=255)
                                                                                                                        gCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGAtcatca                                          <  8:89839/132‑1 (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                                                      >  6:422506/1‑110 (MQ=255)
                                                                                                                                  gCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGAt                                                      <  5:422506/110‑1 (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)
                                                                                                                                                 |                                                                                                                                                  
CCGGTACGCCTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGGCGCGA  >  NC_002947/125971‑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: