Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,457,037 A→C 5.1% E22A (GAG→GCG)  PP_2150 → thiamine biosynthesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,457,0370AC 5.1% 339.3 / 4.4 119E22A (GAG→GCG) PP_2150thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (61/50);  new base C (4/3);  total (66/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

ACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATT  >  NC_002947/2456888‑2457183
                                                                                                                                                     |                                                                                                                                                  
  ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgc                                                                                                                                                  >  2:1780796/1‑150 (MQ=255)
  ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgc                                                                                                                                                  <  1:779343/150‑1 (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)
            aGGATTATCGCTGCCCCTGTGGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtgg                                                                                                                                        >  2:1022439/1‑150 (MQ=255)
                ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTgg                                                                                                                                            >  1:1070543/1‑142 (MQ=255)
                ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCg                                                                                                                                    >  1:386342/1‑150 (MQ=255)
                   tCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAcc                                                                                                                                 >  1:101827/1‑150 (MQ=255)
                      cTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATg                                                                                                                              <  1:580204/150‑1 (MQ=255)
                       tGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATgg                                                                                                                             <  1:1418541/150‑1 (MQ=255)
                         ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCAACGGCAACCAATAAACCCTACTGTGCGCACAGCCATACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgc                                                                                                                           <  1:2271147/150‑1 (MQ=255)
                          cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCAt                                                                                                                               >  2:904514/1‑145 (MQ=255)
                          cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCAt                                                                                                                               <  1:904514/145‑1 (MQ=255)
                          cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCAACACACCTCCTGCACCCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgca                                                                                                                          >  2:1716353/1‑150 (MQ=255)
                            cTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgcagc                                                                                                                        <  1:1320262/150‑1 (MQ=255)
                             tGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg                                                                                                                            <  1:2022116/145‑1 (MQ=255)
                             tGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg                                                                                                                            >  2:2022116/1‑145 (MQ=255)
                             tGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgcagca                                                                                                                       <  1:713017/150‑1 (MQ=255)
                               tAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCagcagc                                                                                                                     >  1:251282/1‑150 (MQ=255)
                                     cGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGc                                                                                                               <  2:1684816/150‑1 (MQ=255)
                                     cGGGGTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGCTCATTCTGCTGCCACTCACCGCCTGCAACCAACGCCCCACCCTGGCGCGCCTGGGGGGTCCCACCATGGCCACCAGCGACAGc                                                                                                               >  2:716792/1‑150 (MQ=255)
                                          ttACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCa                                                                                                          <  2:697951/150‑1 (MQ=255)
                                             cccGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCGACTTTGCGCACATCCCTACTGTTCATTCTGCTGCGACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACGGCATCCAGTa                                                                                                       >  1:1332624/1‑150 (MQ=255)
                                              ccGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCAt                                                                                                                               >  1:2114159/1‑125 (MQ=255)
                                              ccGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCAt                                                                                                                               <  2:2114159/125‑1 (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                                                                                                                   <  1:1434311/127‑1 (MQ=255)
                                                        gCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTa                                                                                                                   >  2:1434311/1‑127 (MQ=255)
                                                        gCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTcc                                                                                                  <  2:258414/144‑1 (MQ=255)
                                                        gCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTcc                                                                                                  >  1:258414/1‑144 (MQ=255)
                                                          cAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATc                                                                                                            <  2:1546948/132‑1 (MQ=255)
                                                          cAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATc                                                                                                            >  1:1546948/1‑132 (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)
                                                                cGGCAACCGGGAAACCCTACTGTGCGCACAGCGCTACTGTTCATCCTGCTGCTACTCACCGCCTGCAGCGAAGGCCCCACCCTGGGGCGCCTGGGGGGTACGAGCATGGACTTGTGCTACAGCAGGCGGTAGGTGCGCGAGCCtggtgac                                                                                    >  2:491976/1‑148 (MQ=255)
                                                                 ggCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGcc                                                                                   <  1:1997936/150‑1 (MQ=255)
                                                                  gCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGccc                                                                                  >  2:1980688/1‑150 (MQ=255)
                                                                       cGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCACGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGc                                                                             >  2:2311109/1‑150 (MQ=255)
                                                                       cGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAg                                                                                                         >  2:675498/1‑122 (MQ=255)
                                                                       cGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAg                                                                                                         <  1:675498/122‑1 (MQ=255)
                                                                       cGAGAAACCCTACTGTGCGCACAGCCCTACGTTTCATTCTGCTGCTACGCACCGCCTGCAACCAAGGCCCCACCCTTGATCGCCTGGGTCGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCTCGCTGc                                                                             <  2:620381/150‑1 (MQ=255)
                                                                        gagaAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTCAATCATCCCGTACGTCCGCGAGCCTGGTGGCCCGCTGcc                                                                            >  1:123977/1‑150 (MQ=255)
                                                                               ccTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCcgc                                                                                                <  1:1040609/123‑1 (MQ=255)
                                                                               ccTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCcgc                                                                                                >  2:1040609/1‑123 (MQ=255)
                                                                               ccTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAg                                                                     <  2:1957132/150‑1 (MQ=255)
                                                                               ccTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAg                                                                     <  2:2329793/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:101827/150‑1 (MQ=255)
                                                                                             aGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGt                                                       >  2:687288/1‑150 (MQ=255)
                                                                                              gCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTg                                                      <  1:687288/150‑1 (MQ=255)
                                                                                                 cTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGAAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGaa                                                   <  1:1478782/150‑1 (MQ=255)
                                                                                                       ttCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATc                                             <  2:1758195/150‑1 (MQ=255)
                                                                                                         cATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCACGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCt                                           <  2:2151436/150‑1 (MQ=255)
                                                                                                          aTTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCgg                                                                          <  2:2158640/118‑1 (MQ=255)
                                                                                                          aTTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCgg                                                                          >  1:2158640/1‑118 (MQ=255)
                                                                                                          aTTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGaaa                                                  <  1:922579/142‑1 (MQ=255)
                                                                                                          aTTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGaaa                                                  >  2:922579/1‑142 (MQ=255)
                                                                                                          aTTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTg                                          >  2:1933437/1‑150 (MQ=255)
                                                                                                          aTTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTg                                          >  2:2317557/1‑150 (MQ=255)
                                                                                                          aTTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAACATCCAGTCCGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGCGCAGAGGCCGGTGGAAACCATCCTg                                          >  1:41050/1‑150 (MQ=255)
                                                                                                                ctgctACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAg                                                                                                         <  2:894970/81‑1 (MQ=255)
                                                                                                                     aCTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTg                                                      <  2:335030/127‑1 (MQ=255)
                                                                                                                        cACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGAAGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGCCCCGCTGCCGGCCCAGGTGCAGATCGCGGTGGGAACCATCCTGCACGACATCCaccc                            >  2:1189730/1‑149 (MQ=255)
                                                                                                                         aCCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGc                                                                         >  2:452725/1‑104 (MQ=255)
                                                                                                                         aCCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGc                                                                         <  1:452725/104‑1 (MQ=255)
                                                                                                                         aCCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGaccac                           >  2:1607619/1‑150 (MQ=255)
                                                                                                                         aCCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGaccac                           <  2:1451592/150‑1 (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)
                                                                                                                                    ccAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACAACTACTCgtc                <  2:807056/150‑3 (MQ=255)
                                                                                                                                      aaGGCCCCACCCTGGAGCGCCTGGGTGGCCCGACCATGGCCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCCGCTCGCGGACGACACCACCCCGCACGCCATCGCCCCCCACTACTCCACCt              >  2:1774973/1‑150 (MQ=255)
                                                                                                                                         ggaccACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTAcc           >  1:397699/4‑150 (MQ=255)
                                                                                                                                          gaccACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTc                   >  2:1637203/3‑141 (MQ=255)
                                                                                                                                          gaccACCCAGGCGCTCCAGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTc                   <  1:1637203/139‑1 (MQ=255)
                                                                                                                                           cccACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATc                                             <  2:882940/114‑1 (MQ=255)
                                                                                                                                           cccACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATc                                             >  1:882940/1‑114 (MQ=255)
                                                                                                                                           cccACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACcgc         >  2:927643/1‑150 (MQ=255)
                                                                                                                                            ccACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACcgcg        >  2:1454937/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                                                                                                            <  1:636275/46‑1 (MQ=255)
                                                                                                                                                ccTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATc                                                                                                            >  2:636275/1‑46 (MQ=255)
                                                                                                                                                  tGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCAAGACATCGACCACCACTACTCGACCTACCGCGGCGAtt  >  1:1956600/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGCAGCAGCTACAGCATCCAGTACGTCCGCGAGCCTGGTGGCCCGCTGCCGGCCCAGGTGCAGATGGCGGTGGAAACCATCCTGCACGACATCGACCACCACTACTCGACCTACCGCGGCGATT  >  NC_002947/2456888‑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: ‑