Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,458,832 C→A 9.2% intergenic (+118/‑63) 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,8320CA 9.2% 278.7 / 17.4 106intergenic (+118/‑63)gapB/PP_2150glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (47/47);  new base A (6/4);  total (53/53)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01

Only 100 of 109 total aligned reads displayed.

CTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATG  >  NC_002947/2458683‑2458981
                                                                                                                                                      |                                                                                                                                                     
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)
        aaCCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTg                                                                                                                                               <  1:2828765/150‑1 (MQ=255)
          cccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTa                                                                                                                                             <  1:1442288/150‑1 (MQ=255)
          cccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTa                                                                                                                                             <  2:3117626/150‑1 (MQ=255)
          cccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTa                                                                                                                                             <  1:3676477/150‑1 (MQ=255)
          cccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTa                                                                                                                                             <  2:4565558/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)
             gCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGGTTATCGCTGCCCCAGTAGGa                                                                                                                                          >  2:251077/1‑150 (MQ=255)
                 gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                      <  2:1424029/150‑1 (MQ=255)
                 gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                      <  2:174100/150‑1 (MQ=255)
                 gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                      <  2:215979/150‑1 (MQ=255)
                  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)
                                                                         ccGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTATCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCCAGAAACCCTACAGTGCGCACAGc                                                                              >  2:3820725/1‑150 (MQ=255)
                                                                           gTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTACTACAGGATTATCGCAGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCAAATTTGCGCACAcccc                                                                            >  2:4410635/1‑150 (MQ=38)
                                                                              ttttttttttTATTATTTTGTTTTGTCTGTGCGGGCCCTTTGGGGGGAAAACCCGCTCCTACAGGTTTTTCGGTGCCCGGGTTGGTGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTa                                                                          <  1:1971126/150‑1 (MQ=255)
                                                                               ttttttttttCTGATTGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGTTTTATCGCTCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTAc                                                                         <  1:2835905/150‑1 (MQ=38)
                                                                                ttttttATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACt                                                                        <  2:4594105/150‑1 (MQ=255)
                                                                                 tttttATTCTGATTGTGTGTTTTCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTATCGTATTATCGCTGCCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTg                                                                       <  1:5123258/150‑1 (MQ=38)
                                                                                  ttttATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCCCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCCAGAAACCCTACAGTGCGCACACCCCAACAGt                                                                      >  1:4652693/1‑150 (MQ=38)
                                                                                        tCTGATTGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGAAAACCCGCTCCTACCGTTTTATCGCTGCCCCTGTAGGAGGGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTc                                                                <  1:3228878/150‑1 (MQ=39)
                                                                                          tGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGATTATCGCTGCCCCAGTAGGAGCGGGTTTACCCGCGAAGAGGCCACCACCGGCAACCGAGAAACCCGACTGTGCGCACAGCCCTACAGTTCATTatg                                                              >  2:3012536/1‑150 (MQ=38)
                                                                                           gAATGTGTGTTGTCTGTGCCGGCCCTTGCGCGGGTAAACCCGCTCATACAGGATTATCGCTGCCCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTctgc                                                             <  2:1700386/150‑1 (MQ=255)
                                                                                              tgtgttTTCACAGTCCCGGCCCATCCCCGGTAAAACCCCCTCCAACAGCATTAACCCTCCCCCTCTACGAGCGGTTTAACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCtgctgc                                                          <  1:2702380/150‑1 (MQ=255)
                                                                                              tgtgtgTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGCTTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACCGGCAACCGAGAAACCCGACAGTGCGCACACCCCAACTGGTCATTCtgctgc                                                          >  1:4235339/1‑150 (MQ=38)
                                                                                              tgtgtgTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCCACCGAGAACCCCTACGGTGCGCACAGCCCTACTGTTCATTGtgctgc                                                          >  1:4235400/1‑150 (MQ=255)
                                                                                              ggtgtgtTGTCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGTATTATCGCTGCCCCTGTTGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCtgctgc                                                          <  1:3745283/149‑1 (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)
                                                                                                                                  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:3998943/1‑150 (MQ=255)
                                                                                                                                            cggtttATCGCTGCCCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtgg            <  2:1303715/146‑1 (MQ=39)
                                                                                                                                            aGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCACCCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCAGTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtgg            >  1:660654/1‑150 (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)
                                                                                                                                                     gCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCAt   >  2:4952833/1‑150 (MQ=255)
                                                                                                                                                      cTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAcacag  >  2:1978726/1‑146 (MQ=255)
                                                                                                                                                      cTGCCCCTGTAGGAGCAGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATg  <  1:251077/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATG  >  NC_002947/2458683‑2458981

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: