Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,461,862 C→A 5.4% intergenic (+121/‑60) 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,461,8620CA 5.4% 678.9 / 18.6 211intergenic (+121/‑60)gapB/PP_2150glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (99/99);  new base A (2/10);  total (102/109)
Fisher's exact test for biased strand distribution p-value = 3.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 221 total aligned reads displayed.

GTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGG  >  NC_002947/2461713‑2462010
                                                                                                                                                         |                                                                                                                                                    
  gtgAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGccc                                                                                                                                                    <  2:862074/150‑1 (MQ=255)
  gtgAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGccc                                                                                                                                                    <  1:662453/150‑1 (MQ=255)
  gtgAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGccc                                                                                                                                                    <  1:840892/150‑1 (MQ=255)
      aCCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGATTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCt                                                                                                                                                  <  2:8385975/150‑1 (MQ=255)
      aCCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGt                                                                                                                                                <  1:1025613/150‑1 (MQ=255)
      aCCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGt                                                                                                                                                <  1:30311/150‑1 (MQ=255)
      aCCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGAGCCGGCCCTTTCGCGGGTTAACCCGCGCCAACACGCTTATAGCTGccccagg                                                                                                                                                >  1:455375/1‑147 (MQ=255)
       cccGCCAGCGTTTCCACGCTGATTCGCTTTTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCACG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTa                                                                                                                                               <  1:7234646/150‑1 (MQ=255)
       cccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTa                                                                                                                                               <  1:6904843/150‑1 (MQ=255)
       cccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCGCCAACACGAATAAAGCGGCCCCGGTa                                                                                                                                               >  1:1774196/1‑150 (MQ=255)
        ccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAg                                                                                                                                              <  1:4649885/150‑1 (MQ=255)
        ccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAg                                                                                                                                              <  1:7388541/150‑1 (MQ=255)
           ccAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCGGTAGGAg                                                                                                                                           >  2:3431026/1‑150 (MQ=255)
             aGCGTTTCCACGCTTATTCGCTTTTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTGTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGc                                                                                                                                          <  2:2733915/150‑1 (MQ=255)
             aGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCGGTAGGAGCg                                                                                                                                         >  2:1185346/1‑150 (MQ=255)
              gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGGATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                        <  1:4762639/150‑1 (MQ=255)
              gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                        <  2:754047/150‑1 (MQ=255)
              gCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCGGTAGGAGCgg                                                                                                                                        >  2:2640390/1‑150 (MQ=255)
               cGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCgg                                                                                                                                        <  1:6328386/150‑1 (MQ=25)
               cGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCggg                                                                                                                                       <  2:499934/150‑1 (MQ=255)
                     cACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCGGGTTTAcc                                                                                                                                 <  2:7686959/150‑1 (MQ=21)
                      aCGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTAccc                                                                                                                                <  1:5732737/150‑1 (MQ=255)
                      aCGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCACCAACAGGATTATAGCTGCCCCAGTAGGAGCCGGTTTAacc                                                                                                                                >  1:1238399/1‑150 (MQ=39)
                        gctgctGCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCcgc                                                                                                                              <  1:6477947/150‑1 (MQ=255)
                           gATTCGCTTGTACGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                           <  2:1991929/150‑1 (MQ=255)
                            aTTCGCTTGTAAGGAAATGAAAAAACGGTAACCTTCGGTTTCCCT‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCGGGTTTACCCGCGaaa                                                                                                                          <  2:1611594/150‑2 (MQ=18)
                             tgcGCTTGTAAGGCATTGAAAAAACGGGAACCTTCGGGTACCCG‑‑TTTTTTTTTATTCTGACTGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTTGGTGCGGGTTTACCCGCGAaga                                                                                                                         <  2:6576902/148‑1 (MQ=39)
                                gCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCAGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                                      >  2:5234524/1‑150 (MQ=255)
                                   tGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTGAAACCCCTACTAAAGGATTATCGCTGGCCCGGGTGGAGCGGGGTTAACCGCGAAGTgggcc                                                                                                                    >  1:7721360/1‑147 (MQ=37)
                                     tAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCa                                                                                                                 >  2:8542981/1‑150 (MQ=255)
                                     tAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCGACAGGGTTATCGCTGCCCCTGGAGGAGCGGGTTTACCCGCGCACGGGCCAGCa                                                                                                                 >  2:5526407/1‑150 (MQ=38)
                                      acggCAATGAAAAAACGGCCCCCTCCGGGTACCCG‑‑ATTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTGCGCGGGTAAACCCGCCCCGACAGGATTATCGCTGCCCCTGCAGGAGCGGGTTTACCCGCGAAGAGGCCAGCAc                                                                                                                <  1:6421720/148‑1 (MQ=38)
                                       aGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCCCTCCGACAGGCTTATCGCTGCCCCGGTAGGAGCGGGTTTTCCCGCGCAGAGGCCAGCAcc                                                                                                               >  1:2036955/1‑150 (MQ=39)
                                         gCAATTAAAAAACGGGAACCTTCGGGTTCCCTT‑‑TTTTTTTTTATTCTTATTGTTTGTTGTCTGTGCCGGCCCTTTCGCGGTTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCg                                                                                                              <  2:332793/150‑1 (MQ=37)
                                           aaTGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCGCCTACAGGATTATCGCTGCCCCCGTAGGGGCGGGTTTACCCGCGCAGAGGCCCGCACCGGCa                                                                                                           >  2:440981/1‑150 (MQ=39)
                                            aaaaaaaaaCGGGAACCCGAAGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTTAAACCGCGCCGACAGGGGTATAGCTGCCCCGGGAGGGGGGGGGTTTCCCCCGGAGCGGCCCGCCCCgggca                                                                                                          >  1:1116782/4‑147 (MQ=255)
                                            aTGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCCCACCTAAAGGATTATCGCGGCCCCAGTAGGAGCGGGTTTAACCGCGCAGAGGCCAGCACCCGcca                                                                                                          >  2:3174650/1‑148 (MQ=255)
                                             aaaaaaaaCGGGAACCCGAAGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCCCGCCTACAGGATTATCGCTTCCCCGGTAGGAGCGGGGTTAACCCCGAAGAGGCCAGCACCCGCAAc                                                                                                         >  1:6834829/3‑150 (MQ=255)
                                              aaaaaaaCGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTAAAGGCTTATCGCTGCCCCGGTAGGAGCGGGTGTACCCGCGAAGAGGCCAGCACCGGCAAcc                                                                                                        >  1:8014381/2‑150 (MQ=39)
                                                aaaaaCGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCga                                                                                                      >  1:2895297/1‑150 (MQ=255)
                                                  aaaCGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACCGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCgaga                                                                                                    <  2:654074/150‑1 (MQ=255)
                                                    aCGGGAACCCGAAGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGATTATCGCTGCCCCGGTAGGAGAGGGTTTAACCCCGAAGAGGCCAGCACCGGCAAccccgaac                                                                                                  >  2:2929919/1‑144 (MQ=38)
                                                    aCGGCAACCTACGGGTACCCT‑‑TATTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCCCCTACAGGATTAACGCTCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGaaa                                                                                                  <  1:8776027/150‑1 (MQ=38)
                                                      gggACCCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCGACGACAGGATTATCGCTGCCCCGGTCGGCGCGGGGTTACCCGCGAAGAGGCCCGCACCGGCCACCGCGAAccc                                                                                                >  2:184213/1‑150 (MQ=37)
                                                      gggAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAAcc                                                                                                >  2:2112695/1‑150 (MQ=255)
                                                      gggAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAAcc                                                                                                >  2:7115340/1‑150 (MQ=255)
                                                         aaCCTTCGGGTTCCTT‑‑TTTTTTTTTTTTTTGATTGTGTGTTGTCTGGGCGGGCCCTTTGGGGGTTAAACCCGCTCCTACAGGTTTTTGGCGGCCCCTGTAGGTGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTa                                                                                             <  1:2606872/150‑1 (MQ=255)
                                                            ccgaaggTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGTTGCCCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCACCCGCACCCGAGCAACCCGAcgg                                                                                          >  2:7653870/6‑148 (MQ=37)
                                                            cTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCACCTACAGGATTATCGCTGCCCCGGTAGGAGCGGGTTTACCCCCCAAGAGGCCAGCAACCGCAACCGAGAAAACCAAcgg                                                                                          >  2:511240/1‑148 (MQ=255)
                                                               ggggTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTgc                                                                                       <  2:4726546/149‑1 (MQ=255)
                                                                  gTTCCCT‑‑TTTTTTTTAATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTCCAGGATTATCGCGGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGca                                                                                    <  2:3766217/150‑1 (MQ=255)
                                                                             ttttttttttttTGTTTGTGTGTTGTCTGGGCCGGCCGTTGGGCGGTTAAACCCGCCCATACCGTTTTATCGCTGCCCCTTTTGGTGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTa                                                                            <  1:2243263/150‑1 (MQ=255)
                                                                             ttttttttATTCTGATTGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTa                                                                            <  2:359238/150‑1 (MQ=255)
                                                                               ttttttATACTGAAGGTGGGTTGTCAGTGCCGGCCATTACGCGGGTAAACCCGCTCCAACCGTTTTATCGCTGCCCCGGTAGGCGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACt                                                                          <  1:4338145/150‑1 (MQ=37)
                                                                                tttttATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCAACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTAACCCCGAAGAGGCCAGCACCGGCAACCCAGAAAACCTACTGTGCGCACCCCCCGACTg                                                                         >  2:6968900/1‑150 (MQ=38)
                                                                                         ttGATTGTGTGTTTTCTGTGCCGGCCGTTTCGCGGTTAAACCCGCTCCTACCGGTTTCTCGCTGCCCGTGTTGGTGGGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTct                                                                 <  1:5563014/149‑1 (MQ=37)
                                                                                              tgtgtgTTGTCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCATCCTGTTTTATCGCTGCCCCTGTTGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCtgctgc                                                            <  1:3395158/150‑1 (MQ=38)
                                                                                               gtgtgtTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGATTTTAGCTGCCCCTGGAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGGAACCCGACTGTGCGGCCAGCCCTACTGGTCCTTCTgctgct                                                           >  2:7479641/1‑150 (MQ=37)
                                                                                                tgtgtTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTa                                                          <  1:5861383/150‑1 (MQ=255)
                                                                                                   tttttCAGGGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACTGTTTTATCGCGGCCCCTGTTGTAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTc                                                       <  1:1681083/146‑1 (MQ=38)
                                                                                                    ttgtCTGTGCGGGCCCTTTCGGGGTTAACCCCGCACCTACTGTTTTTGCGGTGCCCCTTTTGGCGGGGTTTAACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCa                                                      <  1:213915/150‑1 (MQ=255)
                                                                                                    ttgtCAGTGCCGGCCCTTTCGCGGGTAAAACCCCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTAATGTTCATTCTGCTGCTACTCa                                                      >  2:2582498/1‑150 (MQ=255)
                                                                                                        ctgTGCCGGCCCTTTCGGGGGTAAACCCGCTCCTACAGGTTTTTCGCTGCCCGTGTTGGTGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGc                                                  <  1:7035238/148‑1 (MQ=38)
                                                                                                          gTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGTTTTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCt                                                <  2:1053422/150‑1 (MQ=255)
                                                                                                          gTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCGGCCCCTGTAGGAGCGGGTTTACCCGCGAAGTGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCGGCTACTCACCGccg                                                >  2:1241710/1‑149 (MQ=255)
                                                                                                          gTGCCGGCCCTTTCGCGGGTAAACCCGCTCCAACAGGATTATCGCATCCCCAGTAGGAGCGGGTTTACCCGCCGAGAGGCCAGCACCGGCAACCCAGAAAACCAACGGTGAGCCAAGCCCTACTGTTCAATCTGCAGCGACACACCGcca                                                >  2:4018655/1‑149 (MQ=21)
                                                                                                                  ccTTTCGCGGGTAAACCCGCTCCTACAGTATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCaa                                        <  2:6804515/150‑1 (MQ=255)
                                                                                                                  ccTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCaa                                        <  2:4535471/150‑1 (MQ=255)
                                                                                                                    tttCGCGGGTAAACCCCCTCCAACACGAATATCGCATCCCCAGTAGGAGCCGGTTTACCCGCGAAGAGGCCAGCACCGGCAAACGAGAAACCCAACAGTGCGCACAGCCCGACTGTTCATTCTGCTGCTACTCACCCCCTGCAACCaacg                                      >  2:107274/1‑148 (MQ=37)
                                                                                                                     ttCGCGGTTAAACCCGCTCCTACAGTTTTAGCGCTGCCCGTGTTGGAGGGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGc                                     <  1:6056512/150‑1 (MQ=39)
                                                                                                                     ttCGCGGGTAACCCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGc                                     <  2:1614951/150‑1 (MQ=255)
                                                                                                                     ttCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGGAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGc                                     >  1:3577959/1‑150 (MQ=255)
                                                                                                                         cgGGCAACCCCGCTCCTCCTCGATTATCGCCGCCCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCa                                 <  1:1958739/150‑1 (MQ=39)
                                                                                                                                aCCCGCTCCTACTGGTTTAGCGCGCCCCCTTTAGGCGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGa                          <  1:3054626/150‑1 (MQ=255)
                                                                                                                                aCCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGa                          <  2:4328637/150‑1 (MQ=255)
                                                                                                                                  cggCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgc                        >  1:3095831/3‑150 (MQ=255)
                                                                                                                                  ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgc                        >  1:8475246/1‑150 (MQ=255)
                                                                                                                                    gCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAgcgc                      >  1:8617858/1‑150 (MQ=255)
                                                                                                                                     cTCCTACCGGATTATCGCGGCCCCTGTCGGTGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGcc                     <  1:859436/150‑1 (MQ=255)
                                                                                                                                      tCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGa                          >  2:8800355/1‑144 (MQ=255)
                                                                                                                                        cTACTGGATTATCGCTGCCCCTTTTGGTGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTgg                  <  1:3431026/150‑1 (MQ=255)
                                                                                                                                         taccgttttATCGTTGCCCGTGTCGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                 <  1:7653870/143‑1 (MQ=39)
                                                                                                                                         tACCGGTTTCTCGCGGCCCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                 <  1:2795570/150‑1 (MQ=39)
                                                                                                                                         tACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                 >  1:5141439/1‑150 (MQ=255)
                                                                                                                                         tACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                 >  2:4133490/1‑150 (MQ=255)
                                                                                                                                         tACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                 >  2:3961548/1‑150 (MQ=255)
                                                                                                                                            aGGATTATCGCTGCCCCTGTAGGAGCGGGTTTATCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtgg              <  2:7234767/150‑1 (MQ=255)
                                                                                                                                              gattcgcgcgGCCCCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTCCGGTGCGCACAGCCCTACGGTTCAGTCTGCTGCGACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTc            <  2:4691654/144‑1 (MQ=38)
                                                                                                                                              gatacgcgcgGCCCCTGTTGGAGGGGGTTAACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCAGTCTGCTGCTACTCCCCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTc            <  1:856752/144‑1 (MQ=255)
                                                                                                                                              gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTc            >  1:3300707/1‑150 (MQ=255)
                                                                                                                                              gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTc            >  2:6877774/1‑150 (MQ=255)
                                                                                                                                                ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCAGCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCg          >  2:5482849/1‑150 (MQ=255)
                                                                                                                                                ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCg          >  2:6735661/1‑150 (MQ=255)
                                                                                                                                                ttATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCg          >  2:7310790/1‑150 (MQ=255)
                                                                                                                                                  aTCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAc        <  2:2036955/150‑1 (MQ=255)
                                                                                                                                                   tCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGACACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAcc       >  2:611728/1‑150 (MQ=255)
                                                                                                                                                   tCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAcc       >  2:5365168/1‑150 (MQ=255)
                                                                                                                                                   tCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAcc       >  2:8691955/1‑150 (MQ=255)
                                                                                                                                                       tGCCCCTTTATGCGGGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATgg   <  1:440981/150‑1 (MQ=255)
                                                                                                                                                       tGCCCCTGTCGGTGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATgg   <  1:2309729/150‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                    
GTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCG‑‑TTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGG  >  NC_002947/2461713‑2462010

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: