Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,458,827 T→C 9.5% intergenic (+113/‑68) gapB → / → PP_2150 glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,458,8270TC 9.5% 299.5 / 16.2 105intergenic (+113/‑68)gapB/PP_2150glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (47/48);  new base C (6/4);  total (53/52)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

Only 100 of 108 total aligned reads displayed.

TGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGA  >  NC_002947/2458680‑2458976
                                                                                                                                                    |                                                                                                                                                     
tGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTAt                                                                                                                                                     <  1:4859620/150‑1 (MQ=255)
 ggCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGGGAACCCCCACCGACCGGAGTAtg                                                                                                                                                    >  1:1303715/1‑149 (MQ=255)
   cTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGc                                                                                                                                                  <  2:4104099/150‑1 (MQ=255)
    tGGTGTGAGCCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCt                                                                                                                                                 <  1:4592365/150‑1 (MQ=255)
     ggTGTGAGCCCGCCAGCGTTTCCACACTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCt                                                                                                                                                 <  2:5098524/150‑1 (MQ=255)
     ggTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGGAAGGCAATGAAAAAACGGGAACCCTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTg                                                                                                                                                <  2:4373812/150‑1 (MQ=255)
      gtgtGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGc                                                                                                                                               >  1:4519635/1‑150 (MQ=255)
           aaCCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTATAACCGCTCCTACAGGATTATCGCTGCCCCTg                                                                                                                                          >  2:2126491/1‑150 (MQ=255)
           aaCCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTg                                                                                                                                          <  1:2828765/150‑1 (MQ=255)
             cccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTa                                                                                                                                        <  1:3676477/150‑1 (MQ=255)
             cccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTa                                                                                                                                        <  2:3117626/150‑1 (MQ=255)
             cccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTa                                                                                                                                        <  2:4565558/150‑1 (MQ=255)
             cccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTa                                                                                                                                        <  1:1442288/150‑1 (MQ=255)
              ccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAg                                                                                                                                       <  2:1949759/150‑1 (MQ=255)
              ccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCGACAGGATTATCGCTGCCCCTGTAg                                                                                                                                       >  2:770752/1‑150 (MQ=255)
                gCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCTGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGa                                                                                                                                     <  1:3570742/150‑1 (MQ=255)
                gCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGa                                                                                                                                     <  2:3549130/150‑1 (MQ=255)
                gCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGa                                                                                                                                     <  1:3570716/150‑1 (MQ=255)
                    gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                 <  2:215979/150‑1 (MQ=255)
                    gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                 <  2:1424029/150‑1 (MQ=255)
                    gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                 <  2:174100/150‑1 (MQ=255)
                     cGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCggg                                                                                                                                <  1:562923/150‑1 (MQ=37)
                     cGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCggg                                                                                                                                <  1:4906362/150‑1 (MQ=37)
                     cGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCggg                                                                                                                                <  2:2998549/150‑1 (MQ=255)
                      gTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTATTTTTTTATTCTGAATGTGTGTTGTCGGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCGGGt                                                                                                                               <  2:1473730/150‑1 (MQ=18)
                       tttCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgggat                                                                                                                              <  1:95862/150‑3 (MQ=255)
                       tttCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCGGGtt                                                                                                                              <  1:601788/150‑1 (MQ=25)
                            aCGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTGAAACCGCACCTACCGGAGTATGGCTGACCCAGGAGGGGCGGGTTTTcc                                                                                                                          >  1:2198102/1‑150 (MQ=37)
                              gCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTTAAACCGCTCCTACAGGATTATAGCTGCCCCTGTAGGAGCGGGTTTACCcgc                                                                                                                       >  2:3970283/1‑150 (MQ=255)
                                 gtttCGCTTGTATGGAAATTAAAAAACGGGAACCTTGGGTTTCCCTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGGGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                    <  1:2002228/148‑1 (MQ=38)
                                 gatgcccTTGTAAGGCAATGAAAAAACGGGAACCGCCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTGGCGCGGGTAAACCCGCTCCCCCAGGATTATCGCGGCCCCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                    <  2:4296996/144‑1 (MQ=38)
                                  aTTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTGGGGTTCCTGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACCGTATTATCGCTGCCCCTTTTGGAGCGGTTTTACCCGCGAag                                                                                                                   <  1:525832/150‑1 (MQ=38)
                                  aTTCGCTTGTAAGGCAATGAAAAAACGGCAACATTCGGGTTCCTGTTTTTTTTTATACTGAATGTGTGTTGTCAGTGCCGGCCCTTACGCGGGTAAACCCGATCCTACAGGATTATCGCTCCCCCTGTAGGAGCGGGTTTACCCGCGAag                                                                                                                   <  1:1910486/150‑1 (MQ=38)
                                   gtCGCTTGTACGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAaga                                                                                                                  <  2:5038176/149‑1 (MQ=255)
                                    tcgcaagTAAGGCAATGAAAAAACGGGAACCATCGCGATCCCCTTTTTTTATATTCGGAATGTGTGGTGACAGAGCCGGCCCTTTCGCGGGTAAAACCGCTCCGACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                                 >  2:5078960/7‑150 (MQ=37)
                                    tCGCTTGTACGGCAATTAAAAAAGGGTAACCTTCGGGTTCCCTTTTTTTTTTATTCTGACTGTGTGTTGTCAGTGCGGGCCCTTTCGCGGTTAAACCCGCTCCTACCGTTTTATCGCTGCCCCTGTAGGAGGGGTTTTACCCGCGAagag                                                                                                                 <  2:1501884/150‑1 (MQ=255)
                                    tCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                                 >  2:5078986/1‑150 (MQ=255)
                                        ttGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCGACAGGCATATCGCTGCCCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCa                                                                                                             >  2:1365054/1‑150 (MQ=255)
                                             aGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGGTTACCCGCGAAGAGGCCCGCAcc                                                                                                        >  2:822934/1‑150 (MQ=255)
                                             aGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGGTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCAACAGGATTATCGCTGCCCCGGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGCAcc                                                                                                        >  2:2230204/1‑150 (MQ=38)
                                               gCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCCCACCAACAGGATTAAAGCTGCCCCAGTAGGAGCGGGTTTAACCGCGAAGAGGCCCGCACCgg                                                                                                      >  1:4902527/1‑150 (MQ=255)
                                                        ataCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGAGGGTAAACACGAGCCTAAAGGATTATCGCTGACCCTGTAGGAGCGGGTTTAACCGCGAAGAGGCCAGCACCGGCAACCgaga                                                                                             >  2:3639184/3‑150 (MQ=38)
                                                        aaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCGACAGGCTTATAGCGGCCCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCcaga                                                                                             >  2:1346952/1‑150 (MQ=38)
                                                          aCGGGACCCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCGCCGCCAGGCGTATCGCGGCCCCGGGAGGAGCCGGGTTACCCGCGCCGAGGCCCGCACCGGCACCCgcgcca                                                                                           >  2:4376351/1‑145 (MQ=255)
                                                           cGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAAc                                                                                          <  2:3380872/150‑1 (MQ=255)
                                                           cGGGAACCTTCGGGTGCCCTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTAACGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAAc                                                                                          <  2:228353/150‑1 (MQ=255)
                                                            gggAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCCCTCCTACAGGAGTAACGCAGCCCCAGTAGGCGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCCACCCAGAAAcc                                                                                         >  2:892300/1‑150 (MQ=38)
                                                               aaCCTTCGGGTTCCCGTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAAccca                                                                                       >  2:3745283/1‑149 (MQ=255)
                                                                   ttggggTTCCTTTTTTTTTTTTTTTTTATTGTGTGTTGCCTGGGCCGGCCGTTTCGGGGGTAACCCCGCTCCTACCTTTTTACGGGTGCCCGTGTCGGGGGGGTTTTCCCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACtgt                                                                                  <  1:4345695/147‑1 (MQ=255)
                                                                   ttCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACagt                                                                                  >  2:1927645/1‑150 (MQ=255)
                                                                   ttCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCACCAACACGATTAAAGCTGCCCCAGTATGAGCCGGTTTAACCGCGAACAGGCCACCACCCGCAACCCAGAAACCCAACtgt                                                                                  >  2:1927636/1‑150 (MQ=255)
                                                                       ggTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGAATAACGCTGCCCCAGTAGGAGCCGGTTTACACGCCAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTgcgc                                                                              >  2:281590/1‑150 (MQ=39)
                                                                              gTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTACTACAGGATTATCGCAGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCAAATTTGCGCACAcccc                                                                       >  2:4410635/1‑150 (MQ=38)
                                                                                 ttttttttttTATTATTTTGTTTTGTCTGTGCGGGCCCTTTGGGGGGAAAACCCGCTCCTACAGGTTTTTCGGTGCCCGGGTTGGTGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTa                                                                     <  1:1971126/150‑1 (MQ=255)
                                                                                  ttttttttttCTGATTGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGTTTTATCGCTCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTAc                                                                    <  1:2835905/150‑1 (MQ=38)
                                                                                    tttttATTCTGATTGTGTGTTTTCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTATCGTATTATCGCTGCCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTg                                                                  <  1:5123258/150‑1 (MQ=38)
                                                                                     ttttATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCCCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCCAGAAACCCTACAGTGCGCACACCCCAACAGt                                                                 >  1:4652693/1‑150 (MQ=38)
                                                                                          ttCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTTGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCAtt                                                            <  2:4802567/150‑1 (MQ=255)
                                                                                           tCTGATTGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGAAAACCCGCTCCTACCGTTTTATCGCTGCCCCTGTAGGAGGGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTc                                                           <  1:3228878/150‑1 (MQ=39)
                                                                                             tGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGATTATCGCTGCCCCAGTAGGAGCGGGTTTACCCGCGAAGAGGCCACCACCGGCAACCGAGAAACCCGACTGTGCGCACAGCCCTACAGTTCATTatg                                                         >  2:3012536/1‑150 (MQ=38)
                                                                                                 tgtgttTTCACAGTCCCGGCCCATCCCCGGTAAAACCCCCTCCAACAGCATTAACCCTCCCCCTCTACGAGCGGTTTAACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCtgctgc                                                     <  1:2702380/150‑1 (MQ=255)
                                                                                                 tgtgtgTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGCTTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACCGGCAACCGAGAAACCCGACAGTGCGCACACCCCAACTGGTCATTCtgctgc                                                     >  1:4235339/1‑150 (MQ=38)
                                                                                                 ggtgtgtTGTCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGTATTATCGCTGCCCCTGTTGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCtgctgc                                                     <  1:3745283/149‑1 (MQ=255)
                                                                                                   tgtgtTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCGACAGGATTATCGCTGCCCCGGTTGGAGCGGGTTTACCCGCGAGGAGGCCCGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTa                                                   >  2:3062449/1‑150 (MQ=255)
                                                                                                     tgttgtCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTCGGAGCGGGTTTACCCGCGAAGAGGCCCGCACCCGCAACCGAGAAACCCTACGGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACt                                                 >  2:4768540/1‑150 (MQ=255)
                                                                                                      gttgtCTGTGCCGGCCTTTTGGGGGGTAAACCCGCTCCTACTGTTTTTGGGCTCCCCGTTTTGGTGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTc                                                <  1:684326/150‑1 (MQ=255)
                                                                                                      gttgcCAGTGCCGGCCCTTTCGCGGTAAACCCCGCTCCAACCGCTTTATCGCGGCCCCTGTTGGTGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTc                                                <  1:684230/150‑1 (MQ=255)
                                                                                                         gtCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCCACCCAGAAACCCTCCAGTGCGCACAGCCCACCAGTTCATTCTGCTGCTAAACAcc                                             >  2:3846858/1‑150 (MQ=37)
                                                                                                                 cGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGATTATCGCTGCCCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCACCCGCAACCGAGAAACCCGACTGTGCGCACAGCCCTACGGTTCATTCTGCTGCTACTCACCCCCAGCaa                                     >  2:2237657/1‑150 (MQ=38)
                                                                                                                 cGGCCCTTTCGCGGGTAAACCCCCGCCGACAGGATTATAGCTGCCCCTGTTGGAGCGGGGTTACCCCCCAAGAGGCCAGCACCGGCCACCGCGAAACCCAACGCTGCGCCCAGCCCTACTGTTCATTCTGCAGCTACTCACCGCCAGCaa                                     >  2:2237448/1‑150 (MQ=255)
                                                                                                                  ggCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCAACAGTTCATTCTGCTGCTACTCACCGCCTGCAAc                                    >  2:2477067/1‑150 (MQ=255)
                                                                                                                  ggCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTAATTCTGCTGCTACTCACCGCCTGCAAc                                    >  2:3389305/1‑150 (MQ=255)
                                                                                                                     ccTTTCGCGGGTAAACCCGCTCCTACAGTTTTATCGCGGCCCCTTTTGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCaa                                 <  1:619954/150‑1 (MQ=39)
                                                                                                                        ttCGCGGGTAAACCCGCTCCTACAGGAATATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAACCGAGAAACCCTACTGTGCGCACAGCCCAAAAGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGc                              >  2:4204411/1‑150 (MQ=38)
                                                                                                                                  aaCCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTgg                    <  1:1998887/150‑1 (MQ=255)
                                                                                                                                   aCCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGACAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGa                   >  1:1550590/1‑150 (MQ=18)
                                                                                                                                   aCCCGCCCATACTGTTTTCTCGCTGCCCCTGTTGGTGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGa                   <  1:770752/150‑1 (MQ=255)
                                                                                                                                    cccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAggcg                  >  2:4013980/1‑148 (MQ=255)
                                                                                                                                     ccGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAACCCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgc                 >  2:3674791/1‑150 (MQ=37)
                                                                                                                                     ccGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgc                 >  2:2563303/1‑150 (MQ=37)
                                                                                                                                     ccGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgc                 >  2:182468/1‑150 (MQ=25)
                                                                                                                                     ccGCTCCTACAGGGGCAGCGATAATCCTGTACGAGCGGGTTTACCCGCGAAAGGGCCAGCACCGGCAACCCAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATACTGCTGCTACTCACCGCCTCCAACCAAGGCCCCACCCTGGAgc                 >  2:182483/1‑150 (MQ=21)
                                                                                                                                     ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgc                 >  1:2433071/1‑150 (MQ=255)
                                                                                                                                     ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCAGCAACCAAGGCCCCACCCTGGAgc                 >  2:4504271/1‑150 (MQ=255)
                                                                                                                                      cGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTctg                                                         >  2:872027/1‑109 (MQ=255)
                                                                                                                                      cGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTctg                                                         <  1:872027/109‑1 (MQ=255)
                                                                                                                                       gCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgcgc               >  1:1505636/1‑150 (MQ=255)
                                                                                                                                       gCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgcgc               >  1:1505597/1‑150 (MQ=255)
                                                                                                                                         tCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCAGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCt             >  2:671551/1‑150 (MQ=255)
                                                                                                                                           cTACAGGTTTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTgg           <  2:4902527/150‑1 (MQ=255)
                                                                                                                                             aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggt         >  2:1419659/1‑150 (MQ=255)
                                                                                                                                             aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggt         >  2:3998943/1‑150 (MQ=255)
                                                                                                                                               cggtttATCGCTGCCCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtgg       <  2:1303715/146‑1 (MQ=39)
                                                                                                                                               aGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtgg       <  1:4013980/150‑1 (MQ=255)
                                                                                                                                                ggATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCGACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtggg      <  2:4519635/150‑2 (MQ=39)
                                                                                                                                                 gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTc     >  1:2162172/1‑150 (MQ=255)
                                                                                                                                                 gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTc     >  1:2162167/1‑150 (MQ=255)
                                                                                                                                                  tttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTcc    <  1:4560953/149‑1 (MQ=255)
                                                                                                                                                  aTTATCGCTGCCCCTGTTGGAGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTcc    <  1:1052839/150‑1 (MQ=255)
                                                                                                                                                    tATCGCTGCCCCTGTTGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGa  <  1:3970283/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGA  >  NC_002947/2458680‑2458976

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