Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,990,335 G→C 5.1% intergenic (+35/+59) ycgR → / ← PP_4398 flagellar brake protein YcgR/MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,990,3350GC 5.1% 516.3 / 12.5 155intergenic (+35/+59)ycgR/PP_4398flagellar brake protein YcgR/MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (61/86);  new base C (4/4);  total (65/90)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 164 total aligned reads displayed.

ACCTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCTGCATC  >  NC_002947/4990187‑4990484
                                                                                                                                                    |                                                                                                                                                     
accTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGa                                                                                                                                                      <  1:4874181/150‑1 (MQ=255)
accTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGa                                                                                                                                                      <  1:2927529/150‑1 (MQ=255)
accTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGa                                                                                                                                                      <  1:275883/150‑1 (MQ=255)
   tttGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGt                                                                                                                                                   <  2:573981/150‑1 (MQ=255)
   tttGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGt                                                                                                                                                   >  2:5284464/1‑150 (MQ=255)
   tttGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGt                                                                                                                                                   <  1:1355919/150‑1 (MQ=255)
     tGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTgg                                                                                                                                                 <  1:2975344/150‑1 (MQ=255)
      gCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGc                                                                                                                                                <  2:790747/150‑1 (MQ=255)
       ccGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCg                                                                                                                                               >  1:2450742/1‑150 (MQ=255)
        cGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCGGAGGTGGCGt                                                                                                                                              >  2:3403237/1‑150 (MQ=255)
          gcgtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGttt                                                                                                                                            <  1:981971/150‑1 (MQ=255)
          gcgtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGttt                                                                                                                                            <  2:606370/150‑1 (MQ=255)
          gcgtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGttt                                                                                                                                            <  1:982082/150‑1 (MQ=255)
            gtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTGGGTCGGTAAAACAGAACACGCCACCTCAGCAGGTGGCGTTTTc                                                                                                                                          <  1:1379211/150‑1 (MQ=255)
            gtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTc                                                                                                                                          <  1:3618749/150‑1 (MQ=255)
            gtgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTCAGCAGGTGGCGTTTTc                                                                                                                                          <  1:4534998/150‑1 (MQ=255)
             tgcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCa                                                                                                                                         >  1:5049025/1‑150 (MQ=255)
              gcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCAt                                                                                                                                        <  1:2543104/150‑1 (MQ=255)
              gcgtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCAt                                                                                                                                        >  2:467709/1‑150 (MQ=255)
               cgtTTCCACACCCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCAtt                                                                                                                                       <  2:1530076/150‑1 (MQ=255)
                gtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGTTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCAttt                                                                                                                                      <  2:490710/150‑1 (MQ=255)
                gtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCAttt                                                                                                                                      <  2:488270/150‑1 (MQ=255)
                gtTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCAttt                                                                                                                                      <  1:3391414/150‑1 (MQ=255)
                 tttCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGtttt                                                                                                                                           <  1:636371/144‑1 (MQ=255)
                 tttCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTtg                                                                                                                                     >  2:4824489/1‑150 (MQ=255)
                 tttCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTGTTCAGTtg                                                                                                                                     >  2:4485140/1‑150 (MQ=255)
                  ttCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTtgt                                                                                                                                    >  2:1684249/1‑150 (MQ=255)
                  ttCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTtgt                                                                                                                                    <  2:3474910/150‑1 (MQ=255)
                  ttCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACCATTATCGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGATTTCATTtgt                                                                                                                                    <  2:2079560/150‑1 (MQ=255)
                    ccACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTgtgt                                                                                                                                  >  1:4496360/1‑150 (MQ=255)
                     cacaACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGtt                                                                                                                                 >  1:1008196/1‑150 (MQ=255)
                      acaACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGttt                                                                                                                                >  1:2108413/1‑150 (MQ=255)
                           cTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTctgct                                                                                                                           >  1:3138474/1‑150 (MQ=255)
                            tGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTctgctg                                                                                                                          >  2:2673817/1‑150 (MQ=255)
                            tGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTctgctg                                                                                                                          <  1:3168768/150‑1 (MQ=255)
                            tGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTGTCATTTGTGTTTctgctg                                                                                                                          >  2:2673575/1‑150 (MQ=255)
                                    cAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGcc                                                                                                                  <  1:1198142/150‑1 (MQ=255)
                                      ggCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGc                                                                                                                <  2:4304933/150‑1 (MQ=255)
                                            gcgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTcgcg                                                                                                          <  1:1684249/150‑1 (MQ=255)
                                             cgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCgg                                                                                                         >  1:3028531/1‑150 (MQ=255)
                                             cgcaAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCgg                                                                                                         <  1:811191/150‑1 (MQ=255)
                                               aaaGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTAGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCAGGGATGGCAGCTTCGCggga                                                                                                       >  2:2185694/2‑149 (MQ=255)
                                                   aTCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAAc                                                                                                   <  1:1133966/150‑1 (MQ=255)
                                                     cGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAAccc                                                                                                 >  2:3567683/1‑150 (MQ=255)
                                                      gagaGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCg                                                                                                >  2:3739559/1‑150 (MQ=255)
                                                      gagaGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTCAGCAGGTGGCGTTTTct                                                                                                                                         >  2:3089538/1‑108 (MQ=38)
                                                       agagTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGCTTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGc                                                                                               <  2:4005104/150‑1 (MQ=255)
                                                       agagTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTGTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGc                                                                                               >  2:3994295/1‑150 (MQ=255)
                                                         agTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTc                                                                                             <  2:1140559/150‑1 (MQ=255)
                                                         agTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTc                                                                                             <  2:327391/150‑1 (MQ=255)
                                                           tttCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCt                                                                                           <  2:1965779/150‑1 (MQ=255)
                                                           tttCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCt                                                                                           <  2:4704537/150‑1 (MQ=255)
                                                           tttCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCt                                                                                           <  2:1696445/150‑1 (MQ=255)
                                                             tCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTac                                                                                         >  2:1019664/1‑150 (MQ=255)
                                                                tgtACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACAcc                                                                                      <  2:5059329/150‑1 (MQ=255)
                                                                 gtACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCg                                                                                     >  1:5105006/1‑150 (MQ=255)
                                                                  tACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGa                                                                                    <  1:3409238/150‑1 (MQ=255)
                                                                     cAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACtt                                                                                 >  1:1994319/1‑150 (MQ=255)
                                                                     cAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACtt                                                                                 >  2:1399880/1‑150 (MQ=255)
                                                                        cTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCt                                                                              <  1:2217461/150‑1 (MQ=255)
                                                                         tGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTg                                                                             <  1:4751008/150‑1 (MQ=255)
                                                                          gCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGa                                                                            <  2:3665399/150‑1 (MQ=255)
                                                                           cAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGAc                                                                           <  2:4473675/150‑1 (MQ=255)
                                                                           cAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTAGCGGGTAAACCCGCTCCTACACCGACTTGCTGAc                                                                           >  1:1041188/1‑150 (MQ=255)
                                                                                  aGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCa                                                                    <  1:3948657/150‑1 (MQ=255)
                                                                                       cggcggTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGgcgc                                                               <  1:2875046/150‑1 (MQ=255)
                                                                                             ttCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtg                                                         <  1:41306/150‑1 (MQ=255)
                                                                                             ttCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTtg                                                         <  2:3466625/150‑1 (MQ=255)
                                                                                             ttCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTCAGCAGGTGGCGtt                                                                                                                                             <  2:2944107/66‑1 (MQ=25)
                                                                                                     ggACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTgcgcgc                                                 <  2:4261516/150‑1 (MQ=255)
                                                                                                         gATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGCCGTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGcccgt                                             >  1:1913982/1‑148 (MQ=255)
                                                                                                           ttattGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTtct                                           <  2:2445643/150‑1 (MQ=255)
                                                                                                            tattGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTtctc                                          <  1:1481885/150‑1 (MQ=255)
                                                                                                              ttGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCtt                                        >  1:3621873/1‑150 (MQ=255)
                                                                                                              ttGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCtt                                        <  1:1313914/150‑1 (MQ=255)
                                                                                                                gATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCGTCGCGGGTAAAACCGCTCCTACACCCACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCGTCTCTTCa                                      >  1:1957832/1‑150 (MQ=255)
                                                                                                                  tcgtcgGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGTCCTCTTCAGGCGCCTGTTGCTGCGCGCCCGGCTCTTCATc                                    >  1:3486623/1‑150 (MQ=255)
                                                                                                                     tcggtcggTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGATGCTCGCCCTTCTCTTCATCCTc                                 <  2:5159367/150‑1 (MQ=255)
                                                                                                                        gtcggtAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTc                                                                                             <  1:2996589/87‑1 (MQ=255)
                                                                                                                           ggtAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTCAACCCGCGCCGACACCGACTTGCTGACCGCTTCAGGCGCCGGTTGCTGCGCGCCCGTCTCTTCATCCTCGGTGTc                           >  1:5134848/1‑150 (MQ=255)
                                                                                                                            gtAAAAGTGTAAACGGCACCTGCTGTGGTGGAGTTTTCATTTGTGTTTCTGGTGGGCTGGCCGCTTCGCGGGTAAACAGGCTCGTACTCGGAGTTGGTGTCCTCTTCAGGCGCGGGTTGCTGCGCGCCCTTCTCTTCTGCGTGGGTGTCt                          <  2:1332438/150‑1 (MQ=255)
                                                                                                                            gtAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCt                          <  1:467709/150‑1 (MQ=255)
                                                                                                                               aaaCAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCa                       >  1:4285939/1‑150 (MQ=255)
                                                                                                                               aaaCAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCa                       <  2:5049025/150‑1 (MQ=255)
                                                                                                                                aaCAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAg                      >  2:3883349/1‑150 (MQ=255)
                                                                                                                                 aCAGAAATCGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAgc                     >  2:2128676/1‑150 (MQ=255)
                                                                                                                                 aCAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAgc                     >  1:4421021/1‑150 (MQ=255)
                                                                                                                                  cAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAgcc                    >  2:4388591/1‑150 (MQ=255)
                                                                                                                                  atgAAAACGCCACCTCAGCAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAgcc                    >  2:4758350/3‑150 (MQ=255)
                                                                                                                                      aaaCGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCa                >  2:221289/1‑150 (MQ=255)
                                                                                                                                       aaCGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAAACCGCTCCTACAACGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCGTCTCTTCATCCTCGGTGTCTTCAGCCGCCAc               >  1:1162165/1‑150 (MQ=255)
                                                                                                                                        aCGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCAcc              <  1:1627317/150‑1 (MQ=255)
                                                                                                                                         cGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCAACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCGTCTCTTCATACTCGGTGTCTTCAGCCGCCACCg             >  2:1604767/1‑150 (MQ=255)
                                                                                                                                             aCCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGgcgc         <  1:5297005/150‑1 (MQ=255)
                                                                                                                                               ctgctgAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCt       <  1:3991283/150‑1 (MQ=255)
                                                                                                                                                 gctgAGGTGGCGTTTTCATTTGTGTTTATGCTGGGCTGGCCGCTTCGCGGGTAAAACCCCGCCAACAACCACATGCAGACCTCGTCAGGCGCCGGTTGCTTCGCGCCCCTCACTTAATACTCGGTGTCTTCAGCCGCCACCCGCCCCAGc     >  1:596068/1‑150 (MQ=255)
                                                                                                                                                   tgAGGTGGCGTTTTCATTTGTTTTTATGCTGGGCTGGCCGCGTCGCCGGTAAACCCGCACCCACACCCACAAGCTGACCTAATAAGGCCCCAGTTGCTGCGCGCCCATCAATTCAAACACGGGGTAATAATCCGCCACCCGCCCCAGCa  >  1:2951158/1‑149 (MQ=255)
                                                                                                                                                   tgAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACtt                                                                                 >  1:4523886/1‑72 (MQ=255)
                                                                                                                                                   tgAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACtt                                                                                 <  2:4523886/72‑1 (MQ=255)
                                                                                                                                                    gAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCTGCATc  <  2:476634/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACCTTTGCCGGCGTGCGTTTCCACAACCTGAGCGGGCAGGCGCAGCGCAAGATCGAGAGTTTCGTGTACCAGCTGCAGCGTGAGGCGCGGCGGTTCGACAAGGACGATTATTGATCGTCGGTCGGTAAAACAGAAAACGCCACCTGCTGAGGTGGCGTTTTCATTTGTGTTTCTGCTGGGCTGGCCGCTTCGCGGGTAAACCCGCTCCTACACCGACTTGCTGACCTCTTCAGGCGCCTGTTGCTGCGCGCCCTTCTCTTCATCCTCGGTGTCTTCAGCCGCCACCGGCGCCTGCATC  >  NC_002947/4990187‑4990484

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