Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,728,698 A→G 6.7% intergenic (+50/+73) cspA‑I → / ← PP_1523 major cold shock protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,728,6980AG 6.7% 413.2 / 11.2 134intergenic (+50/+73)cspA‑I/PP_1523major cold shock protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (59/65);  new base G (4/5);  total (63/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

AGGGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCG  >  NC_002947/1728554‑1728846
                                                                                                                                                |                                                                                                                                                       
aGGGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTg                                                                                                                                                    <  1:3968965/150‑1 (MQ=255)
aGGGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTg                                                                                                                                                    <  1:2439848/150‑1 (MQ=255)
    gCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTg                                                                                                                                                >  1:3907434/1‑150 (MQ=255)
    gCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTg                                                                                                                                                <  2:3916165/150‑1 (MQ=255)
     cATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGc                                                                                                                                               <  1:3994931/150‑1 (MQ=255)
     cATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGc                                                                                                                                               <  2:3167955/150‑1 (MQ=255)
      aTCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGc‑‑c                                                                                                                                            >  2:2075675/1‑150 (MQ=255)
       tCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑Ca                                                                                                                                           >  2:450869/1‑150 (MQ=255)
           aCCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCGTGGTGGCTGcttcgcg                                                                                                                                         <  2:2916838/150‑7 (MQ=255)
           aCCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCGTGGTGGCTGcttcgcg                                                                                                                                         <  2:2916861/150‑7 (MQ=255)
           aCCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACggg                                                                                                                                       <  1:4800792/150‑1 (MQ=255)
            cccTGGTCGCAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCGTGCCCGCGAAGCAGCCACCGTGGTGGCTGcttcgcgg                                                                                                                                        <  2:3140567/150‑8 (MQ=255)
             ccTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCa                                                                                                                                     >  2:4180439/1‑150 (MQ=255)
               tGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATa                                                                                                                                   <  1:288101/150‑1 (MQ=255)
                 gTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACTTATTGGGCATCTGAGCATCTTGTTTCAAGCATTAGGAGCGGGCTTGCCAGCTAAGCAGCCAACACGGTGGCTGC‑‑CACGGCCATAGc                                                                                                                                 <  1:3012353/150‑1 (MQ=255)
                   cGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGccc                                                                                                                               <  2:2521254/150‑1 (MQ=255)
                   cGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGccc                                                                                                                               <  2:2521489/150‑1 (MQ=255)
                    gAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCg                                                                                                                              >  1:2225338/1‑150 (MQ=255)
                           aGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTcgc                                                                                                                       <  2:193785/150‑1 (MQ=255)
                              gggTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCggg                                                                                                                    >  2:3370921/1‑150 (MQ=255)
                               ggTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGt                                                                                                                   >  1:295869/1‑150 (MQ=255)
                                   gAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAAc                                                                                                               >  1:3227509/1‑150 (MQ=255)
                                      tACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCa                                                                                                            <  1:1444545/150‑1 (MQ=255)
                                       aCGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCAt                                                                                                           <  2:1713035/150‑1 (MQ=255)
                                           cTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCt                                                                                                       >  2:4366470/1‑150 (MQ=255)
                                           cTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCt                                                                                                       >  2:2827209/1‑150 (MQ=255)
                                           cTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCt                                                                                                       <  1:3914794/150‑1 (MQ=255)
                                              cGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACa                                                                                                    <  1:2504152/150‑1 (MQ=255)
                                                tGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGa                                                                                                  >  1:209876/1‑150 (MQ=255)
                                                   aGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCgg                                                                                                                     >  2:425920/1‑128 (MQ=255)
                                                   aGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCgg                                                                                                                     <  1:426600/128‑1 (MQ=255)
                                                   aGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAAc                                                                                               >  2:3346223/1‑150 (MQ=255)
                                                   aGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAAc                                                                                               <  1:303732/150‑1 (MQ=255)
                                                   aGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAAc                                                                                               <  2:1302559/150‑1 (MQ=255)
                                                       ccAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAATCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGc                                                                                           >  1:3458576/1‑150 (MQ=255)
                                                       ccAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGc                                                                                           >  2:1998860/1‑150 (MQ=255)
                                                       ccAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGc                                                                                           >  1:1011280/1‑150 (MQ=255)
                                                         aGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCaa                                                                                         <  1:3665394/150‑1 (MQ=255)
                                                           aaaGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGc                                                                                       <  2:2041888/150‑1 (MQ=255)
                                                             aGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCt                                                                                     >  2:1858362/1‑150 (MQ=255)
                                                                ccTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGcc                                                                                 <  1:391848/150‑1 (MQ=255)
                                                                ccTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGcc                                                                                 <  1:4913736/150‑1 (MQ=255)
                                                                    cAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGACGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTTCCGttt                                                                             <  2:3659175/150‑1 (MQ=255)
                                                                      aGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCa                                                                           <  2:4840657/150‑1 (MQ=255)
                                                                       gCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCat                                                                          <  1:4509497/150‑1 (MQ=255)
                                                                       gCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCat                                                                          <  2:4202163/150‑1 (MQ=255)
                                                                       gCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCat                                                                          <  2:4803269/150‑1 (MQ=255)
                                                                        ccGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCatg                                                                         <  2:1316377/150‑1 (MQ=255)
                                                                        ccGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGCTTCatg                                                                         <  1:1148478/150‑1 (MQ=255)
                                                                        ccGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGGCAGGCTTGCTCGTTTctgt                                                                         <  1:3854792/150‑4 (MQ=39)
                                                                         cGAAGACGTTGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTTGGAGCGGGCTTGCCCGCGAAGCAGCCACCTCGGTGGCTGC‑‑CACGGGCTTAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGGCAGGCTTGCTCGTTTctgta                                                                        <  2:3734038/150‑5 (MQ=37)
                                                                          gAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCatgat                                                                       >  2:896318/1‑150 (MQ=255)
                                                                           aaGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCatgatg                                                                      <  1:1298293/150‑1 (MQ=255)
                                                                              aCGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTAc                                                                   <  2:1415111/150‑1 (MQ=255)
                                                                              aCGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTAc                                                                   <  2:2261154/150‑1 (MQ=255)
                                                                               cGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACg                                                                  <  2:2769385/150‑1 (MQ=255)
                                                                                gTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCGtggtgg                                                                                                                                                   <  2:302708/71‑1 (MQ=255)
                                                                                gTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCGtggtgg                                                                                                                                                   >  1:302708/1‑71 (MQ=255)
                                                                                 tAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGcc                                                                >  2:5073714/1‑150 (MQ=255)
                                                                                 tAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGcc                                                                >  2:5073640/1‑150 (MQ=255)
                                                                                   gTGGGCCTCTGAGCCTCGGGTTTCAAGCCGTAGGAGCGGGCGTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAg                                                              <  1:156826/150‑1 (MQ=255)
                                                                                     gggCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGt                                                            <  2:3820611/150‑1 (MQ=255)
                                                                                        ccTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGAt                                                         >  1:766943/1‑150 (MQ=255)
                                                                                         ctctGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGAtt                                                        >  2:3604160/1‑150 (MQ=255)
                                                                                           ctGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTc                                                      <  1:669039/150‑1 (MQ=255)
                                                                                            tGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTcc                                                     >  2:2101728/1‑150 (MQ=255)
                                                                                            tGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTcc                                                     >  2:3062354/1‑150 (MQ=255)
                                                                                            tGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTcc                                                     >  2:506740/1‑150 (MQ=255)
                                                                                                ccTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCtt                                                 <  1:4274093/150‑1 (MQ=255)
                                                                                                 ctctGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTc                                                <  1:3824557/150‑1 (MQ=255)
                                                                                                 ctctGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTc                                                >  2:1435341/1‑150 (MQ=255)
                                                                                                  tctGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCa                                               <  1:2863057/150‑1 (MQ=255)
                                                                                                  tctGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCa                                               >  1:2931433/1‑150 (MQ=255)
                                                                                                    tGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCAcc                                             >  2:2684801/1‑150 (MQ=255)
                                                                                                      gTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGt                                           >  2:3789229/1‑150 (MQ=255)
                                                                                                       gggcAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTGCAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTc                                          >  1:4827286/4‑150 (MQ=255)
                                                                                                        ttCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCg                                         <  1:3519570/150‑1 (MQ=255)
                                                                                                         tCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCgg                                        >  1:1021288/1‑150 (MQ=255)
                                                                                                          cAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGgc                                       <  1:4180439/150‑1 (MQ=255)
                                                                                                          cAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGgc                                       <  1:1998860/150‑1 (MQ=255)
                                                                                                           aaGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATGTCCTCTTCACCGTCGgcg                                      >  2:855026/1‑150 (MQ=255)
                                                                                                                gTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGAt                                 >  1:514256/1‑150 (MQ=255)
                                                                                                                    gAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCgg                             >  1:2977425/1‑150 (MQ=255)
                                                                                                                     aGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGa                            >  1:5032438/1‑150 (MQ=255)
                                                                                                                      gCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGAc                           <  1:450869/150‑1 (MQ=255)
                                                                                                                           tttacccGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGc                      >  1:488056/5‑150 (MQ=255)
                                                                                                                           tttacccGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGc                      >  1:3042927/5‑150 (MQ=255)
                                                                                                                              gCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGt                   <  2:3058281/150‑1 (MQ=255)
                                                                                                                                ccGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCa                 >  1:4109748/1‑150 (MQ=255)
                                                                                                                                 cgcgAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCACGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAg                >  1:2733271/1‑150 (MQ=255)
                                                                                                                                 cgcgAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAg                >  2:1396471/1‑150 (MQ=255)
                                                                                                                                  ccgtggcagcCACCGTGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGc               >  2:1782609/6‑150 (MQ=255)
                                                                                                                                      ggcagcCACCGTGGTGGCTGCTTCGCGGGCAAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTGTCATGGTGTACGCCCGGTGATTTCCTCTTCACCGGCGGCGCTGATGCGGACCCAGCGGTCAGCAtc            >  1:4354493/2‑150 (MQ=37)
                                                                                                                                        cagcCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCAtcttc         <  2:1122344/150‑1 (MQ=255)
                                                                                                                                          gcCACCGTGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATcttctt       >  2:1645492/1‑150 (MQ=255)
                                                                                                                                           ccaccaCGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCttcttc      >  2:629805/1‑150 (MQ=255)
                                                                                                                                            caccacGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGAGGTATGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTcc     >  2:597087/1‑150 (MQ=255)
                                                                                                                                             accacGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCatg                                                                         >  1:2612887/1‑81 (MQ=255)
                                                                                                                                              ccacGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTcccc   >  1:2868233/1‑150 (MQ=255)
                                                                                                                                               cacGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCg  >  2:3910187/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                       
AGGGGCATCGTACCCTGGTCGAAGGGCAGCGGGTGGAGTACGCCTGCGTGCAGGGCCAGAAAGGCCTGCAAGCCGAAGACGTAGTGGGCCTCTGAGCCTCTGGTTTCAAGCCGTAGGAGCGGGCTTGCCCGCGAAGCAGCCACCACGGTGGCTGC‑‑CACGGGCATAGCCCGTATTCGCGGGTAAACCCATTCCTACAGAAACGAGCAAGCCTGCCGTTTCATGATGTACGCCAGGTGATTTCCTCTTCACCGTCGGCGCTGATGCGGACCCAGCGGTCAGCATCTTCTTCCCCG  >  NC_002947/1728554‑1728846

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