Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,557,574 T→A 7.7% intergenic (‑80/+94) PP_4887 ← / ← PP_4888 hypothetical protein/methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,557,5740TA 7.7% 330.9 / 8.2 127intergenic (‑80/+94)PP_4887/PP_4888hypothetical protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (48/66);  new base A (8/3);  total (57/70)
Fisher's exact test for biased strand distribution p-value = 6.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

CGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGC  >  NC_002947/5557427‑5557722
                                                                                                                                                   |                                                                                                                                                     
  tACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCg                                                                                                                                                   >  2:5263744/1‑150 (MQ=255)
  tACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCg                                                                                                                                                   >  2:273862/1‑150 (MQ=255)
      cTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTTTGTGGCGCAAGGGGCCCATAGCGGccc                                                                                                                                               >  1:3331482/1‑150 (MQ=255)
      cTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGccc                                                                                                                                               >  1:2916402/1‑150 (MQ=255)
       tGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGcccc                                                                                                                                              >  1:903960/1‑150 (MQ=255)
         tCACTGTGCTGGATATTGAGGTCGGGCTGGGCAGCGATTAGTCGGGCGGACATGAGCATCTCCTTTTGTCTTGTGTCGCCGACAGCATAGCAAGAGCGCTCTGGTGTGGGAGCGGGCTTGTGACGCGTAAGGGCCGCTATGCGGcccttt                                                                                                                                            <  1:613082/150‑4 (MQ=255)
           aCTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCCCCTAGCGGCCCCGGcc                                                                                                                                          >  1:3135061/1‑149 (MQ=255)
            cTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCCTAGCGGCCCCCGCaa                                                                                                                                         >  1:1774017/1‑150 (MQ=255)
                gCAGGAAATTCAGGTCGGCCTGGGCAGCTATTAGTCGGGCGGAAATGAGAAACTCCTTTTGTCATGTGTCGCCTACAGAATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTc                                                                                                                                     <  1:3384967/150‑1 (MQ=38)
                 ctggAAATTCTGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGAAACTCCTTTTGTCTTGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGTAGCGGCCTTGTGTCGCTAAAGGGCCGCATAGCGGCCCCGGCAATTTcc                                                                                                                                    <  1:1641760/148‑1 (MQ=255)
                 cAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGGAAGGGCCGCATAGCGGCCCCGGCAATTTcc                                                                                                                                    >  1:4979245/1‑150 (MQ=255)
                   ggAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGGAAGGGCCCCATAGCGGCCCCGGCCAATTAcgc                                                                                                                                  >  2:2313984/1‑150 (MQ=39)
                   ggAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCcgc                                                                                                                                  <  1:4811252/150‑1 (MQ=255)
                    gAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCcg                                                                                                                                   >  1:1823428/1‑148 (MQ=255)
                    gAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCcg                                                                                                                                   <  2:1823428/148‑1 (MQ=255)
                       aTTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCcgc                                                                                                                                  <  2:2109219/146‑1 (MQ=255)
                         tCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAAGTTCCGCGTGTAt                                                                                                                            >  2:4860289/1‑150 (MQ=255)
                          cAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATg                                                                                                                           <  1:4905974/150‑1 (MQ=255)
                             gTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGcaa                                                                                                                        <  1:2145506/150‑1 (MQ=255)
                             gTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGcaa                                                                                                                        <  1:2145852/150‑1 (MQ=255)
                                    tGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCt                                                                                                                 <  2:2916402/150‑1 (MQ=255)
                                      ggCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTgg                                                                                                               <  1:4862322/150‑1 (MQ=255)
                                       gCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTggg                                                                                                              <  2:3571818/150‑1 (MQ=255)
                                          gCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGgct                                                                                                           <  1:751543/150‑1 (MQ=255)
                                           cGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGgctg                                                                                                          <  1:1758528/150‑1 (MQ=255)
                                           cGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGgctg                                                                                                          >  1:4626498/1‑150 (MQ=255)
                                            gATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGgctgc                                                                                                         >  2:5419660/1‑150 (MQ=255)
                                            gATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAAGTTCCGCGGGGTTGCAACAAACCAGGGgctgc                                                                                                         >  1:3131722/1‑150 (MQ=38)
                                             aTTAGTCGGGCGGAAATTAGAAACTCCTTTTTTCTTGTGTCGCCGACAGCATAGCAACACCGCTCTGGTGTGGGTGCGGCCTTGTGTCGCGAATGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGgctgct                                                                                                        <  1:2427009/150‑1 (MQ=38)
                                                 gTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGCAAGGGCCCCATAGCGGCCCCGCCCATTACCGCGGGTATGCCACAAAACTGGGGCTGCTATGc                                                                                                    >  2:4749641/1‑150 (MQ=38)
                                                 gTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGc                                                                                                    <  2:502237/150‑1 (MQ=255)
                                                 gTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGc                                                                                                    <  2:4626498/150‑1 (MQ=255)
                                                 gTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAAACTGGGGCTGCTATGc                                                                                                    >  2:2120219/1‑150 (MQ=255)
                                                  tCGGGGGGACATGCGCACCTCCTTTTGTCATGTGTCGCCGACAGCATCGCAACACCGCGCTGGTGTGGGAGCGGCCTTGTGTCGCGAACGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCa                                                                                                   <  1:1509162/150‑1 (MQ=38)
                                                  tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCCTAGCGGCCCCCGCAAATTCCGCGTGTATGCAACAAACCAGGGGCTGCGATGcc                                                                                                   >  2:4270856/1‑149 (MQ=39)
                                                  tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGc                                                                                                                                                    >  1:4038464/1‑101 (MQ=255)
                                                  tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGc                                                                                                                                                    <  2:4038464/101‑1 (MQ=255)
                                                  tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCa                                                                                                   <  2:4600921/150‑1 (MQ=255)
                                                  tCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCa                                                                                                   <  2:2504482/150‑1 (MQ=255)
                                                    gggCGGACATGAGCAACTCCTTTTGTCATGTGGCGCCGACAGCATCGCAACACCGCTCAGGGGGGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGc                                                                                                 <  1:4568205/150‑1 (MQ=255)
                                                       cGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCTGGTGTGGGAGCGGCCTTGTGTCGCTAAAGGGCCGCATCGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCATAGCAGcccc                                                                                              <  1:2769658/150‑2 (MQ=38)
                                                       cGGACATGAGAAACTCCTTTTGTCTTGTGTCGCCGACAGCATAGAAACACCGCTCCGGTGTGGGAGCGGCCTTGTGTCGCGAATGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCt                                                                                              <  1:936380/150‑1 (MQ=38)
                                                          aCATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTc                                                                                           <  2:903960/150‑1 (MQ=255)
                                                          aCATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTc                                                                                           <  1:546244/150‑1 (MQ=255)
                                                          aCATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTc                                                                                           <  2:335708/150‑1 (MQ=255)
                                                          aCATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCATGGGCTGCTATGCAGCCCATTc                                                                                           >  2:3852277/1‑150 (MQ=255)
                                                            aTGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTcgc                                                                                        <  1:4645597/150‑1 (MQ=255)
                                                              gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGCAAGGGCCCCCTAGCGGCCCCCGCCAATTACGCGGGGATGCAACAAAACAGGGGCTGCGATGCAGCCCATTCGCGa                                                                                      >  2:2959140/1‑150 (MQ=255)
                                                              gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGa                                                                                      <  1:1889608/150‑1 (MQ=255)
                                                              gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGa                                                                                      <  1:1889554/150‑1 (MQ=255)
                                                              gAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCCGCAATTTACGCGGGTATGCAACAAAACTGGGGCTGCTATGCAGCCCGTTCGCGa                                                                                      >  2:4189756/1‑150 (MQ=38)
                                                                gCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCACCACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGaca                                                                                    <  1:115427/150‑1 (MQ=255)
                                                                gCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGaca                                                                                    <  1:4065702/150‑1 (MQ=255)
                                                                  aaCTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAAATTCCGCGTGTATGCAACAAACCAGGGGCTGCTATGCAGCCCATTCGCGAcaca                                                                                  >  2:865718/1‑150 (MQ=255)
                                                                   aCTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGGAAAGGCCGCAAAGCGGCCCCCGCACGTTACGCCGGTAATCAACAAACCTGGGGCAGCTATTCACCCCATTAGCGACACaa                                                                                 >  1:4930718/1‑150 (MQ=25)
                                                                       cTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTGGCGCAATGGGCGCCTATCGGCCCCCGCCAATTCCGCGGGGTTGCAACAAACCTGGGGCTGCGATGCAGCCCTTTGGCGGCACAAgggg                                                                             >  2:1885546/1‑148 (MQ=255)
                                                                        ttttttCATGTGTCGCCGACAGAATCGCAACACCGCTCAGGTGTGGTTGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGccc                                                                            <  2:3801165/150‑1 (MQ=37)
                                                                        ttttGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCCTAGCGGCCCCGGCCATTTCCGCGTGTATGCAACAAACCAGGGGCTGCGATGCAGCCCCTTAGCGACACAAGGccc                                                                            >  2:3411004/1‑150 (MQ=38)
                                                                        ttttGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGccc                                                                            >  2:4481996/1‑150 (MQ=255)
                                                                            gTCTTGTGTCGCCGACAGCATAGCAACCCGGCTCTGGTGTGGGCGCGGCGTTGTGGCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACaa                                                                        <  2:8665/150‑1 (MQ=37)
                                                                               aTGGGTCGCCGACAGCATAGCAACACCGCTCTGGTGTGGGTGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAgg                                                                     <  1:1755533/150‑1 (MQ=255)
                                                                                    tCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTACGCGGGTATGCAACAAACCAGGGGCTGCGATGCAGCCCGTTCGCGACACAAAGCCCACAAAGGACcgc                                                                >  2:387701/1‑150 (MQ=38)
                                                                                        cGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAg                                                            <  2:573038/150‑1 (MQ=255)
                                                                                          aCAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGc                                                          <  1:2227983/150‑1 (MQ=255)
                                                                                            aGCATAGCAACACCGCTCTGGTGGGGGCGGGGCGTTGTGTCGCTAACGGGCCGCATCGCGGCCCCGGCATTTGCCGGGTGTATGCAACAAACCTGGGGCTGATATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCtt                                                        <  1:557860/150‑1 (MQ=255)
                                                                                              cagaGCACCACCGCTCTGGTGTGGGCGCGGCCTTGTGGCGCGAACGGGCCGCATAGCGGCCCCGGCAATTCCCGCGTGTTTGCAACAACCCTGGGGCTGCTATGCAGCCCTTCCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCa                                                      <  1:1492764/147‑1 (MQ=255)
                                                                                              cATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCa                                                      <  2:3616358/150‑1 (MQ=255)
                                                                                               atcccAACCCGGCTCCGGTGTGGGCGCGGCCTTGTGGCGCGAAGGGGCCGCATCGGGGCCCCGGAAATTTCCGGGTGTATGCAACAAACGTGGGGCTGCTATGCAGCCCTTCCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAg                                                     <  1:2314982/146‑1 (MQ=255)
                                                                                                 aGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAAGGCCCCAAAGCGGCCCCCGCACATTACGCGGGTTTGCAACAAAACTGGGGCGGCGATGCAGCCCGTTAGCGAAAAAAGGCCCACAAACGACCCCGCAGGCTTCACGc                                                   >  2:2991594/1‑150 (MQ=21)
                                                                                                   cAACACCGCTCAGGTGTGGGAGCGGCCTTGTGGGGCCGAAGGGGCCCATATCCGCCCCCGCCAGTTACCCGTGGAATCAACAAAACAGGGGCGGCGATTCAGCCCGGTAGCCCAACCAAGGCCCACAAAGACCCCGCAGGCGTCAGGccc                                                 >  2:1410809/1‑150 (MQ=255)
                                                                                                       acCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGAAATTTCCGCTTTTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGaa                                             <  2:2547165/150‑1 (MQ=255)
                                                                                                        ccccTCCGGTGTGGGAGGGGCCTTGTGTCGCGAACGGGCCGCATCGCGCCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAAc                                            <  1:4216208/147‑1 (MQ=38)
                                                                                                        ccGCTCTGGTGTGGGAGCGGCCTGGTGGCGCTAACGGGCCCCATTGCGGCCCCGGCAATTTCCGCGTTTATGCAAAAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAAc                                            <  1:2068172/150‑1 (MQ=38)
                                                                                                             cAGGTGTGGGAGCGGCCTTGTGTCGCGAACGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTTGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGcc                                       <  2:293577/150‑1 (MQ=255)
                                                                                                                    gggCGCGGCCTTGTGCCGCGAACGGGCCGCTTCGCGGCCCCGGCAATTGCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTGCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTg                                <  1:4270856/150‑1 (MQ=255)
                                                                                                                                gtgtCGCGAAAGGGCCGCATCGCGGCCCCGGCAATGTCCGCGTGTATGCAACAACCCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGa                    <  2:3135059/150‑1 (MQ=255)
                                                                                                                                gtgtCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGa                    <  2:3135061/150‑1 (MQ=255)
                                                                                                                                gtgtCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGa                    <  2:1660654/150‑1 (MQ=255)
                                                                                                                                   tCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgcc                 >  1:2123147/1‑150 (MQ=255)
                                                                                                                                   tCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgcc                 >  1:2982204/1‑150 (MQ=255)
                                                                                                                                     gcgAACGGGCCGCACCGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc               <  1:5011003/150‑1 (MQ=255)
                                                                                                                                     gcgAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc               >  2:268957/1‑150 (MQ=255)
                                                                                                                                     gcgAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc               >  1:5376501/1‑150 (MQ=255)
                                                                                                                                     gcgAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc               >  1:5376703/1‑150 (MQ=255)
                                                                                                                                     gcgAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAgccgc               >  2:1891488/1‑150 (MQ=255)
                                                                                                                                       gAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCa             <  2:4859799/150‑1 (MQ=255)
                                                                                                                                         aaGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGa           <  1:4860289/150‑1 (MQ=255)
                                                                                                                                          ggggCCGCTATGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGat          >  2:4216168/2‑150 (MQ=255)
                                                                                                                                          ggggCCGCTATGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGat          >  2:4215876/2‑150 (MQ=255)
                                                                                                                                            ggCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGa                                              <  1:5158267/112‑1 (MQ=255)
                                                                                                                                            ggCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGa                                              >  2:5158267/1‑112 (MQ=255)
                                                                                                                                            ggCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGAAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatca        >  2:164918/1‑150 (MQ=255)
                                                                                                                                             gCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcat       >  1:2256364/1‑150 (MQ=255)
                                                                                                                                             gCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcat       >  2:679902/1‑150 (MQ=255)
                                                                                                                                             gCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcat       >  2:2134787/1‑150 (MQ=255)
                                                                                                                                              ccGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGatcatc      <  1:2016346/150‑1 (MQ=255)
                                                                                                                                               cgaaTAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATcc     >  2:4863820/4‑150 (MQ=255)
                                                                                                                                                gCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCa    >  2:5237543/1‑150 (MQ=255)
                                                                                                                                                gCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCa    <  2:1774017/150‑1 (MQ=255)
                                                                                                                                                  aTAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGc  <  2:677921/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGTACGCTGTCACTGTGCAGGAAATTCAGGTCGGCCTGGGCAGCGATTAGTCGGGCGGACATGAGCAACTCCTTTTGTCATGTGTCGCCGACAGCATAGCAACACCGCTCAGGTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGCATAGCGGCCCCGGCAATTTCCGCGTGTATGCAACAAACCTGGGGCTGCTATGCAGCCCTTTCGCGACACAAGGCCCACAAAGGACCGCGCAGGCTTCAGGCCCGGAACTGGCCCAGGCTGGCGCGCAGCTGAGCCGCCAGATCATCCAGC  >  NC_002947/5557427‑5557722

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