Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,065,150 C→G 11.1% intergenic (+42/‑165) PP_1838 → / → PP_1839 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,065,1500CG11.1% 370.3 / 31.0 117intergenic (+42/‑165)PP_1838/PP_1839membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/70);  new base G (6/7);  total (40/77)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 123 total aligned reads displayed.

CCCGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACT  >  NC_002947/2065003‑2065299
                                                                                                                                                   |                                                                                                                                                     
cccGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCca                                                                                                                                                     >  1:668234/1‑150 (MQ=255)
 ccGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCGCGccccc                                                                                                                                                    >  1:4071673/1‑148 (MQ=255)
  cGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTACACCGGGGCCGCGTGGCGGCCCATTCGCGGGCTTGCCCGCTCCcaca                                                                                                                                                   <  2:4464645/150‑1 (MQ=255)
       cAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAGACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAt                                                                                                                                              <  1:2599517/150‑1 (MQ=255)
        aGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGTACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAtt                                                                                                                                             <  1:3096025/150‑1 (MQ=255)
        aGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCAtt                                                                                                                                             <  2:1602965/150‑1 (MQ=255)
         gCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTTCCCGCTCCCACAGGCAATg                                                                                                                                            >  1:3673930/1‑150 (MQ=255)
         gCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTg                                                                                                                                            <  2:2896876/150‑1 (MQ=255)
         gCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCGCCCACAGGCATTg                                                                                                                                            >  2:3178122/1‑150 (MQ=255)
         gCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCAGTCGCGGGCGGTCCCGCGCCCCCCGGccagg                                                                                                                                            >  1:2881224/1‑146 (MQ=255)
          ccGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGc                                                                                                                                           <  2:4235964/150‑1 (MQ=255)
           cGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGca                                                                                                                                          <  2:638246/150‑1 (MQ=255)
           cGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGca                                                                                                                                          <  1:2781239/150‑1 (MQ=255)
             gAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGcaca                                                                                                                                        <  1:606311/150‑1 (MQ=255)
             gAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGcaca                                                                                                                                        <  1:3139858/150‑1 (MQ=255)
             gAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGcaca                                                                                                                                        <  1:3139659/150‑1 (MQ=255)
               aaTGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGGGTGCCCGCGGCCCCAGGCCGTGCACAgg                                                                                                                                      >  1:1313135/1‑150 (MQ=255)
                   cTCAAGAATCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTc                                                                                                                                  <  1:231402/150‑1 (MQ=255)
                    tCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATGCGCGGGCTTGCCCGCTCCCCCAGGCCATGCCCAGGGTGCg                                                                                                                                 >  2:3461746/1‑150 (MQ=255)
                     cAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCAGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTAGCGGGCTTGCCCGCTCCCACAGGCCTTGCACAGGGTTCgg                                                                                                                                >  2:1244168/1‑150 (MQ=255)
                      aagaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACCGGCATTGCACAGGGTTCGGa                                                                                                                               >  2:1686361/1‑150 (MQ=255)
                       agaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGAc                                                                                                                              <  2:3087883/150‑1 (MQ=255)
                       agaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGAc                                                                                                                              <  1:3292524/150‑1 (MQ=255)
                        gaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGAcc                                                                                                                             <  1:2012485/150‑1 (MQ=255)
                        gaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGAACAGGGTTCGGAcc                                                                                                                             <  1:3885700/150‑1 (MQ=255)
                        gaagCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCGCCCACAGGCATTGCACAGGGTTCGGAcc                                                                                                                             >  2:1648856/1‑150 (MQ=255)
                          agCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCAGTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                           <  2:3887947/150‑1 (MQ=255)
                           gCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtg                                                                                                                          <  1:3098175/150‑1 (MQ=255)
                            cTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgt                                                                                                                         <  1:2535216/150‑1 (MQ=255)
                             tGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgtg                                                                                                                        <  2:524063/150‑1 (MQ=255)
                              ggAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGc                                                                                                                       <  2:2833701/150‑1 (MQ=255)
                              ggAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGc                                                                                                                       <  1:2802589/150‑1 (MQ=255)
                                   cGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGACCTGCACAAGGTCCGAACCCTGTGCaatgc                                                                                                                  <  2:854720/150‑5 (MQ=25)
                                   cGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGACCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTc                                                                                                                  <  1:3739410/150‑1 (MQ=255)
                                   cGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGACCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTc                                                                                                                  <  1:1845967/150‑1 (MQ=255)
                                     tctatCCAGACCGCGACCAAGAGCCTGCTGTATAACAGCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGTCCTTGTGCAGGTCCt                                                                                                                <  2:4070691/148‑1 (MQ=255)
                                       tatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                              <  2:668234/150‑1 (MQ=255)
                                        atCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTACCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgtg                                                                                                             <  1:2004052/150‑1 (MQ=255)
                                           cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGa                                                                                                          <  1:4308978/150‑1 (MQ=255)
                                            aGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAg                                                                                                         >  2:1442299/1‑150 (MQ=255)
                                             gACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGc                                                                                                        <  1:3910117/150‑1 (MQ=18)
                                             gACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGc                                                                                                        <  1:3876886/150‑1 (MQ=255)
                                              aCCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGTGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCg                                                                                                       <  2:598725/150‑1 (MQ=18)
                                               ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                      <  2:260766/150‑1 (MQ=18)
                                               ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                      <  2:2231412/150‑1 (MQ=18)
                                               ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                      <  1:991392/150‑1 (MQ=255)
                                               ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                      <  2:3807883/150‑1 (MQ=255)
                                               ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCgg                                                                                                      <  1:3136471/150‑1 (MQ=255)
                                               ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGCGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                      <  1:2392078/150‑1 (MQ=18)
                                                 gcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                      >  2:1593078/1‑148 (MQ=18)
                                                 gcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                      <  1:1593078/148‑1 (MQ=18)
                                                   gACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCaa                                                                                                  <  1:2324650/150‑1 (MQ=255)
                                                    aCCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAg                                                                                                 <  2:3410016/150‑1 (MQ=255)
                                                       aaGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGccc                                                                                              <  2:3477474/150‑1 (MQ=255)
                                                             ctgctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGGGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAAt                                                                                        <  2:413055/150‑1 (MQ=255)
                                                               gctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCCCCGGGTTCGGACCGTGTGCAGGGCCTGGGCGAGCGGGCCCGCCCGCGAATgg                                                                                      >  2:2552715/1‑150 (MQ=38)
                                                                ctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCCCAGGCCTTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCCAGCCCGCGAATggg                                                                                     >  2:2891354/1‑150 (MQ=255)
                                                                     tAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGGCCCCTCCCCCAGGCAATTCAAACGGTTAGGGCCAATTGCACGTACGGGGCGAGCCGGCCAGCCCCCCAAAGGGccccc                                                                                >  2:1625479/1‑149 (MQ=255)
                                                                       acacCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCCTTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGca                                                                              >  2:3676676/1‑150 (MQ=255)
                                                                         acCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCCCTCCCACAGGCCTTGCACAGGGTTCGGACCTTGTGCAGGGCCAGGGCGAGCGGGCAAGCCCGCGACTGGGCCCCGcaca                                                                            >  2:3241569/1‑150 (MQ=255)
                                                                         acCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCCTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGcaca                                                                            >  1:4080805/1‑150 (MQ=255)
                                                                           cgccgGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAAt                                                                          <  2:4070118/150‑1 (MQ=255)
                                                                               ggCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCCTTGCACAGGGTTCGGGCCATGTGCAGGGCCGGTGCGAGCGGGCAACCCCGCGAAAGGGCCCCGCACAAATTTc                                                                      >  2:3042903/1‑150 (MQ=37)
                                                                                     gtgCACGGCACACCGACGTAAACCGGGGCCGCGTTGCGGCCCATACGCGGGCTGGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCTAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGAc                                                                <  2:2968168/150‑1 (MQ=255)
                                                                                           ggCACACCTTAGTAACCGGGGGCGTGGTGGGGGCCCTTTGGCGGGGTTGCCGGCTCCCCCAGGCTTTGCAACGTGTTGGTACGTTGTGCAGTTCCTGTGGGAGGGGGCAAGCCCGCGAATGGGCCCCGCTCAATTTTCCCGGACTGAACa                                                          <  2:875423/150‑1 (MQ=255)
                                                                                           ggCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACCGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACa                                                          <  2:2619620/150‑1 (MQ=255)
                                                                                             cacaCCGAAGTAACCCGGGGCCGCGTTGCGGCCCATCCGCGGGCTGGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACAtt                                                        <  1:3327062/150‑1 (MQ=255)
                                                                                              acacCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACTGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTa                                                       <  2:4471362/150‑1 (MQ=255)
                                                                                                acCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCCCAGGCCTTGCACCGGGTTAGGACCGTTTGCAGGTCCTGTTCGAGCGGGCAAGCCCCCCAATGGGCCCCGCACAATTTTACCCGACAGAAccatata                                                     >  2:3022113/1‑144 (MQ=25)
                                                                                                acCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATc                                                     <  2:2714708/150‑1 (MQ=255)
                                                                                                     attAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCTAGGGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGc                                                <  2:4159429/148‑1 (MQ=255)
                                                                                                       tAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCCCTCCCACAGGCCTTGCACAGGGTTCGGTCCTTGTGCAGGTCCTGGGCGAGCGGGCAACCCCCCCAATGGGCCCCGCACAATTTTACCGGACTGAAAAGTATACCCGccc                                              >  2:4289029/1‑150 (MQ=255)
                                                                                                         aaCCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgtg                                                                                                                        >  1:3156935/1‑74 (MQ=255)
                                                                                                          ccgggggCCGCGTGGCGGCCCATTCGCGGGCTTGCCCGCTCCCACTGGCTTTTCACTGGGTTCGGACATTGTGCTGGTCGTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCaa                                           <  1:3679358/147‑1 (MQ=38)
                                                                                                           ccGGGGCCGCGTTGCGGCCCATTCGCGGGCGTGCCCCCTCCCCCAGGCCATGGCCAGGGTTCGGACCTTGTGCAGGTGCTGTGCGGGCGGGCAACCCCCCGGAGGGGCCCCCCCCAAATTTACCCGCCGGCACAATATCCCCCCCCCCAg                                          >  2:3392619/1‑150 (MQ=255)
                                                                                                           ccGGGGCCGCGTTGCGGCCCAGTCGCGGGCGTGCCCGCGCCCCCCGGCCGTGCCCCGGGGTCGGCCCGGGTGCCGGGGCGGGGCGCGCGGGCCAGCCCGCGCAGGGGCCCCGCACCATTTTCCCGGCCGGGACAGTATCCCCGccccccg                                          >  2:1662153/1‑147 (MQ=255)
                                                                                                                     gTTGCGGCCCATTCGCGGGCTTGCCCCCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCCCACAATTTTCCCGGACTGAACAATATCCCCCCCCCAAGGGCAGTTAtt                                >  2:1856021/1‑150 (MQ=38)
                                                                                                                      gggCGGCCCTTGCGCGGGGTGGCCCGCCCCCCCCGGCATTGCACTGGGGGGGGACGTTGGGCCGGTCCGGGGCGCGCGGCCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTc                               <  1:1298828/148‑1 (MQ=255)
                                                                                                                       tGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                              <  2:3099534/150‑1 (MQ=255)
                                                                                                                       tGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                              <  2:2425889/150‑1 (MQ=255)
                                                                                                                       tGCGGCCCATTCGCGGGCTTGCCCGCCCCCACAGGCATTGCACAGGGTTCGGACCTGGTGCAGGTCCTGGGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                              <  1:3364499/150‑1 (MQ=255)
                                                                                                                                      ggCTTGCCCGCTCCCACCGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTg               <  1:3010997/150‑1 (MQ=255)
                                                                                                                                      ggCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTTTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTg               <  1:2970857/150‑1 (MQ=255)
                                                                                                                                         ttGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgtg            >  1:1417826/1‑150 (MQ=255)
                                                                                                                                         ttGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgtg            >  2:4350571/1‑150 (MQ=255)
                                                                                                                                             ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg        >  1:1178477/1‑150 (MQ=18)
                                                                                                                                             ccGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg        >  1:564106/1‑150 (MQ=18)
                                                                                                                                             ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCTGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCg        >  1:1950586/1‑150 (MQ=255)
                                                                                                                                             ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGTCATCACACACCTGGTGTTct        >  1:3180912/1‑149 (MQ=255)
                                                                                                                                              cGCTCCCACAGGCATTGCACCGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGt       <  1:3352910/150‑1 (MQ=255)
                                                                                                                                               gCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGtt      >  2:905752/1‑150 (MQ=18)
                                                                                                                                               gCTCCCCCCGGCATTGCACCGGGTTCGGACCTTGTGCAGGTCGTGTGCGTGGGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCCGGCATCACACACCTGGTGTTCGtt      <  2:1978080/150‑1 (MQ=255)
                                                                                                                                               gCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGc                           <  1:4021593/129‑1 (MQ=255)
                                                                                                                                               gCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGc                           >  2:4021593/1‑129 (MQ=255)
                                                                                                                                                cTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTg     >  2:2026872/1‑150 (MQ=18)
                                                                                                                                                 tCCCACCGGCATTGCACTGGGTTCGGACCTTGTGCTGGTCCTGGGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGa    <  2:968924/150‑1 (MQ=255)
                                                                                                                                                   ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt  >  2:3159882/1‑150 (MQ=255)
                                                                                                                                                   ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt  >  2:3159917/1‑150 (MQ=255)
                                                                                                                                                   ccACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATAGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt  >  1:2223846/1‑150 (MQ=255)
                                                                                                                                                   ccACACGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGAGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt  >  2:3356824/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCCGGACCAGCCGGAAATGCTCAAGAAGCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTAAACCGGGGCCGCGTTGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACT  >  NC_002947/2065003‑2065299

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