Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,653,983 T→C 7.7% Q28R (CAG→CGG) ‡ ygjP ← metal‑dependent hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,653,9830TC 7.7% 1088.1 / 54.0 364Q28R (CAG→CGG) ‡ygjPmetal‑dependent hydrolase
Reads supporting (aligned to +/- strand):  ref base T (165/168);  new base C (24/4);  total (191/173)
Fisher's exact test for biased strand distribution p-value = 2.45e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.56e-01

Only 100 of 371 total aligned reads displayed.

GGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAAGCC  >  NC_002947/5653834‑5654132
                                                                                                                                                     |                                                                                                                                                     
ggTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCt                                                                                                                                                       <  1:6562379/150‑1 (MQ=255)
ggTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCt                                                                                                                                                       <  1:644476/150‑1 (MQ=255)
           cTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCggg                                                                                                                                            <  2:239490/150‑1 (MQ=255)
           cTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCggg                                                                                                                                            <  2:5425385/150‑1 (MQ=255)
            ttCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCgggg                                                                                                                                           >  1:4215470/1‑150 (MQ=255)
            ttCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCgggg                                                                                                                                           >  2:1441789/1‑150 (MQ=255)
            ttCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCgggg                                                                                                                                           >  1:542386/1‑150 (MQ=255)
             tCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCggggg                                                                                                                                          <  2:1769794/150‑1 (MQ=255)
              cGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCgggggg                                                                                                                                         <  1:5384324/150‑1 (MQ=255)
                     aCTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCtt                                                                                                                                  <  2:2577651/150‑1 (MQ=255)
                     aCTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCtt                                                                                                                                  <  2:3707629/150‑1 (MQ=255)
                       tGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGc                                                                                                                                >  1:2618475/1‑150 (MQ=255)
                       tGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGAAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGc                                                                                                                                >  2:1182880/1‑150 (MQ=255)
                         ttGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCaa                                                                                                                              <  2:5205383/150‑1 (MQ=255)
                         ttGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCaa                                                                                                                              <  1:6345601/150‑1 (MQ=255)
                            gcgcAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTa                                                                                                                           <  1:7951779/150‑1 (MQ=255)
                             cgcAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTaa                                                                                                                          <  1:6204008/150‑1 (MQ=255)
                              gcAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAAc                                                                                                                         >  1:1073039/1‑150 (MQ=255)
                              gcAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAAc                                                                                                                         >  2:5924722/1‑150 (MQ=255)
                              gcAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTCCCGCACCTGCAGGCGCAGGTGc                                                                                                                                               <  1:1773110/128‑1 (MQ=255)
                              gcAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGACATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAAc                                                                                                                         >  2:5710314/1‑150 (MQ=255)
                                aGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACgt                                                                                                                       <  2:7863604/150‑1 (MQ=255)
                                       gggCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCAGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACgg                                                                                                                >  2:3263819/1‑150 (MQ=255)
                                         gCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTc                                                                                                              >  2:7112079/1‑150 (MQ=255)
                                            tAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATg                                                                                                           <  2:8057284/150‑1 (MQ=255)
                                            tAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATg                                                                                                           <  1:2788079/150‑1 (MQ=255)
                                               ccATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGgcg                                                                                                        >  2:6320473/1‑150 (MQ=255)
                                                cATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGgcgc                                                                                                       <  1:5560323/150‑1 (MQ=255)
                                                cATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGgcgc                                                                                                       <  1:2887795/150‑1 (MQ=255)
                                                 aTACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCt                                                                                                      >  2:4701638/1‑150 (MQ=255)
                                                 aTACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCt                                                                                                      >  1:3111786/1‑150 (MQ=255)
                                                  tACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTg                                                                                                     >  1:4448060/1‑150 (MQ=255)
                                                    cAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTa                                                                                                   >  1:7713115/1‑150 (MQ=255)
                                                    cAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTa                                                                                                   <  2:3846813/150‑1 (MQ=255)
                                                        gcCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCg                                                                                               <  2:6546144/150‑1 (MQ=255)
                                                        gcCTTGACGCTTGGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTAGACACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTGGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCg                                                                                               <  2:2096503/150‑1 (MQ=255)
                                                          cTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGc                                                                                             <  2:700571/150‑1 (MQ=255)
                                                           ttGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCg                                                                                            >  1:1275057/1‑150 (MQ=255)
                                                              tCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTaa                                                                                         <  2:5035790/150‑1 (MQ=255)
                                                              tCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTaa                                                                                         >  1:4593535/1‑150 (MQ=255)
                                                                  tttGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGaa                                                                                     <  2:2206559/150‑1 (MQ=255)
                                                                  tttGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGaa                                                                                     >  2:2484389/1‑150 (MQ=255)
                                                                       aCGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCg                                                                                >  2:4889774/1‑150 (MQ=255)
                                                                            aTGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCt                                                                           >  2:8174358/1‑150 (MQ=255)
                                                                             tGGCGGTCAGGGTAACGGCGCTGCAGTTATTCACCCAGGCGTTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTa                                                                          <  1:4889017/150‑1 (MQ=255)
                                                                             tGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTa                                                                          <  1:1954936/150‑1 (MQ=255)
                                                                                  gTCTGGGTAACGGCGCTGCTGGTATTCACCCCGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCaa                                                                     <  1:3815727/150‑1 (MQ=255)
                                                                                    cAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATc                                                                   >  1:3615634/1‑150 (MQ=255)
                                                                                     aGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTACCGAATGCCGATTCTAGCCAATCa                                                                  >  2:5139206/1‑150 (MQ=255)
                                                                                      gggTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCAc                                                                 <  1:8427839/150‑1 (MQ=255)
                                                                                       ggTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACg                                                                >  1:226851/1‑150 (MQ=255)
                                                                                             ggCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTg                                                          <  2:6403233/150‑1 (MQ=255)
                                                                                             ggCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTg                                                          >  2:6605792/1‑150 (MQ=255)
                                                                                                  tGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAg                                                     >  2:4311008/1‑150 (MQ=255)
                                                                                                    cAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCCATCACGGCCCTGGCCAggc                                                   >  1:6807646/1‑149 (MQ=255)
                                                                                                    cAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGga                                                   <  2:195311/150‑1 (MQ=255)
                                                                                                      ggTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGgaga                                                 <  2:6106228/150‑1 (MQ=255)
                                                                                                          ttCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAAt                                             >  1:8761166/1‑150 (MQ=255)
                                                                                                          ttCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAAt                                             <  1:6028485/150‑1 (MQ=255)
                                                                                                           tCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTAGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAAtt                                            >  1:2369659/1‑150 (MQ=255)
                                                                                                            cACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCAGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTg                                           >  1:6901658/1‑150 (MQ=255)
                                                                                                             aCCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGc                                          >  2:8586524/1‑150 (MQ=255)
                                                                                                             aCCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGc                                          >  2:8382945/1‑150 (MQ=255)
                                                                                                             aCCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGc                                          >  1:2562483/1‑150 (MQ=255)
                                                                                                             aCCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCGTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGc                                          >  1:3346365/1‑150 (MQ=255)
                                                                                                              cccAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGcc                                         <  1:949271/150‑1 (MQ=255)
                                                                                                              cccAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGcc                                         >  2:3902301/1‑150 (MQ=255)
                                                                                                               ccAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCg                                        >  1:3687262/1‑150 (MQ=255)
                                                                                                                     ggTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGc                                  >  2:4134582/1‑150 (MQ=255)
                                                                                                                       tCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCg                                >  1:2123590/1‑150 (MQ=255)
                                                                                                                         gCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCt                              >  2:6855545/1‑150 (MQ=255)
                                                                                                                          cTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTa                             <  1:3263819/150‑1 (MQ=255)
                                                                                                                           tGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAg                            >  1:3178546/1‑150 (MQ=255)
                                                                                                                            gTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGc                           <  2:2850042/150‑1 (MQ=255)
                                                                                                                            gTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGc                           >  2:4444314/1‑150 (MQ=255)
                                                                                                                             tCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCg                          >  2:2690300/1‑150 (MQ=255)
                                                                                                                             tCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCg                          >  1:1003162/1‑150 (MQ=255)
                                                                                                                               gATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTc                        >  1:8618622/1‑150 (MQ=255)
                                                                                                                                atcatcTGCCGCACCTGCGCCTGCCGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCa                       <  1:2210439/150‑1 (MQ=255)
                                                                                                                                atcatcTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCa                       <  1:7531148/150‑1 (MQ=255)
                                                                                                                                atcatcTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCa                       <  1:3431731/150‑1 (MQ=255)
                                                                                                                                    tcTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTcg                   >  2:1097119/1‑150 (MQ=255)
                                                                                                                                    tcTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTcg                   <  2:1557941/150‑1 (MQ=255)
                                                                                                                                    ccccccGCACCTGCAGGCGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTcg                   >  1:4121253/5‑150 (MQ=255)
                                                                                                                                     cTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTcgc                  >  2:7666280/1‑150 (MQ=255)
                                                                                                                                      tGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTcgcg                 <  1:4632968/150‑1 (MQ=255)
                                                                                                                                      tGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTcgcg                 <  1:5733990/150‑1 (MQ=255)
                                                                                                                                       cccGCACCTGCAGGCGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCgg                >  2:2558508/2‑150 (MQ=255)
                                                                                                                                       cccGCACCTGCAGGCGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCgg                >  1:8080845/2‑150 (MQ=255)
                                                                                                                                        ccGCACCTGCAGGCGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGt               >  1:3940509/1‑150 (MQ=255)
                                                                                                                                        ccGCACCTGCAGGCGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATCGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGt               >  1:7070251/1‑150 (MQ=255)
                                                                                                                                            aCCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATc           <  2:4114/150‑1 (MQ=255)
                                                                                                                                             ccTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCg          <  1:2493891/150‑1 (MQ=255)
                                                                                                                                                  gcCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATaa     <  1:5924722/150‑1 (MQ=255)
                                                                                                                                                   ccTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAAg    >  2:7334816/1‑150 (MQ=255)
                                                                                                                                                    cTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAAGc   >  2:358428/1‑150 (MQ=255)
                                                                                                                                                    cTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAAGc   >  2:2934862/1‑150 (MQ=255)
                                                                                                                                                    cTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAAGc   <  2:5557463/150‑1 (MQ=255)
                                                                                                                                                     tGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAAGcc  >  1:1596993/1‑150 (MQ=255)
                                                                                                                                                     tGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAAGcc  >  1:5034060/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTTTGGCGCACTTCGCAGGTACTGTTGGCGCAGATCCTGGGCATAGCCATACAGCGCCTTGTCGCTTTGCACGTCATGGCGGTCAGGGTAACGGCGCTGCAGGTATTCACCCAGGCGGTCGCTGTCGATCATCTGCCGCACCTGCGCCTGCAGGTGCGGGGGGTAGGCTTGCAAGTAACGTAGTACGGTCATGGCGCTGTATTCGGCGTAACGAATGCCGATTCTAGCCAATCACGGCCCTGGCCAGGAGAAAATTGCCGGGAAGCCGCTAGCGTCAGTCGCGGTCATCGGATAAGCC  >  NC_002947/5653834‑5654132

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