Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 268,746 G→C 6.7% intergenic (‑77/+23) PP_0216 ← / ← desA sensory box/GGDEF family protein/delta‑9 fatty acid desaturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947268,7460GC 6.7% 848.0 / 38.4 252intergenic (‑77/+23)PP_0216/desAsensory box/GGDEF family protein/delta‑9 fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base G (107/128);  new base C (8/9);  total (115/137)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 272 total aligned reads displayed.

CGGGGAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGGCCGTTG  >  NC_002947/268597‑268895
                                                                                                                                                     |                                                                                                                                                     
cGGGGAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTgg                                                                                                                                                       >  1:5430264/1‑150 (MQ=255)
cGGGGAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTgg                                                                                                                                                       <  2:3376168/150‑1 (MQ=255)
  gggAAAATCGTTGGAATGTTCTTGGCGATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTgggg                                                                                                                                                     <  2:749038/150‑1 (MQ=255)
  gggAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGTACGCCCCGCTgggg                                                                                                                                                     <  1:4245379/150‑1 (MQ=255)
       aaTCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGt                                                                                                                                                <  1:1184859/150‑1 (MQ=255)
        aTCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTg                                                                                                                                               <  1:1780080/150‑1 (MQ=255)
        aTCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCCCAGCGGGGCGTg                                                                                                                                               <  1:4129663/150‑1 (MQ=255)
         tCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGCAGTTCCGAATGGCACGCCCCGCTTGGGGGGGTGc                                                                                                                                              <  2:1301247/150‑1 (MQ=255)
            ttGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGa                                                                                                                                           >  2:3234428/1‑150 (MQ=255)
             tGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGAt                                                                                                                                          <  1:558385/150‑1 (MQ=255)
              ggAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATc                                                                                                                                         >  1:1384069/1‑150 (MQ=255)
              ggAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCCCAGCGGGGCGTGccattc                                                                                                                                         <  1:807988/150‑5 (MQ=255)
              ggAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGCGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGCGTGccattcg                                                                                                                                        <  1:5350851/150‑6 (MQ=255)
               gAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCa                                                                                                                                        >  1:4930151/1‑150 (MQ=255)
               gAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCa                                                                                                                                        >  1:4542249/1‑150 (MQ=255)
                aaTGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAACGGCACGCCCCCCAGCGGGGCGTGccattcgg                                                                                                                                       <  2:351592/150‑7 (MQ=255)
                   gTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAcgg                                                                                                                                    <  1:3720288/150‑1 (MQ=255)
                    ttCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAcggc                                                                                                                                   <  2:3324354/150‑1 (MQ=255)
                     tCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAcggcg                                                                                                                                  <  1:3823164/150‑1 (MQ=255)
                      cTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAAcggcgg                                                                                                                                 <  2:2218331/150‑1 (MQ=255)
                        tGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTc                                                                                                                               <  1:2399318/150‑1 (MQ=255)
                        tGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTc                                                                                                                               <  1:4032551/150‑1 (MQ=255)
                         ggCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCa                                                                                                                              >  1:2391849/1‑150 (MQ=255)
                           ccATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAgg                                                                                                                            <  1:1319281/150‑1 (MQ=255)
                                      gAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCggg                                                                                                                 >  1:2500737/1‑150 (MQ=255)
                                          aGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCa                                                                                                             <  1:5580902/150‑1 (MQ=255)
                                            gTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGc                                                                                                           >  1:4659992/1‑150 (MQ=255)
                                             tctcAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCCGACCGgcgcctga                                                                                                          <  1:850438/150‑5 (MQ=255)
                                               tcAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAg                                                                                                        >  1:3611541/1‑150 (MQ=255)
                                               tcAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAg                                                                                                        <  2:2973020/150‑1 (MQ=255)
                                                cAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGt                                                                                                       >  2:4806417/1‑150 (MQ=255)
                                                 aaCGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTa                                                                                                      >  1:2564321/1‑150 (MQ=255)
                                                 aaCGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTa                                                                                                      >  1:2852002/1‑150 (MQ=255)
                                                  aCGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTaa                                                                                                     <  2:5115220/150‑1 (MQ=255)
                                                  aCGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTaa                                                                                                     <  1:3234428/150‑1 (MQ=255)
                                                  aCGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTaa                                                                                                     <  1:5491216/150‑1 (MQ=255)
                                                  aCGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTaa                                                                                                     >  2:4368552/1‑150 (MQ=255)
                                                    gTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGt                                                                                                   >  2:775882/1‑150 (MQ=255)
                                                     tGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTc                                                                                                  <  2:3589137/150‑1 (MQ=255)
                                                     tGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTc                                                                                                  <  2:4025852/150‑1 (MQ=255)
                                                         aGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGGTATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGACTGGCACGCCCCGCTGGGGGGCGTGCCGATCTACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGa                                                                                              <  2:3195564/150‑1 (MQ=255)
                                                         aGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGa                                                                                              >  2:339905/1‑150 (MQ=255)
                                                         aGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGa                                                                                              <  2:4548838/150‑1 (MQ=255)
                                                         aGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGa                                                                                              <  1:3558062/150‑1 (MQ=255)
                                                            gTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGtt                                                                                           <  2:2814786/150‑1 (MQ=255)
                                                                ggCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGtgc                                                                                       <  2:415391/150‑1 (MQ=255)
                                                                 gCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGgcgc                                                                                                                         <  1:530711/115‑1 (MQ=255)
                                                                 gCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGtgct                                                                                      <  1:4418520/150‑1 (MQ=255)
                                                                 gCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGtgct                                                                                      >  1:520859/1‑150 (MQ=255)
                                                                 gCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGtgct                                                                                      >  1:4060442/1‑150 (MQ=255)
                                                                  ccAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGtgctg                                                                                     >  1:1101145/1‑150 (MQ=255)
                                                                   cAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGtgctgc                                                                                    >  2:1927519/1‑150 (MQ=255)
                                                                   cAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGtgctgc                                                                                    <  2:1533977/150‑1 (MQ=255)
                                                                   cAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGtgctgc                                                                                    <  1:379446/150‑1 (MQ=255)
                                                                      cATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGaa                                                                                 <  2:2640/150‑1 (MQ=255)
                                                                      cATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGaa                                                                                 <  2:3852286/150‑1 (MQ=255)
                                                                      cATCCGGCCCGCCACGCAGGCGATAGTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCCCGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGaa                                                                                 <  2:3609748/150‑1 (MQ=255)
                                                                        tCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGt                                                                               >  2:1545017/1‑150 (MQ=255)
                                                                                gCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGt                                                                       >  2:2070200/1‑150 (MQ=255)
                                                                                     gCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTg                                                                  <  1:3943337/150‑1 (MQ=255)
                                                                                      cAGGCGATATGCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTgg                                                                 <  1:3133553/150‑1 (MQ=255)
                                                                                          cGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGc                                                             <  2:873114/150‑1 (MQ=255)
                                                                                           gATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGcc                                                            >  1:2119660/1‑150 (MQ=255)
                                                                                               gtCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTg                                                        <  2:4963458/149‑1 (MQ=255)
                                                                                                       tCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGc                                                <  1:1655713/150‑1 (MQ=255)
                                                                                                       tCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGc                                                <  2:1418581/150‑1 (MQ=255)
                                                                                                            tGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCAt                                           <  2:4174970/150‑1 (MQ=255)
                                                                                                              gagaGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGt                                         <  1:4638048/150‑1 (MQ=255)
                                                                                                                gagGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTAc                                       <  1:1979300/150‑1 (MQ=255)
                                                                                                                   gTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTc                                             >  2:127754/1‑141 (MQ=255)
                                                                                                                   gTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCa                                    >  2:685641/1‑150 (MQ=255)
                                                                                                                   gTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCa                                    <  1:4065740/150‑1 (MQ=255)
                                                                                                                    tGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAg                                   >  2:1134955/1‑150 (MQ=255)
                                                                                                                    tGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAg                                   >  2:3776538/1‑150 (MQ=255)
                                                                                                                      ggggCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTc                                 <  1:204349/150‑1 (MQ=255)
                                                                                                                       gggCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCt                                <  2:4702274/150‑1 (MQ=255)
                                                                                                                        ggCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCtt                               <  1:42812/150‑1 (MQ=255)
                                                                                                                         gCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCg                                                                                                          >  1:1827029/1‑74 (MQ=255)
                                                                                                                           gAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTct                            <  1:2712126/150‑2 (MQ=255)
                                                                                                                           gAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCa                            >  1:5011857/1‑150 (MQ=255)
                                                                                                                              gTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATgg                         <  2:2852002/150‑1 (MQ=255)
                                                                                                                               ttCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGtt                                                     <  2:1754258/121‑1 (MQ=255)
                                                                                                                               ttCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGc                        <  1:1857795/150‑1 (MQ=255)
                                                                                                                                tCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCt                       <  1:238517/150‑1 (MQ=255)
                                                                                                                                     tcggCACGCCCCCCAGCGGGGCGTGCCAATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAg                  >  1:2777652/3‑150 (MQ=255)
                                                                                                                                     aTGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAg                  >  2:3812500/1‑150 (MQ=255)
                                                                                                                                        gCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCAt               >  1:617452/1‑150 (MQ=255)
                                                                                                                                         cACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATg              <  1:2309987/150‑1 (MQ=255)
                                                                                                                                          aCGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGt             >  1:1919435/1‑150 (MQ=255)
                                                                                                                                          aCGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGt             >  2:1121539/1‑150 (MQ=255)
                                                                                                                                            gCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCg           >  2:2948269/1‑150 (MQ=255)
                                                                                                                                             ccccGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGt          <  1:4806417/150‑1 (MQ=255)
                                                                                                                                             ccccGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGt          >  1:3138144/1‑150 (MQ=255)
                                                                                                                                              cccGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTg         >  2:489433/1‑150 (MQ=255)
                                                                                                                                              cccGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTg         <  2:4067262/150‑1 (MQ=255)
                                                                                                                                              accGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTg         >  2:2450563/2‑150 (MQ=255)
                                                                                                                                               ccGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTgg        >  2:2280504/1‑150 (MQ=255)
                                                                                                                                                   tGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGGCCGt    >  1:3348350/1‑150 (MQ=255)
                                                                                                                                                   tGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGGCCGt    >  2:5494264/1‑150 (MQ=255)
                                                                                                                                                    ggggggCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGGCCGtt   <  1:5155270/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGGGAAAATCGTTGGAATGTTCTTGGCCATGTGCACGGAGCAGGTCTCAACGTGAAAGCGTGCGGCCAGCATCCGGCCCGCCACGCAGGCGATATTCCGGATTCGAGTGGAGAGGTGGGGGCGAAGTTCCGAATGGCACGCCCCGCTGGGGGGCGTGCCGATCAACGGCGGTCAGGCGCCGGTCGGGCGCAGCGAGTAAGTCTTGAGTTGTGCTGCGAAGTCGCGCAGTGCCTGGATGCCGCTGGTTTCGGCCTCATGTACCCAGTCTTTCATGGCTGCCAGCATGTCGTGGCCGTTG  >  NC_002947/268597‑268895

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: