Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,415,032 T→C 5.6% intergenic (+141/+27) PP_5644 → / ← PP_3899 membrane protein/hicB family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,415,0320TC 5.6% 270.9 / 4.4 89intergenic (+141/+27)PP_5644/PP_3899membrane protein/hicB family protein
Reads supporting (aligned to +/- strand):  ref base T (42/42);  new base C (3/2);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

GTAGATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATTG  >  NC_002947/4414888‑4415179
                                                                                                                                                |                                                                                                                                                   
gTAGATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCa                                                                                                                                                 <  5:337722/149‑1 (MQ=255)
 tAGATGGCTAAGAAGCACACCTCTTCAAGCCCCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCACCCCTGCACAGTGAAGGATCATCCGCTTCCCTACGTTAGCCAGCATTGCAg                                                                                                                                                >  8:236703/1‑149 (MQ=255)
   gATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTg                                                                                                                                              >  3:364060/1‑149 (MQ=255)
   gATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTg                                                                                                                                              >  7:289857/1‑149 (MQ=255)
        cTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCACTGCAATGCGGGc                                                                                                                                         <  5:223169/149‑1 (MQ=255)
         tATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCt                                                                                                                                        >  5:393848/1‑149 (MQ=255)
           tGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCACTGCAATGCGGGCta                                                                                                                                       <  2:58221/148‑2 (MQ=255)
             aaGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCtttt                                                                                                                                     <  8:52002/148‑1 (MQ=255)
                 aCTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTg                                                                                                                                <  1:282943/149‑1 (MQ=255)
                  ctcctcTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTgg                                                                                                                               <  5:234804/149‑1 (MQ=255)
                   tcctcTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTggg                                                                                                                              >  5:315006/1‑149 (MQ=255)
                    cctcTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGggc                                                                                                                             >  1:61571/1‑149 (MQ=255)
                     ctctTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTggg                                                                                                                              <  8:296124/147‑1 (MQ=255)
                     ctctTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTACCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGggcg                                                                                                                            >  7:329394/1‑149 (MQ=255)
                     ctctTGAAGCCGCCATCAAGTGATCCCACAGGCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGggcg                                                                                                                            <  8:164915/149‑1 (MQ=255)
                            aGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTAcgc                                                                                                                     >  2:211719/1‑149 (MQ=255)
                                gccATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGtga                                                                                                                 >  7:126177/1‑149 (MQ=255)
                                gccATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGtga                                                                                                                 <  4:364060/149‑1 (MQ=255)
                                 ccATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGtgat                                                                                                                >  5:93696/1‑149 (MQ=255)
                                   aTCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCAGGATGCATTCACGCTTTTTTCACGCCCTCCCCTGCAGAGTGAAGGCTCATCCGCTTCCCTACGTATGCCCGCTTTGCAGTGCGGGCTTTTGTTTGGGCGGCTACGCGtgatga                                                                                                              <  7:236703/149‑1 (MQ=255)
                                   aTCAAGTGATCCCACAGTCCAGCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGGTCATCCGCTTACCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGttatga                                                                                                              <  5:132135/149‑1 (MQ=255)
                                     cAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGAcc                                                                                                            >  3:170343/1‑149 (MQ=255)
                                        gTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAg                                                                                                         <  6:315006/149‑1 (MQ=255)
                                        gTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAg                                                                                                         >  7:281523/1‑149 (MQ=255)
                                         tGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAgg                                                                                                        >  1:375145/1‑149 (MQ=255)
                                          gATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTACCCAGCACAGAGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCTGCTACGCGTGATGACCCAgg                                                                                                        >  5:266142/1‑148 (MQ=255)
                                              ccACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGCCCCAGGTACt                                                                                                    >  8:68160/1‑148 (MQ=255)
                                              ccACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACt                                                                                                    >  5:375612/1‑148 (MQ=255)
                                               cacaGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCa                                                                                                  <  8:99792/149‑1 (MQ=255)
                                               cacaGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCa                                                                                                  >  6:9786/1‑149 (MQ=255)
                                                acaGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCa                                                                                                                                                 <  5:266495/101‑1 (MQ=255)
                                                acaGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGCCCCAGGTACTCAc                                                                                                 >  6:115603/1‑149 (MQ=255)
                                                acaGTCCGCCTCGGCGGGCTTAACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCa                                                                                                                                                 >  6:266495/1‑101 (MQ=255)
                                                   gTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTaa                                                                                              <  6:375612/149‑1 (MQ=255)
                                                     ccGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAgc                                                                                            >  3:70664/1‑149 (MQ=255)
                                                      cGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCGTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAgcg                                                                                           >  2:146180/1‑149 (MQ=255)
                                                       gCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAgcgc                                                                                          >  6:132461/1‑149 (MQ=255)
                                                        ccTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCt                                                                                         >  6:326878/1‑149 (MQ=255)
                                                          tCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCttt                                                                                       >  1:241520/1‑149 (MQ=255)
                                                            ggcggGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCttt                                                                                       <  8:289857/147‑1 (MQ=255)
                                                            ggcggGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGc                                                                                     <  6:103176/149‑1 (MQ=255)
                                                              cggGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCgt                                                                                   <  5:326878/149‑1 (MQ=255)
                                                                       ccATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCa                                                                          >  3:9126/1‑149 (MQ=255)
                                                                         aTGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCAtt                                                                        <  2:91783/149‑1 (MQ=255)
                                                                            cATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCaa                                                                     <  8:126177/149‑1 (MQ=255)
                                                                            cATTCACGCGTTTATCACCCCCTCCCCTGCACAGGGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGTGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCaa                                                                     <  8:329394/149‑1 (MQ=255)
                                                                             aTTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAg                                                                    >  4:50005/1‑149 (MQ=255)
                                                                             aTTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAg                                                                    >  8:86075/1‑149 (MQ=255)
                                                                             aTTAACGCTTTTTTCACGCCCTCCCCTGCACAGTGAACGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAg                                                                    >  8:224198/1‑149 (MQ=255)
                                                                                cACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGa                                                                  <  5:9786/148‑1 (MQ=255)
                                                                                 aCGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGTGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGaaa                                                                >  6:243868/1‑149 (MQ=255)
                                                                                 aCGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGaaa                                                                >  6:162621/1‑149 (MQ=255)
                                                                                  cGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAAt                                                               <  1:211719/149‑1 (MQ=255)
                                                                                   gCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTAcgc                                                                                                                     >  8:84926/1‑94 (MQ=255)
                                                                                   gCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATa                                                              >  5:309199/1‑149 (MQ=255)
                                                                                    cTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATAc                                                             >  2:2083/1‑149 (MQ=255)
                                                                                          tCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGc                                                        <  2:241520/148‑1 (MQ=255)
                                                                                          tCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGc                                                        <  4:62598/148‑1 (MQ=255)
                                                                                              gCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTtggt                                                   <  1:146180/149‑1 (MQ=255)
                                                                                              gCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTtggt                                                   >  7:50135/1‑149 (MQ=255)
                                                                                               cccTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTtggtg                                                  <  4:265946/149‑1 (MQ=255)
                                                                                                ccTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTtggtgg                                                 <  1:350303/149‑1 (MQ=255)
                                                                                                  tCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCt                                               >  6:237327/1‑149 (MQ=255)
                                                                                                    cccTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCt                                               <  8:281523/147‑1 (MQ=255)
                                                                                                    cccTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGc                                             <  2:61571/149‑1 (MQ=255)
                                                                                                            aGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTgg                                     <  7:224198/149‑1 (MQ=255)
                                                                                                             gTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGc                                    <  7:86075/149‑1 (MQ=255)
                                                                                                                aaGGCTCATCCGCTTCCCTACGTTAGCCCTCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCg                                  <  4:220100/148‑1 (MQ=255)
                                                                                                                  ggCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTAGTGACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTa                               <  3:303556/149‑1 (MQ=255)
                                                                                                                  ggCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTa                               <  3:308769/149‑1 (MQ=255)
                                                                                                                       aTCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGt                            <  6:309199/147‑1 (MQ=255)
                                                                                                                           gCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGt                      >  2:286037/1‑149 (MQ=255)
                                                                                                                           gCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGt                      >  2:37992/1‑149 (MQ=255)
                                                                                                                               cccTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGt                   <  3:121378/148‑1 (MQ=255)
                                                                                                                                ccTACGTTAGCCCGGATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGc                 <  6:393848/149‑1 (MQ=255)
                                                                                                                                ccTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGTGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGc                 <  5:243868/149‑1 (MQ=255)
                                                                                                                                ccTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTg                  <  3:211107/148‑1 (MQ=255)
                                                                                                                                ccTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGc                 >  5:326719/1‑149 (MQ=255)
                                                                                                                                ccTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGc                 <  7:68160/149‑1 (MQ=255)
                                                                                                                                 cTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGcc                <  5:115603/149‑1 (MQ=255)
                                                                                                                                 cTACGTTAGCCCGCATTGCAGTGCGGGCTTATCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGcc                >  1:34782/1‑149 (MQ=255)
                                                                                                                                 aaagaaaaGCCCGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGcc                >  7:125580/8‑149 (MQ=255)
                                                                                                                                      aaaGCCCGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTAcc           >  6:34292/3‑149 (MQ=255)
                                                                                                                                         gCCCGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCt         >  3:255094/1‑148 (MQ=255)
                                                                                                                                          cccGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGt                   >  3:290397/1‑137 (MQ=255)
                                                                                                                                          cccGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGt                   <  4:290397/137‑1 (MQ=255)
                                                                                                                                          cccGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTca       >  6:244826/1‑149 (MQ=255)
                                                                                                                                             gCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACAt    >  1:363990/1‑149 (MQ=255)
                                                                                                                                               aTTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATTg  <  5:132461/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GTAGATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCTCACATTG  >  NC_002947/4414888‑4415179

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