Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,076,188 A→T 6.2% intergenic (+148/+79) PP_1846 → / ← PP_1847 group II intron‑encoding maturase/TonB‑dependent ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,076,1880AT 6.2% 547.6 / 4.3 175intergenic (+148/+79)PP_1846/PP_1847group II intron‑encoding maturase/TonB‑dependent ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base A (73/88);  new base T (3/11);  total (76/99)
Fisher's exact test for biased strand distribution p-value = 9.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 182 total aligned reads displayed.

GAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGT  >  NC_002947/2076039‑2076337
                                                                                                                                                     |                                                                                                                                                      
gAACCGCCGTGTACGGAACCGTACGCACGGTGATGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGa                                                                                                                                                        <  2:5209312/150‑1 (MQ=255)
 aaCCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGAt                                                                                                                                                       <  2:4673525/150‑1 (MQ=255)
 aaCCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGAt                                                                                                                                                       >  1:2165993/1‑150 (MQ=255)
  aCCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATg                                                                                                                                                      <  2:604187/150‑1 (MQ=255)
    cgccgTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGc                                                                                                                                                    <  2:2087977/150‑1 (MQ=255)
    cgccgTGTACGGAACCGTACGCACGGTGGGGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGc                                                                                                                                                    >  1:4077631/1‑150 (MQ=255)
          gtACGGAACCGTACGCACGCTGGAGAGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCCTTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCATCACCGTGGTGGATGGCACGGGc                                                                                                                                              <  1:2303907/150‑1 (MQ=37)
           tACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCt                                                                                                                                             >  2:995113/1‑150 (MQ=255)
           tACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCt                                                                                                                                             >  2:4834680/1‑150 (MQ=255)
             cGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTc                                                                                                                                           <  1:5543074/150‑1 (MQ=255)
             cGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTc                                                                                                                                           <  2:4179756/150‑1 (MQ=255)
               gAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGc                                                                                                                                         <  2:5507207/150‑1 (MQ=255)
               gAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGc                                                                                                                                         <  1:679494/150‑1 (MQ=255)
                 aCCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGccc                                                                                                                                       <  2:310663/150‑1 (MQ=255)
                 aCCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGccc                                                                                                                                       <  2:4416163/150‑1 (MQ=255)
                 aCCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGccc                                                                                                                                       >  1:3238794/1‑150 (MQ=255)
                      aCGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCGAAGCCCgtgcc                                                                                                                                  <  1:2239926/150‑3 (MQ=255)
                       cGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCgtgtca                                                                                                                                 <  1:3069254/150‑3 (MQ=255)
                          aCGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTcgcg                                                                                                                              >  2:2548760/1‑150 (MQ=255)
                           gtaaagattgtagtcactggatgccgtttgctttttttgggagcaactgtcttgctcaatttctaacagctggcgcgatccctgggggagcgggCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCgg                                                                                                                             <  2:2974239‑M2/56‑1 (MQ=255)
                           gtaaagattgtagtcactggatgccgtttgctttttttgggagcaactgtcttgctcaatttctaaaagctggcgcgatccctgggggagcgggCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCgg                                                                                                                             <  2:1757265‑M2/56‑1 (MQ=255)
                           gtaaagattgtagtcactggatgccgtttgctttttttgggagcaactgtcttgctcaatttctaaaagctgacgcgacccgtgtgggagcggcCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCgg                                                                                                                             <  2:3133391‑M2/56‑1 (MQ=255)
                           cGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACCCCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCgg                                                                                                                             <  2:4726031/150‑1 (MQ=255)
                           cGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCgg                                                                                                                             <  1:5094484/150‑1 (MQ=255)
                            ggtggtGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCggg                                                                                                                            <  2:1165599/150‑1 (MQ=255)
                               ggtgTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCaa                                                                                                                        <  2:387597/150‑1 (MQ=255)
                                ggcaccgcgtaccgctcggctttcggcttcgatcagcttctagcttgccggggccgctctgcggccctttcgcgacctgtggcAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAg                                                                                                                       >  1:5379836‑M1/84‑150 (MQ=255)
                                   tgtagtcactggatgccgtttgctttttttgggagcaactgtcttgctcaatttctaaaagctggcgcgacccgtgtgggagcggcCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGccc                                                                                                                    <  1:4023417‑M2/64‑1 (MQ=255)
                                    gagaGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCg                                                                                                                   >  2:4785587/1‑150 (MQ=255)
                                    gagaGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCg                                                                                                                   <  1:4194005/150‑1 (MQ=255)
                                     agagGACGGCGTGTGTGAATCCGCCTCCTACTCGTTTTTCCGGGCCTTTGGGGCTGCTGTGCAGCCCTTCGCGACAGACGCGTGCTTGGCGGCGAAGACAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCTAGCCCGc                                                                                                                  <  2:2695734/150‑1 (MQ=25)
                                     agagGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTGCCCGCGAACACGGGCGAAGCCCgt                                                                                                                  <  2:5046752/150‑2 (MQ=38)
                                       agGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCtt                                                                                                                <  1:4713029/150‑2 (MQ=255)
                                         cacgggatgccgtttgctttttttgggagcaactgtcttgctcaatttctaaaagctggcgcgatccctgggggagcggcCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTccc                                                                                                              <  2:439017‑M2/70‑1 (MQ=255)
                                             gcggGTGTTAACCCGCCTCCTACTCGTTTTGCTGGGCCATGGGGGCTGCTGTGCAGCCCACCGCGACACACGGGGGCTTGCCCGCGAAGCCAACCCCGGGGTGGAGGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAg                                                                                                          <  2:2356523/150‑1 (MQ=21)
                                                  tgtgAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTcgc                                                                                                     >  1:2870291/1‑150 (MQ=255)
                                                    tttcggcttcgatcagcttctagattgccggggccgctctgcggccctttcgcgacctgtgggAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTcgcg                                                                                                                              <  1:2762771‑M1/61‑1 (MQ=255)
                                                    tgAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGa                                                                                                   >  2:3369905/1‑150 (MQ=255)
                                                      aaTCCGCCTCCTACTCTTTTTTCCGGGCCATGGGGGCTGCTGTGCAGCCCTGCGCGACACACGGGGGCTTGCCCGCGAAGCCAACACCGTGGTGGTTGGCAGGGGCTTCGCCCGTGTTCGCGGGCATGCCCGCTCCCACAGGTCGCGaaa                                                                                                 <  1:1720122/150‑1 (MQ=21)
                                                       aTCCGCCTCCTCCTCGATTTTCCGGGCCTTTGGGGCTGCTGGGCAGCCCATCGCGACACACGCGGGCTTGCCCGCGAAGCAAACACCGGGGTGGTGGGCACGGGCTCCGCCGGTGTTCGCGGGAACGCCCGCTCCCACAGGTCGCGAAAg                                                                                                <  1:387659/150‑1 (MQ=25)
                                                       aTCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGACGCCAACCCCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAg                                                                                                <  2:2794625/150‑1 (MQ=255)
                                                        tCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGCGCAAGCCCGCTCCCACAGGTCGCGAAAgg                                                                                               <  1:4514252/150‑1 (MQ=255)
                                                            cctcctACTCTATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACATGCGGGCTTGCCCGTGAAGCCAACACCGTGGTGGTTGGCACGGGCTTGGCCCGTGTTCGCGGGAAAGCCCGCTCCCACAGGTCGCGAAAGGGCCg                                                                                           <  2:2513095/150‑1 (MQ=255)
                                                            cctcctACTCTATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGGTTGCCCGCGAAGCCAACACCGTGGTGGATGGCAGGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCg                                                                                           <  2:3067536/150‑1 (MQ=255)
                                                            cctcctACTCGATTTTCCGGGCCATGGGGGCTGCTGTGCAGCCCATCGCGACACACGCGGGCTTGCCCGCGAAGCCAACCCCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCg                                                                                           <  2:2197924/150‑1 (MQ=255)
                                                              tcctACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCa                                                                                         >  1:2870999/1‑150 (MQ=255)
                                                                ctACTCGATTTTCTGGGCCTTTGGGGCTGCTGTGCAGCCCATCGCGACACAAGGGGGCTTGCCCGCGAAGCCAACACGGTGGTGGATGGCACGGGGTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCaga                                                                                       <  2:3853315/150‑1 (MQ=37)
                                                                 tACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCagag                                                                                      <  2:2722547/150‑1 (MQ=255)
                                                                   cTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGTCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCg                                                                                    >  2:5599553/1‑150 (MQ=255)
                                                                    tCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCg                                                                                                                                          <  1:4611388/96‑1 (MQ=255)
                                                                    tCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGaa                                                                                                  <  2:5485600/135‑1 (MQ=255)
                                                                      gATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGcc                                                                                 >  2:5292540/1‑150 (MQ=255)
                                                                        ttttCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACATGCGGGCTTGCCCGCGAAGCCAACCCCGTGGTGTATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGcccc                                                                               <  2:135711/150‑1 (MQ=255)
                                                                        ttttCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACACGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCTCTGCGGccct                                                                               <  2:5143576/150‑2 (MQ=255)
                                                                           tCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGc                                                                            >  1:5561355/1‑150 (MQ=255)
                                                                               ggCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATc                                                                        >  1:2650921/1‑150 (MQ=255)
                                                                               ggCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATc                                                                        >  2:4878260/1‑150 (MQ=255)
                                                                                gCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCt                                                                       >  2:5499080/1‑150 (MQ=255)
                                                                                  ctttGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGACGCCACCCCCGTGGGGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAg                                                                     <  2:4488975/148‑1 (MQ=255)
                                                                                  cATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAg                                                                     >  1:1240695/1‑150 (MQ=255)
                                                                                  cATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAg                                                                     >  2:1710248/1‑150 (MQ=255)
                                                                                   aTTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCACTCTAGa                                                                    >  1:104658/1‑150 (MQ=255)
                                                                                   aTTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGa                                                                    >  1:2622297/1‑150 (MQ=255)
                                                                                    tgggggCTGCTGTGCAGCCCATCGCGACACACGGGGGCTTGCCGGCGAAGCCAACCCCTTGGTGTTTGGCAGGGGTTTCGCCCGTGTTCGCGGGCACGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGaa                                                                   <  1:2251501/148‑1 (MQ=37)
                                                                                    tgggggCTGCGGGGCGGCCCATCGCGACACACGGGGGCTGGCCCGCGACGCCAACCCGGTGGTGGATGGCAGGGGCTTCGCCCGTGTCCGCGGGCACGCCCGCTCCCCCAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGaa                                                                   <  1:5400042/148‑1 (MQ=255)
                                                                                         gctgctGTGCAGCCCATCGCGACACACGGGGGCTTGCCCGCGACGCCAACACCGTGGTGTTTGGCCCGGGCTTCGCCCGTGTTCGCGGGCACGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGa                                                              <  1:3732997/150‑1 (MQ=37)
                                                                                                 gCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaagcc                                                      >  1:4822240/1‑150 (MQ=255)
                                                                                                     cccATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaagccgaaa                                                  >  2:3448881/1‑150 (MQ=255)
                                                                                                      ccATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGa                                                                                                                                                        >  1:1256583/1‑48 (MQ=255)
                                                                                                      ccATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGa                                                                                                                                                        <  2:1256583/48‑1 (MQ=255)
                                                                                                        aTCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaagccgaaagcc                                               >  2:1058268/1‑150 (MQ=255)
                                                                                                        aTCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaagccgaaagcc                                               >  1:1784520/1‑150 (MQ=255)
                                                                                                          cgctccacaAGCGGGCTTGCCCGCGAAGCAACCACGGTGGTGGTTGGCTGGGGCTTCGCCCTTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaagccgaaagccga                                             <  1:5016434/145‑1 (MQ=37)
                                                                                                          cgccACACATGCGGGCTGGCCCGCCAACCCAACACCTTGGTGGATGCCACGGGCTGCCCCCGGGTGCGCGGCCACGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGaagccgaaagccga                                             <  1:5016348/150‑1 (MQ=255)
                                                                                                           gcgACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAg                                            >  1:1995260/1‑150 (MQ=255)
                                                                                                            cgACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGc                                           >  2:4758232/1‑150 (MQ=255)
                                                                                                             gaaacaAGCGGGCTTGCCCGCTAAGCAAAAACCGTTTTGTATGGCACGGGCTTCGCCCGTTTTCGCGGCAAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCg                                          <  1:1777520/147‑1 (MQ=255)
                                                                                                              acacAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGc                                                                                                                      <  1:5486262/73‑1 (MQ=255)
                                                                                                                      GGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGt                                 >  2:4601912‑M1/1‑150 (MQ=255)
                                                                                                                          TTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGccc                             <  2:4613958‑M1/150‑1 (MQ=255)
                                                                                                                             CCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCt                          <  1:3496896‑M1/150‑1 (MQ=255)
                                                                                                                               cgctacgCCAACCGCGTGGTGTTTGGCTCGGGGTCCGCCGTTGTTGGCGGGCACGCCCGCTCCCACCGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTgg                        <  1:4296585/144‑1 (MQ=255)
                                                                                                                               CGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTgg                        <  1:892715‑M1/150‑1 (MQ=255)
                                                                                                                               CGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTgg                        <  1:5335748‑M1/150‑1 (MQ=255)
                                                                                                                                GCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGa                       <  1:839077‑M1/150‑1 (MQ=255)
                                                                                                                                  GAAGCCAACCCCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAAc                     >  2:3276498‑M1/1‑150 (MQ=255)
                                                                                                                                   AAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACa                    <  1:5302242‑M1/150‑1 (MQ=255)
                                                                                                                                     gCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAgg                  <  1:3280613/150‑1 (MQ=255)
                                                                                                                                       caccccggtggtggtgGGCAGGGGCTCCGCCGGGTTGCGGGGGAACCCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTa                <  1:4220130/143‑1 (MQ=255)
                                                                                                                                       cAACCCCGTGGTGGATGGCAGGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTa                <  2:4020393/150‑1 (MQ=255)
                                                                                                                                       cAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTa                >  1:1508538/1‑150 (MQ=255)
                                                                                                                                        aaCACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCACTCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAg               >  2:4341611/1‑150 (MQ=255)
                                                                                                                                        aaCACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGCAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAg               >  2:5091827/1‑150 (MQ=255)
                                                                                                                                        aaCACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACTGGTAg               >  1:2461644/1‑150 (MQ=255)
                                                                                                                                          cacCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGtt             >  1:563591/1‑150 (MQ=255)
                                                                                                                                           cccGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTa            <  2:1690351/149‑1 (MQ=255)
                                                                                                                                             cGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATc          >  1:5428424/1‑150 (MQ=255)
                                                                                                                                                    gATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCCgg   <  2:5563506/150‑1 (MQ=255)
                                                                                                                                                    gATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAgg   <  2:2742398/150‑1 (MQ=255)
                                                                                                                                                     aTGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGt  <  2:2165993/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTTTCAGGGCCATTGGGGCTGCTGTGCAGCCCATCGCGACACAAGCGGGCTTGCCCGCGAAGCCAACACCGTGGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCGCTCCCACAGGTCGCGAAAGGGCCGCAGAGCGGCCCCGGCAATCTAGAAGCTGATCGAAGCCGAAAGCCGAGCGGTACGCGGTGCCCCCTGGAACAGGTAGTTATCCCCCAGGT  >  NC_002947/2076039‑2076337

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: