Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,457,034 T→A 5.4% L21Q (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,0340TA 5.4% 383.5 / 7.6 124L21Q (CTG→CAG) PP_2150thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (64/53);  new base A (3/4);  total (67/57)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

TAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATT  >  NC_002947/2456885‑2457183
                                                                                                                                                     |                                                                                                                                                     
taaccacGCTTCTCCAGGCTTATAGCAGCACTCGTAGCATCGTGTTTACTCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCt                                                                                                                                                       <  2:1191075/144‑1 (MQ=255)
tAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCt                                                                                                                                                       >  2:580204/1‑150 (MQ=255)
tAAACCCGCTCCTACAGGATCATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCt                                                                                                                                                       >  1:1957132/1‑150 (MQ=255)
     ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTg                                                                                                                                                      <  2:880709/146‑1 (MQ=255)
     ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgc                                                                                                                                                  >  2:1780796/1‑150 (MQ=255)
      cGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgcg                                                                                                                                                 >  1:1758195/1‑150 (MQ=255)
           cTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGCCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTgg                                                                                                                                            <  1:1780796/150‑1 (MQ=255)
             aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggt                                                                                                                                          >  1:215711/1‑150 (MQ=255)
              cAGGATTATTGTTTTTCCTGTACGAGCGTATTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtg                                                                                                                                         <  2:295467/150‑1 (MQ=255)
              cAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGCGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtg                                                                                                                                         >  1:620381/1‑150 (MQ=255)
               aGGATTATCGCTGCCCCTGTGGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtgg                                                                                                                                        >  2:1022439/1‑150 (MQ=255)
                   ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCg                                                                                                                                    >  1:386342/1‑150 (MQ=255)
                      tCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTg                                                                                                                                                      >  1:1021260/1‑129 (MQ=255)
                      tCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTg                                                                                                                                                      <  2:1021260/129‑1 (MQ=255)
                      tCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAcc                                                                                                                                 >  1:101827/1‑150 (MQ=255)
                      tCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCACCACCCAgg                                                                                                                                                     <  2:237419/130‑1 (MQ=255)
                        gCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtggg                                                                                                                                       >  2:173423/1‑141 (MQ=255)
                         cTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATg                                                                                                                              <  1:580204/150‑1 (MQ=255)
                          tGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATgg                                                                                                                             <  1:1418541/150‑1 (MQ=255)
                            ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCAACGGCAACCAATAAACCCTACTGTGCGCACAGCCATACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgc                                                                                                                           <  1:2271147/150‑1 (MQ=255)
                             cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCAt                                                                                                                               <  1:904514/145‑1 (MQ=255)
                             cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCAt                                                                                                                               >  2:904514/1‑145 (MQ=255)
                             cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCAACACACCTCCTGCACCCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgca                                                                                                                          >  2:1716353/1‑150 (MQ=255)
                               cTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgcagc                                                                                                                        <  1:1320262/150‑1 (MQ=255)
                                tGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg                                                                                                                            >  2:2022116/1‑145 (MQ=255)
                                tGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgcagca                                                                                                                       <  1:713017/150‑1 (MQ=255)
                                  tAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCagcagc                                                                                                                     >  1:251282/1‑150 (MQ=255)
                                   aGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTGCTGTGCATGCTGCTGCTACTCACCGCCTGAAACCAAGCCCCCACCCTGGAGCGCCTGGGTGATCCTACCATCGGCAGCAGCt                                                                                                                    >  2:346761/1‑150 (MQ=38)
                                      aGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCCGCTACa                                                                                                                 <  1:1986869/150‑1 (MQ=255)
                                        cGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGc                                                                                                               <  2:1684816/150‑1 (MQ=255)
                                        cGGGGTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGCTCATTCTGCTGCCACTCACCGCCTGCAACCAACGCCCCACCCTGGCGCGCCTGGGGGGTCCCACCATGGCCACCAGCGACAGc                                                                                                               >  2:716792/1‑150 (MQ=255)
                                            tttACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATcc                                                                                                           >  1:1158574/1‑150 (MQ=255)
                                             ttACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCa                                                                                                          <  2:697951/150‑1 (MQ=255)
                                                 ccGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCAt                                                                                                                               >  1:2114159/1‑125 (MQ=255)
                                                 ccGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTAc                                                                                                      >  1:687301/1‑150 (MQ=255)
                                                 ccGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTAc                                                                                                      >  1:807056/1‑150 (MQ=255)
                                                        gagGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                                                                                                                                           >  2:1858671/1‑106 (MQ=255)
                                                        gagGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                                                                                                                                           <  1:1858671/106‑1 (MQ=255)
                                                           gCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTa                                                                                                                   >  2:1434311/1‑127 (MQ=255)
                                                           gCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTa                                                                                                                   <  1:1434311/127‑1 (MQ=255)
                                                           gCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTcc                                                                                                  <  2:258414/144‑1 (MQ=255)
                                                             cAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATc                                                                                                            >  1:1546948/1‑132 (MQ=255)
                                                             cAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATc                                                                                                            <  2:1546948/132‑1 (MQ=255)
                                                              aGCACCGGCACCCGAGAAACCCTACTGTGCGCACAGCCCGACGGTCCATGCTGCTGCTACTAACCGCCTGCAACCGAGGCCCCCCCCTGGAGCGCCTGGGTGGTCCGACCATGGCCCGCAGCTACCGCCTCCATGCCGTCGCCGAGCCtg                                                                                         >  2:836930/1‑150 (MQ=255)
                                                                 aCCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCtggtg                                                                                      <  1:1022439/150‑1 (MQ=255)
                                                                  ccGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCcgc                                                                                                >  1:779357/1‑139 (MQ=255)
                                                                  ccGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCcgc                                                                                                <  2:779357/139‑1 (MQ=255)
                                                                    ggCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGcc                                                                                   <  1:1997936/150‑1 (MQ=255)
                                                                     gCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGccc                                                                                  >  2:1980688/1‑150 (MQ=255)
                                                                          cGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAg                                                                                                         >  2:675498/1‑122 (MQ=255)
                                                                           gagaAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTCAATCATCCCGTACGTCCGCGAGCCTGGTGGCCCGCTGcc                                                                            >  1:123977/1‑150 (MQ=255)
                                                                                aCCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGAAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGccc                                                                       >  2:1478782/1‑150 (MQ=255)
                                                                                  ccTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCcgc                                                                                                <  1:1040609/123‑1 (MQ=255)
                                                                                  ccTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCcgc                                                                                                >  2:1040609/1‑123 (MQ=255)
                                                                                  ccTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAg                                                                     <  2:2329793/150‑1 (MQ=255)
                                                                                  ccTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAg                                                                     <  2:1957132/150‑1 (MQ=255)
                                                                                    tACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGt                                                                   >  1:860088/1‑150 (MQ=255)
                                                                                      ctttgCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGc                                                                 <  1:1980688/147‑1 (MQ=255)
                                                                                      cTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTCCGGCAGCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGc                                                                 <  2:327160/150‑1 (MQ=255)
                                                                                        gtgCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAg                                                               >  1:1140243/1‑150 (MQ=255)
                                                                                         tgCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGa                                                              >  2:668230/1‑150 (MQ=255)
                                                                                         tgCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGGACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGa                                                              >  1:2273445/1‑150 (MQ=255)
                                                                                             cacaGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGc                                                          >  1:1603888/1‑150 (MQ=255)
                                                                                              acaGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCg                                                         <  2:386342/150‑1 (MQ=255)
                                                                                                aGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGt                                                       >  2:687288/1‑150 (MQ=255)
                                                                                                 gCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTg                                                      <  1:687288/150‑1 (MQ=255)
                                                                                                  cccTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCGGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTgg                                                     <  2:215711/150‑1 (MQ=255)
                                                                                                     tACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGaaa                                                  >  2:608550/1‑150 (MQ=255)
                                                                                                            cATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCACGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCt                                           <  2:2151436/150‑1 (MQ=255)
                                                                                                             aTTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCgg                                                                          >  1:2158640/1‑118 (MQ=255)
                                                                                                             aTTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCgg                                                                          <  2:2158640/118‑1 (MQ=255)
                                                                                                             aTTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGaaa                                                  <  1:922579/142‑1 (MQ=255)
                                                                                                             aTTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGaaa                                                  >  2:922579/1‑142 (MQ=255)
                                                                                                             aTTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTg                                          >  2:2317557/1‑150 (MQ=255)
                                                                                                             aTTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTg                                          >  2:1933437/1‑150 (MQ=255)
                                                                                                             aTTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAACATCCAGTCCGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGCGCAGAGGCCGGTGGAAACCATCCTg                                          >  1:41050/1‑150 (MQ=255)
                                                                                                                   ctgctACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAg                                                                                                         >  1:894970/1‑81 (MQ=255)
                                                                                                                   ctgctACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGAc                                    >  1:1451592/1‑150 (MQ=255)
                                                                                                                        aCTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTg                                                      >  1:335030/1‑127 (MQ=255)
                                                                                                                           cACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGacca                            <  2:860088/150‑1 (MQ=255)
                                                                                                                           cACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGAAGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGCCCCGCTGCCGGCCCAGGTGCAGATCGCGGTGGGAACCATCCTGCACGACATCCaccc                            >  2:1189730/1‑149 (MQ=255)
                                                                                                                            aCCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGc                                                                         >  2:452725/1‑104 (MQ=255)
                                                                                                                            aCCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGaccac                           <  2:1451592/150‑1 (MQ=255)
                                                                                                                            aCCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGaccac                           >  2:1607619/1‑150 (MQ=255)
                                                                                                                             ccGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGaccacc                          >  2:314208/1‑150 (MQ=255)
                                                                                                                                     aaCCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCg                  >  1:1661603/1‑150 (MQ=255)
                                                                                                                                      aCCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGa                 >  2:1470041/1‑150 (MQ=255)
                                                                                                                                      aCCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGa                 >  1:2210577/1‑150 (MQ=255)
                                                                                                                                         aaGGCCCCACCCTGGAGCGCCTGGGTGGCCCGACCATGGCCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCCGCTCGCGGACGACACCACCCCGCACGCCATCGCCCCCCACTACTCCACCt              >  2:1774973/1‑150 (MQ=255)
                                                                                                                                           ggCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTAc            <  1:757411/150‑1 (MQ=255)
                                                                                                                                            ggaccACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTAcc           >  1:397699/4‑150 (MQ=255)
                                                                                                                                             gaccACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTc                   <  1:1637203/139‑1 (MQ=255)
                                                                                                                                              cccACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATc                                             >  1:882940/1‑114 (MQ=255)
                                                                                                                                              cccACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATc                                             <  2:882940/114‑1 (MQ=255)
                                                                                                                                              cccACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACcgc         >  2:927643/1‑150 (MQ=255)
                                                                                                                                               ccACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACcgcg        >  2:5984/1‑150 (MQ=255)
                                                                                                                                                cACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCgcgg       >  2:2170727/1‑150 (MQ=255)
                                                                                                                                                cACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCgcgg       <  2:397699/150‑1 (MQ=255)
                                                                                                                                                   ccTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATc                                                                                                            >  2:636275/1‑46 (MQ=255)
                                                                                                                                                   ccTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATc                                                                                                            <  1:636275/46‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATT  >  NC_002947/2456885‑2457183

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