Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,457,000 T→A 5.7% L24Q (CTG→CAG)  PP_2150 → thiamine biosynthesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,457,0000TA 5.7% 311.8 / 6.9 105L24Q (CTG→CAG) PP_2150thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (50/49);  new base A (4/2);  total (54/51)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

Only 100 of 105 total aligned reads displayed.

ACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATTCCATCG  >  NC_002947/2456855‑2457146
                                                                                                                                                   |                                                                                                                                                  
aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                                                                                                                                                 >  7:38934/1‑149 (MQ=255)
aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                                                                                                                                                 >  2:159448/1‑149 (MQ=255)
aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                                                                                                                                                 >  8:142231/1‑149 (MQ=255)
 cAGGATTTACGCTGCCCCTGTAGGAGCTGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTgg                                                                                                                                                  <  5:189237/147‑1 (MQ=255)
 cAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                                                                                                                                                 <  8:38934/148‑1 (MQ=255)
   ggATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtgg                                                                                                                                              <  1:91108/149‑1 (MQ=255)
    gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtg                                                                                                                                               >  6:288463/1‑147 (MQ=255)
    gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtg                                                                                                                                               <  5:288463/147‑1 (MQ=255)
    gATTATCGCTGCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtggt                                                                                                                                             <  3:294619/149‑1 (MQ=255)
        aTCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGa                                                                                                                                         >  4:38520/1‑149 (MQ=255)
         tCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAc                                                                                                                                        >  8:246213/1‑149 (MQ=255)
         tCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAc                                                                                                                                        <  8:135918/149‑1 (MQ=255)
             taatccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATg                                                                                                                                    <  6:258136/145‑1 (MQ=255)
               ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg                                                                                                                                  <  3:38520/149‑1 (MQ=255)
               ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg                                                                                                                                  <  1:159448/149‑1 (MQ=255)
               ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg                                                                                                                                  >  6:87516/1‑149 (MQ=255)
               ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTACAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATg                                                                                                                                    <  5:140707/147‑1 (MQ=255)
                cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgc                                                                                                                                 <  5:87516/149‑1 (MQ=255)
                 ccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg                                                                                                                                  >  4:53694/1‑147 (MQ=255)
                  cTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGa                                                                                                                                         <  5:346464/139‑1 (MQ=255)
                  cTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGa                                                                                                                                         >  6:346464/1‑139 (MQ=255)
                  cTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgcag                                                                                                                               <  8:83676/149‑1 (MQ=255)
                  cTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgcag                                                                                                                               <  7:142231/149‑1 (MQ=255)
                     tcggAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgcagca                                                                                                                             <  1:80400/146‑1 (MQ=255)
                     tAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgcagca                                                                                                                             <  5:145450/148‑1 (MQ=255)
                     tAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgcagca                                                                                                                             <  3:168010/148‑1 (MQ=255)
                     tAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgcagca                                                                                                                             <  3:53694/148‑1 (MQ=255)
                      tcggcggggTTTACCCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgcagca                                                                                                                             <  5:54266/142‑1 (MQ=255)
                         aGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtg                                                                                                                                               >  8:304147/1‑126 (MQ=255)
                         aGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtg                                                                                                                                               <  7:304147/126‑1 (MQ=255)
                           cGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCGCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACa                                                                                                                       <  1:73039/149‑1 (MQ=255)
                           cGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                                                                                                                                                 >  2:49887/1‑122 (MQ=255)
                           cGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                                                                                                                                                 <  1:49887/122‑1 (MQ=255)
                           cGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAg                                                                                                                      >  6:355892/1‑149 (MQ=255)
                              gTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGTGTGGTCCGACCAt                                                                                                                                     <  8:200986/131‑1 (MQ=255)
                              gTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCAt                                                                                                                                     >  7:200986/1‑131 (MQ=255)
                               tttACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATc                                                                                                                  <  7:246213/149‑1 (MQ=255)
                               tttACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATc                                                                                                                  >  4:67314/1‑149 (MQ=255)
                                  aCCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCa                                                                                                                <  6:84492/148‑1 (MQ=255)
                                    cccGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAg                                                                                                               <  3:67314/148‑1 (MQ=255)
                                     ccGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTa                                                                                                             >  8:3921/1‑149 (MQ=255)
                                       gcgAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTAGTGTTCTTTCTGCTGGTAATCCCCGCGTGCACCCACGGCCCCCCCCTGGCGCGCCTGGGTGGTCCGCCAATGGGCAGCAGCTACAGCATCCAGTACg                                                                                                           >  6:341757/1‑149 (MQ=255)
                                       gcgAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAc                                                                                                                                        <  2:191551/120‑1 (MQ=255)
                                             agGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTa                                                                                                             <  6:135745/141‑1 (MQ=255)
                                             agGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTa                                                                                                             >  5:135745/1‑141 (MQ=255)
                                               gCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGcagcag                                                                                                                            >  1:216419/1‑124 (MQ=255)
                                                 cAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACa                                                                                                                       <  1:121391/127‑1 (MQ=255)
                                                 cAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACa                                                                                                                       >  2:121391/1‑127 (MQ=255)
                                                  aGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCt                                                                                                >  1:136489/1‑149 (MQ=255)
                                                   gCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCtg                                                                                               >  2:178055/1‑149 (MQ=255)
                                                   gCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCtg                                                                                               <  8:224047/149‑1 (MQ=255)
                                                       cGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCtggt                                                                                             >  7:200452/1‑147 (MQ=255)
                                                       cGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCtggt                                                                                             <  8:200452/147‑1 (MQ=255)
                                                       cGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCACCAGCTACAGCATCCAGTACGTCCGCGACCCtggtcg                                                                                           >  7:262886/1‑147 (MQ=255)
                                                                agaAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGcc                                                                                  <  7:9099/149‑1 (MQ=255)
                                                                agaAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGcc                                                                                  <  7:3921/149‑1 (MQ=255)
                                                                 gaAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCg                                                                                 <  5:362967/149‑1 (MQ=255)
                                                                               cgcACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGa                                                                    <  5:355892/148‑1 (MQ=255)
                                                                                   catcccTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCg                                                               <  5:341757/146‑1 (MQ=255)
                                                                                     gCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGt                                                             <  4:104637/149‑1 (MQ=255)
                                                                                             gTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGt                                                                         <  2:33717/129‑1 (MQ=255)
                                                                                             gTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGt                                                                         >  1:33717/1‑129 (MQ=255)
                                                                                              ttCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCtggt                                                                                             <  1:341691/108‑1 (MQ=255)
                                                                                              ttCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGCGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCtggt                                                                                             >  2:341691/1‑108 (MQ=255)
                                                                                                cATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATc                                                   <  1:178055/148‑1 (MQ=255)
                                                                                                      gctgctACTCACCGCCTGCAACCAAGGCCCCACCCTCGAGCGCCTGCGTGGTCCGACCATGGGCAGCAGCAACAGCAACCAGTAAGTCCCCGAGCCTGGTCGCCCGCTCCCGCCCCACGAGCAGATCGCGGTGTAACCCAACCTGCAc                                             >  5:148767/1‑148 (MQ=255)
                                                                                                           tACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATc                                       >  5:201346/1‑149 (MQ=255)
                                                                                                            aCTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATc                                                   <  3:251055/136‑1 (MQ=255)
                                                                                                             cTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGa                                     >  3:76182/1‑149 (MQ=255)
                                                                                                              tCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTcc                                                                                                        >  3:270272/1‑81 (MQ=255)
                                                                                                              tCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTcc                                                                                                        <  4:270272/81‑1 (MQ=255)
                                                                                                              tCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGac                                    >  7:302998/1‑149 (MQ=255)
                                                                                                              tCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGac                                    >  6:169151/1‑149 (MQ=255)
                                                                                                                aCCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACGGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGacca                                  >  1:119025/1‑149 (MQ=255)
                                                                                                                aCCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGacca                                  >  6:56720/1‑149 (MQ=255)
                                                                                                                aCCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGac                                    <  8:302998/147‑1 (MQ=255)
                                                                                                                 ccGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGACCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAAc                                                       >  8:43355/1‑127 (MQ=255)
                                                                                                                 ccGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGACCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAAc                                                       <  7:43355/127‑1 (MQ=255)
                                                                                                                  cGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGAccacca                               >  3:307981/1‑149 (MQ=255)
                                                                                                                      tGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGa                                                                    >  8:306519/1‑109 (MQ=255)
                                                                                                                       gCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATc                                                   >  2:121132/1‑125 (MQ=255)
                                                                                                                         aaCCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGaccac                                 <  3:221581/141‑1 (MQ=255)
                                                                                                                         aaCCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGaccac                                 >  4:221581/1‑141 (MQ=255)
                                                                                                                            cAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGAc                      >  4:37639/1‑149 (MQ=255)
                                                                                                                            cAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGAc                      >  7:275646/1‑149 (MQ=255)
                                                                                                                             aaGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGAcc                     >  1:333603/1‑149 (MQ=255)
                                                                                                                                 ggaccACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTAc                  >  3:153641/4‑149 (MQ=255)
                                                                                                                                  gaccACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTAcc                 >  8:187622/3‑149 (MQ=255)
                                                                                                                                  gaccACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTAcc                 >  8:184897/3‑149 (MQ=255)
                                                                                                                                  cccACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAAGCATCCTGCACGAc                                          >  3:227824/1‑123 (MQ=255)
                                                                                                                                  cccACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGAc                                          <  4:227824/123‑1 (MQ=255)
                                                                                                                                    ccACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACcgc               >  2:296194/1‑149 (MQ=255)
                                                                                                                                           ggAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGAtt        >  7:251553/1‑149 (MQ=255)
                                                                                                                                            gAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACGGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGAtt        <  2:119025/148‑1 (MQ=255)
                                                                                                                                            gAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATTc       >  3:271412/1‑149 (MQ=255)
                                                                                                                                              gcgcCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATTCCa     <  8:275646/149‑1 (MQ=255)
                                                                                                                                               cgcCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATTCCAt    <  5:56720/149‑1 (MQ=255)
                                                                                                                                                  cTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATTCCATCg  >  5:56291/1‑148 (MQ=255)
                                                                                                                                                  cTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATTCCATCg  >  5:51116/1‑148 (MQ=255)
                                                                                                                                                  cTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATTCCATCg  >  1:271749/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATTCCATCG  >  NC_002947/2456855‑2457146

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: