Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,503,385 T→A 13.5% intergenic (+26/‑193) PP_4834 → / → PP_4835 SPFH domain‑containing protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,503,3850TA13.5% 330.8 / 3.5 127intergenic (+26/‑193)PP_4834/PP_4835SPFH domain‑containing protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/69);  new base A (16/2);  total (54/73)
Fisher's exact test for biased strand distribution p-value = 8.09e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

GGGCGAACGCTATGTGGGCGCGATGGAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATA  >  NC_002947/5503239‑5503533
                                                                                                                                                  |                                                                                                                                                    
gggCGAACGCTATGTGGGCGCGATGGAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATC                                                                                                                                                   <  1:3021311‑M1/150‑1 (MQ=255)
 ggCGAACGCTATGTGGGCGCGATGGAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCG                                                                                                                                                  >  2:2004589‑M1/1‑150 (MQ=255)
  gCGAACGCTATGTGGGCGCGATGGAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCGTTGCGCCCCAATCGC                                                                                                                                                 >  2:1560220‑M1/1‑150 (MQ=255)
    gAACGCTATGTGGGCGCGATGGAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGA                                                                                                                                               <  2:4657410‑M1/150‑1 (MQ=255)
         cTATGTGGGCGCGATGGAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAA                                                                                                                                          <  2:1289290‑M1/150‑1 (MQ=255)
            tgtgGGCGCGATGGAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGC                                                                                                                                       <  1:4514949‑M1/150‑1 (MQ=255)
             gtgGGCGCGATGGAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCC                                                                                                                                      <  2:1577193‑M1/150‑1 (MQ=255)
              tgGGCGCGATGGAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGccc                                                                                                                                     <  1:1701395/150‑1 (MQ=255)
               gggCGCGATGGAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGcccc                                                                                                                                    <  1:3972364/150‑1 (MQ=255)
                  cgcgATGGAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCcaca                                                                                                                                 <  2:4701506/150‑1 (MQ=255)
                     gATGGAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGa                                                                                                                              <  2:2784795/150‑1 (MQ=255)
                     gATGGAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGa                                                                                                                              <  2:2784825/150‑1 (MQ=255)
                        ggAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAg                                                                                                                           >  2:4317642/1‑150 (MQ=255)
                        ggAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCGGTGCGCCCCAATCGCGACACAAAGCCCCACACGACAg                                                                                                                           >  2:4317705/1‑150 (MQ=255)
                         gagaACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCGATCGCGCCACAAGGCCCCCCAGGACAgc                                                                                                                          >  2:563463/1‑150 (MQ=255)
                            aaCCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCa                                                                                                                       <  1:1708351/150‑1 (MQ=255)
                            aaCCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCAAACGCCACACAAAGCCCCACAGGACAGCGCa                                                                                                                       >  2:2483484/1‑150 (MQ=255)
                              ccTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTTCGCCCCTATAGCGACACAAGGCCCCACAGGACACCGCAGt                                                                                                                     >  2:4775868/1‑150 (MQ=255)
                                tCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTAt                                                                                                                   >  2:2396673/1‑150 (MQ=255)
                                 cgcgGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTTGGCCCCGAAAGCGACACAAGGCCCCCCACGCCACCGCAGTAtt                                                                                                                  >  1:3686362/1‑150 (MQ=255)
                                    ggggAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGa                                                                                                               <  1:5239325/150‑1 (MQ=255)
                                    ggggAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGa                                                                                                               <  2:2907005/150‑1 (MQ=255)
                                    ggggAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGa                                                                                                               <  1:4806168/150‑1 (MQ=255)
                                    ccgatctagcaATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGa                                                                                                               <  1:949818/143‑1 (MQ=255)
                                      ggAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAgg                                                                                                             <  2:1159154/150‑1 (MQ=255)
                                       gAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGGGGTTTTGGGCCCTTTGGGGACACAATGGGCCGCAGGGCAGGGGGGTAttggttgg                                                                                                            >  1:2063232/1‑144 (MQ=255)
                                          cagcaATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCGTTGCGCCCCAATCGCGACACAAGGCCCCACACGACAGCGCAGTATTAGAGGCTGc                                                                                                         >  1:392341/1‑150 (MQ=255)
                                            gcaATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGGGCGTTTGGCCCCCATCGCGACACAAAGCCCCACACGCCAGCGCCCTATTAGGGGCTGACg                                                                                                       >  1:4610942/1‑150 (MQ=255)
                                            gcaATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTTCGCCCCAATAGCGACACAAGGCCCCACACGACAGCGCCCTATTagagagatagg                                                                                                       >  2:949818/1‑143 (MQ=255)
                                               aTGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCaa                                                                                                    <  1:658535/150‑1 (MQ=255)
                                                tGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAAt                                                                                                   <  1:5297313/150‑1 (MQ=255)
                                                tGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAAt                                                                                                   <  1:5297340/150‑1 (MQ=255)
                                                 gCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCGTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAAtt                                                                                                  >  2:5116582/1‑150 (MQ=255)
                                                   cAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCGTTGCGCCCCAATAGCGACACAAGGCCCCCCAGGACAGCGCCGTATTAGAGGCATCCGCaaatag                                                                                                >  2:1070510/1‑146 (MQ=255)
                                                     aGGTGGTGGTGTTGCCGGCTGACCGGCAAGAGACCGTGCGCGGGGTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCtgt                                                                                              <  2:165818/150‑1 (MQ=255)
                                                       gtggtggtgTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTgg                                                                                            <  1:4006808/150‑1 (MQ=255)
                                                       gtggtggtgTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTgg                                                                                            <  2:2938184/150‑1 (MQ=255)
                                                       gtggtggtgTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCGGTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTgg                                                                                            <  2:5160538/150‑1 (MQ=255)
                                                           tggtgTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCGTTGCGCCCCTATCGCCACACAAGGCCCCACAGGACAGCGCCGTATTCGAGGCTGCCGCAAATCAGTGGGAGc                                                                                        >  2:3923558/1‑150 (MQ=255)
                                                                 tGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCggcctt                                                                                  >  1:3900506‑M2/1‑144 (MQ=255)
                                                                  gCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAg                                                                                         <  1:3079680/142‑1 (MQ=255)
                                                                  gCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAg                                                                                         >  2:3079680/1‑142 (MQ=255)
                                                                   ccGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCggccttgt                                                                                >  2:1806790‑M2/1‑142 (MQ=255)
                                                                     ggCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCgctcccacagaattgcggcagcctcgaatactgcgctgtcctgtggggccttgtgt                                                                              <  1:4369904‑M1/150‑57 (MQ=255)
                                                                      gCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCggccttgtgtc                                                                             >  2:3880206‑M2/1‑139 (MQ=255)
                                                                       cTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCCATCGCGACACAAGGCCgctcccacagaattgcggcagcctcgaatactgcgctgtcctgtggggccttgtgtcg                                                                            <  2:4802756‑M1/150‑59 (MQ=255)
                                                                        tGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCgctcccacagaattgcggcagcctcgaatactgcgctgtcctgtggggccttgtgtcgc                                                                           <  1:2197044‑M1/150‑60 (MQ=255)
                                                                        tGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCggccttgtgtcgc                                                                           >  2:3457443‑M2/1‑137 (MQ=255)
                                                                        tGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCGTTGCCCCCCCAAACCCCCACCACGCCCCACACGACACCCCCCGATTCGCGGCGGCCCCAAATCAGTGGGCGCCGCCGGGTGTcgc                                                                           >  2:2835296/1‑150 (MQ=255)
                                                                           ccTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTTCGCCCCAATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCCATTCAGTTGGAGCgggcaggtgtcgcgat                                                                        <  2:3722364‑M2/150‑17 (MQ=255)
                                                                            cTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGGGCTTTTGGCCCCAATAGCGTCACAAGGCCCCACAGGAGAGCGCAGTATTGGTGGCTGCCGCCATTATTTTGGTTGGGGCTTGTGTGGCGATg                                                                       >  1:5280786/1‑150 (MQ=255)
                                                                            cTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATTGCGACACAAGGCCgctcccacagaattgcggcagcctcgaatactgcgctgtcctgtggggccttgtgtcgcgata                                                                       <  2:4874421‑M1/150‑64 (MQ=255)
                                                                             tGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTTCGCCCCAATAGCGACACAAGGCCCCACAGGACAGCGCCGTATTAGAGGCTGCCGCAATTATGTGGGAGCGGCCTTTTGTCGCGATgg                                                                      >  1:4030805/1‑150 (MQ=255)
                                                                             tGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCGTTGGGCCCCAAACGCGCCACAAAGCCCCACCGGACAGCCCACTCTTAGAGGCTGCCGCAAATAATTGGGAGGGGCCGTGTGTggcgcatg                                                                      >  1:599327/1‑146 (MQ=255)
                                                                             tGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCGTTGCGCCCCTATCGCGACACAAGGCCCCACAGGCCAGCGCAGTATTCGAGGCTGCCGCAATTCGGTGGGAGCggccgggtgtcgcgaggg                                                                      <  2:2522560‑M2/150‑19 (MQ=255)
                                                                              gCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCAAAACCGACACAATGCCCCACACGCCCGCGCCGTATTCGTGGCGGCCGCACATCAGTGGGAGCGGCCGGGTGTCGCCGTggg                                                                     >  2:3430445/1‑150 (MQ=255)
                                                                              gCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCGTTGCGCCCCTATCGCGCCCCAACGCCCCACAGGCCCGCGCAGTATTCGAGGCGGCCGCACGTCTGTGGGAGCGGGCGTGTGTCGCGATggg                                                                     >  2:3430508/1‑150 (MQ=255)
                                                                                aaGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCgctcccacagaattgcggcagcctcgaatactgcgctgtcctgtggggccttgtgtcgcgatag                                                                      >  1:3280443‑M1/1‑83 (MQ=255)
                                                                                aaGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCgctcccacagaattgcggcagcctcgaatactgcgctgtcctgtggggccttgtgtcgcgatag                                                                      <  2:3280443‑M1/147‑65 (MQ=255)
                                                                                aaGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCgctcccacagaattgcggcagcctcgaatactgcgctgtcctgtggggccttgtgtcgcgatag                                                                      >  1:3280537‑M1/1‑83 (MQ=255)
                                                                                     aCCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGCCCCAAGGCCCCACAGGACCGCGCAGTATTCGAGGCGGCCGCCATTCTGTGGGAGCGGCCGTGTGTCGCGATGGGCTGcccc                                                              >  2:5358221/1‑147 (MQ=255)
                                                                                     aCCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGcaca                                                              <  2:4582092/150‑1 (MQ=255)
                                                                                       cGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGTGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGctgtgc                                                            <  1:666626/150‑6 (MQ=255)
                                                                                        gTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAcagca                                                           <  1:1905079/150‑1 (MQ=255)
                                                                                        gTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAcagca                                                           <  1:1905326/150‑1 (MQ=255)
                                                                                            gcgGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTTGGCCCCAATAGCGACACAAAGCCCCACAGGACACCGCCGTAATCGAGGCTGCCGCAAATATGTTGGATCGGCCGTTTTTCGCGATGGGGTGCACAGCAGccc                                                       >  1:1268869/1‑150 (MQ=255)
                                                                                              gggTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGGGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAAATCTGTGGGAGCGGCCGTGTGTCGCGCTGGGCTGCACAGCAGccccc                                                     >  2:1187901/1‑150 (MQ=255)
                                                                                                gTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCTGTGCAGcccatcg                                                   <  2:4610942/150‑5 (MQ=255)
                                                                                                 ttGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCTGTGCAGcccatcgc                                                  <  2:2949631/150‑6 (MQ=255)
                                                                                                 ttGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGAGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCAtt                                                  <  1:400741/150‑1 (MQ=255)
                                                                                                      gggCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTTCGCCCCTATAGCGACAAAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAAGTCTGTTGGGGCGGCCTTGTTTGGCGATTGGGTGCACAGCAGCCCCCATTACCAg                                             >  2:78976/1‑150 (MQ=255)
                                                                                                      gggCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTTCGCCCCTATAGCCACACAAAGCCCCACAGGAAAGCGCAGTAATCGAGGCTGCCCCAAATCGGTGGGAGCGGCCGTGTGTAGCGATTGGCTGCACAGCAGCCCCCAATaaaac                                             >  2:3456991/1‑146 (MQ=255)
                                                                                                      gggCCGCAACAAGATTTGAGAGTGTCGGGGGCGCGTTTCGCCCCAATAGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAAATCAGTGGGGGCGGCCTAGTGTCGCGATGGGCTGCACAGCAGCCCCCAAAACCAg                                             >  2:3456956/1‑150 (MQ=255)
                                                                                                        gCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGac                                           <  2:3243910/150‑1 (MQ=255)
                                                                                                        gCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGac                                           <  2:3243819/150‑1 (MQ=255)
                                                                                                             aacaaGATTTGAGAGTGTCGGGGGCGCTTTGGGCCCCAATAGCGGCACAAGGCCCCACAGGGCAGCGCAGTATTAGAGGCTGCCGCAAATCTTTTGGTGCGGCCGGGTGTCGCGATGGGCTGCACAGCAGCCCCCAATAACACaccaccg                                      >  2:4144949/1‑149 (MQ=255)
                                                                                                               caaGATTTGAGAGTGTCGGGGGGGGTTTTGGCCCCTATGGCGTGACAATGGCCCACAGGTCATCGCAGTTTTGGTGGGTGGGGGAATTGTGTGGGTGCGGGCTTGTGTGGGGTTTGGCTGGACAGCGGGCCCCATTAACAGtacaccagg                                    >  2:404001/1‑147 (MQ=255)
                                                                                                               caaGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGc                                    <  2:878318/150‑1 (MQ=255)
                                                                                                                    tttGTGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACACGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCAc                               <  1:1560220/150‑1 (MQ=255)
                                                                                                                              tCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCAc                     <  2:1789675/150‑1 (MQ=255)
                                                                                                                              tCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCAc                     <  2:1789858/150‑1 (MQ=255)
                                                                                                                                        tttGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCggccttgtgtcgcgatgggggcgcaaagcgcccccgacactctcaaatcttgttgcggcccatcaacccgcgcacgg           <  2:253397‑M2/150‑78 (MQ=255)
                                                                                                                                        tttGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCAttt           >  1:4238368/1‑150 (MQ=255)
                                                                                                                                          tGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCAttttt         <  2:4238368/150‑1 (MQ=255)
                                                                                                                                          tGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCAttttt         <  1:5116582/150‑1 (MQ=255)
                                                                                                                                           gcgcCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTTGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTa        >  1:2193541/1‑150 (MQ=255)
                                                                                                                                            cgcCCCTATCGCTACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGAAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTAc       <  1:2429648/150‑1 (MQ=255)
                                                                                                                                            cgcCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTAc       <  1:1429529/150‑1 (MQ=255)
                                                                                                                                             gcCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTAcc      >  1:2263211/1‑150 (MQ=255)
                                                                                                                                             gcCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTAcc      <  1:4317705/150‑1 (MQ=255)
                                                                                                                                               cccTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTTGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCta    <  2:2193541/150‑1 (MQ=255)
                                                                                                                                               cccTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCt                                      >  1:130396/1‑116 (MQ=255)
                                                                                                                                                ccTCTCGCGCCCCACGGCCCCACAGGACAGCGCCGTATTCGAGGCCGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACCGCACCCCCCACGACCAGCCCACCAGCGTCACGCCCCCGCACCCCGGCATTATTACCtat   <  2:302161/150‑1 (MQ=255)
                                                                                                                                                ccTCTCGCGACACAAGGCCCCACTGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGAGATGGGCTGCACAGAAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCtat   <  2:302238/150‑1 (MQ=255)
                                                                                                                                                 cTATCGCGACACAAGGCCCCACAGGACAGCGCAGTTTTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCtata  <  1:2305011/150‑1 (MQ=255)
                                                                                                                                                 cTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCggccttgtgtcgcgataggggcgcaaagcgcccccgacactctcaaatcttgttgcggcccatcaacccgcgcacggtctctt     >  1:3280537‑M2/1‑64 (MQ=255)
                                                                                                                                                 cTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCggccttgtgtcgcgataggggcgcaaagcgcccccgacactctcaaatcttgttgcggcccatcaacccgcgcacggtctctt     <  2:3280443‑M2/147‑84 (MQ=255)
                                                                                                                                                 cTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCggccttgtgtcgcgataggggcgcaaagcgcccccgacactctcaaatcttgttgcggcccatcaacccgcgcacggtctctt     <  2:3280537‑M2/147‑84 (MQ=255)
                                                                                                                                                 cTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCtata  >  1:3958400/1‑150 (MQ=255)
                                                                                                                                                 cTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCtata  <  1:2396673/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGGCGAACGCTATGTGGGCGCGATGGAGAACCTCGCGGGGAGCAGCAATGCCAAGGTGGTGGTGTTGCCGGCTGACCTGCAAGAGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATA  >  NC_002947/5503239‑5503533

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: