Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,963,799 G→C 19.5% intergenic (+73/+69) PP_2589 → / ← PP_2590 aldehyde dehydrogenase family protein/ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,963,7990GC19.5% 665.4 / 148.1 270intergenic (+73/+69)PP_2589/PP_2590aldehyde dehydrogenase family protein/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (142/73);  new base C (20/33);  total (162/108)
Fisher's exact test for biased strand distribution p-value = 2.60e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 280 total aligned reads displayed.

TTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTA  >  NC_002947/2963650‑2963948
                                                                                                                                                     |                                                                                                                                                     
ttGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAggg                                                                                                                                                       <  1:763717/150‑1 (MQ=255)
     cgcgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                                                  <  2:2722591/150‑1 (MQ=255)
     cgcgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                                                  >  2:7407642/1‑150 (MQ=255)
     cgcgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                                                  >  2:121730/1‑150 (MQ=255)
     cgcgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                                                  <  2:5184410/150‑1 (MQ=255)
       cgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCt                                                                                                                                                <  1:5300023/150‑1 (MQ=37)
          tCTGTCACTGCATTCGTTCGCCAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTgaa                                                                                                                                             <  2:4788083/150‑3 (MQ=37)
          tCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAACCCTCAGACTTGTGCTgaa                                                                                                                                             <  2:1770427/150‑3 (MQ=38)
          tCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTgaa                                                                                                                                             <  2:5477575/150‑3 (MQ=38)
           cTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGAc                                                                                                                                            <  2:3922917/150‑1 (MQ=255)
           cTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGAc                                                                                                                                            <  1:326074/150‑1 (MQ=255)
           cTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAAc                                                                                                                                            <  2:3142652/150‑1 (MQ=37)
             gTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCt                                                                                                                                          <  2:6080356/150‑1 (MQ=37)
             gTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCACGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCt                                                                                                                                          <  2:3829391/150‑1 (MQ=37)
              tCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtg                                                                                                                                         >  1:3361902/1‑150 (MQ=255)
               cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCtgt                                                                                                                                        <  1:2842192/150‑1 (MQ=25)
               cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                                        <  2:4687898/150‑1 (MQ=255)
               cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                                        <  1:679117/150‑1 (MQ=255)
               cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                                        <  2:227320/150‑1 (MQ=255)
               cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                                        <  2:1513118/150‑1 (MQ=255)
               cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                                        <  2:1513552/150‑1 (MQ=255)
                aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgtg                                                                                                                                       <  2:2398101/150‑1 (MQ=255)
                aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgtg                                                                                                                                       <  1:6230911/150‑1 (MQ=255)
                aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgtg                                                                                                                                       <  1:5757816/150‑1 (MQ=255)
                 cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCCGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                                      <  2:1189650/150‑1 (MQ=255)
                 cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTGTGAGGGTTGTGCTGGACCTGTgg                                                                                                                                      >  1:1894879/1‑150 (MQ=255)
                 cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                                      <  1:1418757/150‑1 (MQ=255)
                 cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCACGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTgg                                                                                                                                      <  2:6713232/150‑1 (MQ=25)
                   gCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGa                                                                                                                                    <  1:3039207/150‑1 (MQ=255)
                       tcgttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCgg                                                                                                                                <  1:3493979/150‑1 (MQ=25)
                         gttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACCGGTTCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGgc                                                                                                                              <  2:3666049/150‑1 (MQ=37)
                             gACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGgcatgc                                                                                                                          <  2:5321302/150‑5 (MQ=18)
                               catgTATTCGCAACTCACGACCACCTGGTACCTGTGGCGCGGCTGATGTTTTGATCATTGACCGCGTCGTCTGCTTCGGGGGAATGCCCGCACCCACAGGTTCAGCACAAGTCTGAGGGTTGGGCTGGACCGGTGGGAGCGGGCGCGccc                                                                                                                        <  2:2113935/147‑1 (MQ=255)
                               cAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTTATGTTTTGATCATTGACCGCGTCGTCTGATTCGCGGGCAGGCCCGCTCCCACAGGTTCAGCACAATTCTGAGGTTTGTGCTGGACCTGTGGGAGCGGGCGCGccc                                                                                                                        <  2:1455558/150‑1 (MQ=255)
                                aaGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCcg                                                                                                                       >  1:6115826/1‑150 (MQ=255)
                                    atgcGCAACTCACGCCCCCGGGGGACCCGGTGCGCGGCGGCGGTATTGATCATTGACCGCGTGGTCGGCTTCGCGGGCCTGCCCGCGCCCCCCGGTGCCGCCCCCGGCTGAGGGTTGTGCTGGACCGGGGGGAGCGGGCGCGCCCGCGaa                                                                                                                   <  2:5259601/147‑1 (MQ=255)
                                    aTTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATTTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCTCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGaa                                                                                                                   <  1:3928161/150‑1 (MQ=17)
                                    aTTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGGTCAGCACAAGTCTGAGGGTTGGGCTGGACCTGTGGGAGCGGGCGCGCCCGCGaa                                                                                                                   >  2:6761136/1‑150 (MQ=255)
                                        gCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAgg                                                                                                               >  1:6967767/1‑150 (MQ=255)
                                             tCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                                        <  2:3301497/120‑1 (MQ=255)
                                             tCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCGGTGGGAGCGGGCGCGCCCGCGAAGAGGCCgga                                                                                                          >  1:5058192/1‑148 (MQ=255)
                                              cAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAAc                                                                                                                                            <  1:400115/115‑1 (MQ=25)
                                               aaGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATc                                                                                                        >  1:2980285/1‑150 (MQ=255)
                                                aGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCt                                                                                                       >  1:6499384/1‑150 (MQ=255)
                                                 gACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTg                                                                                                      >  2:5713216/1‑150 (MQ=255)
                                                     accTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCCCCCGCGAAGAGGCCGGATATGgcgc                                                                                                  >  1:4384199/1‑150 (MQ=255)
                                                          gTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcgg                                                                                             >  2:6876388/1‑150 (MQ=255)
                                                          gTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcgg                                                                                             >  2:7044685/1‑150 (MQ=255)
                                                           tccccgTTGCGGGGCGGATGTTTTGATCAGTGACCGCGTCGTCTGCTGCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcggc                                                                                            <  2:5874234/145‑1 (MQ=38)
                                                           tcccAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGCCCGGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcggc                                                                                            <  2:4960564/148‑1 (MQ=255)
                                                            aCCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGgcgc                                                                                                                            >  1:1209551/1‑117 (MQ=255)
                                                            aCCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcggcg                                                                                           >  1:6655120/1‑150 (MQ=255)
                                                              cAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCggg                                                                                         >  2:3003398/1‑150 (MQ=255)
                                                                 gtGCGCGGCTTATGTTTTGATCTTTGACCGCGTCGTGTGCTTCGCGGGCATGCCCGCGCCCACAGGTTCCGCACACGTCTGCGGGTTGTGCTGGACCTGTGGGCGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGAcc                                                                                      <  2:819806/149‑1 (MQ=38)
                                                                               tATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGATGGTTGTGCTGGACCTGTGGGCGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAATTccc                                                                        >  1:909091/1‑149 (MQ=255)
                                                                                    aTCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGGTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGGTAGAAGTCCACGCTa                                                                   >  2:2544005/1‑150 (MQ=255)
                                                                                           cccGGGTCGGCTGGTGCGGGGGCTGGCCCGCTCCCACGGGTTCCGCAACAGTCTGTGGGTGGTGGTGGACCTGTGGGCGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGaaa                                                            <  2:2763594/149‑1 (MQ=255)
                                                                                           aCCGCGTCGTCTGCTTCGGGGGAATGCCCGCTCCCACAGGTTCAGCACAAGTCTGTGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGaaa                                                            <  1:4940873/150‑1 (MQ=255)
                                                                                            ccGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCCGCACACGTCTGTGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAg                                                           <  2:1651360/150‑1 (MQ=255)
                                                                                               cgtcgtCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTg                                                        >  1:5478562/1‑150 (MQ=255)
                                                                                                   gtcTGCCTCGCGGCCACGCCCGCTCCCACCGCTCCAGCACAAGTCTGCGGGTTGTGCCGGACCGGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAg                                                    <  1:7315330/150‑1 (MQ=255)
                                                                                                         ttCGCGGGCGCGCCCGCTCCCACAGGTCCAGCACAACCCTCCGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGCCCCGGGTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCcgc                                              >  2:4302309/1‑150 (MQ=17)
                                                                                                         ttCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGGGCGCCCGCGAAGAGGCCGTATATGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTAcgc                                              >  1:4442349/1‑150 (MQ=255)
                                                                                                               ggCAGGCCCGCTCCCCCAGGTTCAGCACACGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGcc                                        <  2:4619824/150‑1 (MQ=255)
                                                                                                                  aTGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGccccc                                     >  2:4001431/1‑150 (MQ=255)
                                                                                                                  aTGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGccccc                                     >  1:2695520/1‑150 (MQ=255)
                                                                                                                           cccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGGGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg                            <  2:6464531/150‑1 (MQ=255)
                                                                                                                            ccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAggg                           >  2:1045998/1‑150 (MQ=255)
                                                                                                                            ccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAggg                           >  1:6430578/1‑150 (MQ=255)
                                                                                                                            ccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAggg                           >  1:1328951/1‑150 (MQ=255)
                                                                                                                            ccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAggg                           >  2:2207657/1‑150 (MQ=255)
                                                                                                                            ccACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAggg                           >  1:2515042/1‑150 (MQ=37)
                                                                                                                             cacaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGt                          >  2:565033/1‑150 (MQ=255)
                                                                                                                              acaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTc                         >  1:320362/1‑150 (MQ=255)
                                                                                                                              acaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTc                         >  1:5224278/1‑150 (MQ=255)
                                                                                                                              acaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTc                         >  1:1053723/1‑150 (MQ=255)
                                                                                                                              acaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCAGATAGGTCCGCGGCGCCCCCAACCCCAGGGTc                         >  2:5363635/1‑150 (MQ=255)
                                                                                                                                 ggtaccgCACAGGGCTGCGGGTTGTGGTGGACCGGTGGGCGGGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTgg                      <  1:4729773/144‑1 (MQ=255)
                                                                                                                                  gtccAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCAc                                                                       >  2:6624427/4‑100 (MQ=18)
                                                                                                                                  gtccAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTggg                     >  1:5107484/4‑150 (MQ=37)
                                                                                                                                  gTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCACGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTggg                     >  2:325455/1‑150 (MQ=255)
                                                                                                                                   ttCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGt                    >  2:4242716/1‑150 (MQ=255)
                                                                                                                                    tCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTc                   >  2:739809/1‑150 (MQ=255)
                                                                                                                                    ccAGCACAACACTCAGACTTGTGCTGAACCTGTGGGAGCGGTCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTc                   >  2:4517341/2‑150 (MQ=38)
                                                                                                                                     cAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCAGGCGCCGCGGCGGGACCCGCGAAGAAGGCCACGCGAGCCCAAAGCTGATAGCTACGCGGCGCCCCCAACCCCCGGGTCTGGGTCa                  >  1:1002722/1‑150 (MQ=38)
                                                                                                                                       gCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGt                >  1:1787053/1‑150 (MQ=255)
                                                                                                                                       ccacaCGTCTGAGGGTTGTGCTGGACCTGGGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGt                <  2:4396900/149‑1 (MQ=255)
                                                                                                                                        cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATCGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTg               >  2:4718728/1‑150 (MQ=255)
                                                                                                                                         acaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTgg              >  1:5474447/1‑150 (MQ=255)
                                                                                                                                         acaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTACGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTgg              >  1:2083337/1‑150 (MQ=255)
                                                                                                                                           aaGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTAAGTGGca            >  1:5590480/1‑150 (MQ=255)
                                                                                                                                              tCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCacac         >  1:6819480/1‑150 (MQ=255)
                                                                                                                                              tCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCacac         >  1:3975706/1‑150 (MQ=255)
                                                                                                                                                tGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTgg                      >  1:2307820/1‑135 (MQ=255)
                                                                                                                                                tGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACAccc       >  1:5915651/1‑150 (MQ=255)
                                                                                                                                                tGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACAccc       >  2:5747372/1‑150 (MQ=255)
                                                                                                                                                   gggTTGGGGTGGACCTGGGGGCGGGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCag    <  1:946186/150‑1 (MQ=255)
                                                                                                                                                    ggTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCagt   >  1:1685672/1‑150 (MQ=255)
                                                                                                                                                    ggTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCagt   >  2:5163889/1‑150 (MQ=255)
                                                                                                                                                    ggTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGGCAGTGGCACACCCCagt   >  2:7073642/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTA  >  NC_002947/2963650‑2963948

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: