Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,963,792 T→C 19.0% intergenic (+66/+76) 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,7920TC19.0% 624.7 / 127.0 277intergenic (+66/+76)PP_2589/PP_2590aldehyde dehydrogenase family protein/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (137/83);  new base C (19/33);  total (158/119)
Fisher's exact test for biased strand distribution p-value = 9.80e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.82e-01

Only 100 of 281 total aligned reads displayed.

TCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACAC  >  NC_002947/2963643‑2963941
                                                                                                                                                     |                                                                                                                                                     
tCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTACGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGt                                                                                                                                                       >  2:4729773/1‑150 (MQ=255)
  ggggTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCt                                                                                                                                                     >  2:4502213/1‑150 (MQ=255)
   gggTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTg                                                                                                                                                    >  2:1928/1‑150 (MQ=255)
    ggTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGa                                                                                                                                                   <  2:4205819/150‑1 (MQ=255)
    ggTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGa                                                                                                                                                   <  2:3332326/150‑1 (MQ=255)
      tttGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAgg                                                                                                                                                 <  2:946467/150‑1 (MQ=255)
         ggCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGtt                                                                                                                                              <  1:775163/150‑1 (MQ=255)
            cgcgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                                           <  2:2722802/150‑1 (MQ=255)
            cgcgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                                           <  1:5333628/150‑1 (MQ=255)
            cgcgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                                           >  2:121730/1‑150 (MQ=255)
            cgcgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                                           >  2:7406360/1‑150 (MQ=255)
            cgcgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                                           >  1:3736241/1‑150 (MQ=255)
             gcgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCATCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGc                                                                                                                                          <  1:5582660/150‑1 (MQ=255)
               gATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTg                                                                                                                                        <  2:4638150/150‑1 (MQ=255)
                aTCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAGGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTga                                                                                                                                       <  1:5630834/150‑2 (MQ=37)
                 tCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAACCCTCAGACTTGTGCTgaa                                                                                                                                      <  2:1770427/150‑3 (MQ=38)
                  cTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGAc                                                                                                                                     <  2:3922886/150‑1 (MQ=255)
                  cTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAAc                                                                                                                                     <  2:7377587/150‑1 (MQ=37)
                    ttCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCt                                                                                                                                   <  1:2022670/149‑1 (MQ=37)
                    gTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCt                                                                                                                                   <  2:6080356/150‑1 (MQ=37)
                     tCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtg                                                                                                                                  >  1:3361902/1‑150 (MQ=255)
                     tCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACAACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtg                                                                                                                                  >  1:5908610/1‑150 (MQ=255)
                      cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCtgt                                                                                                                                 <  1:3952528/150‑1 (MQ=37)
                      cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCtgt                                                                                                                                 <  1:1807508/150‑1 (MQ=37)
                      cACTGCATTCGTTAGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                                 <  2:1023699/150‑1 (MQ=255)
                       aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgtg                                                                                                                                <  1:6230911/150‑1 (MQ=255)
                       aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCtgtg                                                                                                                                <  1:1372883/150‑1 (MQ=37)
                        cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                               <  2:3247051/150‑1 (MQ=255)
                        cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCACGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTgg                                                                                                                               <  2:6713232/150‑1 (MQ=25)
                        cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCGGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACCGGTTCAGCACACGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                               >  1:2763594/1‑150 (MQ=255)
                          gCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGa                                                                                                                             <  1:3039207/150‑1 (MQ=255)
                              tcgttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCgg                                                                                                                         <  1:479829/150‑1 (MQ=25)
                              tcgttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCgg                                                                                                                         <  1:3493979/150‑1 (MQ=25)
                                    gACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGgcgcgc                                                                                                                   <  2:6996073/150‑1 (MQ=255)
                                      cacttATTCGAAACTCCCGCCCACCTGTTACCAGGTGCGCGGCTGATTTTTTTATCATTGACCGCGTGTGCTGCTTCGGGGGCATGCCGGCGCCCACCGGTTCCGCACACGTCTGAGGGTTGTGCTGGACGTGGGGGAGGGGGCGCGccc                                                                                                                 <  2:2525579/146‑1 (MQ=255)
                                      cAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTTATGTTTTGATCATTGACCGCGTCGTCTGATTCGCGGGCAGGCCCGCTCCCACAGGTTCAGCACAATTCTGAGGTTTGTGCTGGACCTGTGGGAGCGGGCGCGccc                                                                                                                 <  2:1455558/150‑1 (MQ=255)
                                           aTTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGGTCAGCACAAGTCTGAGGGTTGGGCTGGACCTGTGGGAGCGGGCGCGCCCGCGaa                                                                                                            >  2:6761136/1‑150 (MQ=255)
                                               gCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAgg                                                                                                        >  1:5185562/1‑150 (MQ=255)
                                               gCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAgg                                                                                                        >  1:6967767/1‑150 (MQ=255)
                                               gCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCGTCGCGGGCAGGCCCGCCCCCCCAGGGTCAGCCCAAGGCGGGGGGGGGGGCTGGGCCGGGGGGAGCGGGCGCCCCCCCGCagcag                                                                                                        >  1:1760683/1‑147 (MQ=255)
                                                    tCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTa                                                                                                   >  2:5297057/1‑150 (MQ=255)
                                                     cAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAAc                                                                                                                                     <  1:400115/115‑1 (MQ=25)
                                                     cAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAAc                                                                                                                                     >  2:400115/1‑115 (MQ=25)
                                                            accTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCCCCCGCGAAGAGGCCGGATATGgcgc                                                                                           >  1:4384199/1‑150 (MQ=255)
                                                              cTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGcgccg                                                                                         >  2:5471184/1‑150 (MQ=255)
                                                                 gTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcgg                                                                                      >  2:6876388/1‑150 (MQ=255)
                                                                  tccccgTTGCGGGGCGGATGTTTTGATCAGTGACCGCGTCGTCTGCTGCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcggc                                                                                     <  2:5874234/145‑1 (MQ=38)
                                                                  tcccAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGCCCGGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcggc                                                                                     <  2:4960564/148‑1 (MQ=255)
                                                                  tcccAGGTGCGCGGCTGATTTATTGATCATTGACCGGGTCGTCGGCTGCGCGGGCATGCCCGCTCCCCCCGGTCCCGCACACGTCTGCGGGTTGTGCTGGACCTGGGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcggc                                                                                     <  2:2124082/148‑1 (MQ=255)
                                                                   aCCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGgcgc                                                                                                                     >  1:1209146/1‑117 (MQ=255)
                                                                   aCCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGgcgc                                                                                                                     <  2:1209146/117‑1 (MQ=255)
                                                                   aCCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcggcg                                                                                    >  1:6655452/1‑150 (MQ=255)
                                                                     cAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCggg                                                                                  >  2:3003398/1‑150 (MQ=255)
                                                                                 tGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGATCTGGCGCCGCGGCGGGACCCGGTTAGaa                                                                      >  1:2150955/1‑150 (MQ=255)
                                                                                      tATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACACGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCCCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGGTAGAAGTccc                                                                 >  1:1240772/1‑149 (MQ=255)
                                                                                       aTTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGGGGGTGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCAc                                                                >  2:3823736/1‑150 (MQ=255)
                                                                                         tGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGc                                                              >  2:5554191/1‑150 (MQ=255)
                                                                                           aTCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTa                                                            >  2:5742060/1‑150 (MQ=255)
                                                                                           aTCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATATGGCGCCGCGGCGGGCCCCGGTTAGAAGTCCACGcaa                                                            >  1:4619824/1‑148 (MQ=255)
                                                                                              aTTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGcc                                                         >  2:7448786/1‑150 (MQ=255)
                                                                                                  aCCGCGTCGTCTGCTTCGGGGGCATGCCCGCTCCCACCGGTTCAGCACACGTCTGTGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGaaa                                                     <  1:7422170/150‑1 (MQ=255)
                                                                                                  aCCGCGTCGTCTGCTTCGGGGGAATGCCCGCTCCCACAGGTTCAGCACAAGTCTGTGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGaaa                                                     <  1:4940873/150‑1 (MQ=255)
                                                                                                   ccGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCCGCACACGTCTGTGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAg                                                    <  2:1651360/150‑1 (MQ=255)
                                                                                                         cgtcTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTGTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATa                                              <  2:2239902/150‑1 (MQ=255)
                                                                                                          gtcTGCCTCGCGGCCACGCCCGCTCCCACCGCTCCAGCACAAGTCTGCGGGTTGTGCCGGACCGGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAg                                             <  1:7315330/150‑1 (MQ=255)
                                                                                                                ttCGCGGGCGCGCCCGCTCCCACAGGTCCAGCACAACCCTCCGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGCCCCGGGTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCcgc                                       >  2:4302309/1‑150 (MQ=17)
                                                                                                                ttCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCcgc                                       >  2:3281548/1‑150 (MQ=255)
                                                                                                                     gggCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGgcgc                                  >  1:4428340/1‑150 (MQ=255)
                                                                                                                      ggCATGCCCGCTCCCACAGGTTCAGCACACGTGTGTGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGcc                                 <  2:52133/150‑1 (MQ=255)
                                                                                                                       gCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGGTGTGCTGGACCGGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGccc                                >  2:114007/1‑150 (MQ=255)
                                                                                                                         aTGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGccccc                              >  1:2695520/1‑150 (MQ=255)
                                                                                                                            cccGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATATGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAAc                           >  1:5071145/1‑150 (MQ=255)
                                                                                                                             ccGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAAcc                          >  1:6434738/1‑150 (MQ=255)
                                                                                                                                 tCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGGTAGAAGTCCACGCTAGCCCAAAGCGGATAGGTCCGCGGCGCCCCCAAcccccc                      >  1:6382782/1‑148 (MQ=25)
                                                                                                                                  cccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg                     <  2:6448745/150‑1 (MQ=255)
                                                                                                                                   ccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAggg                    >  2:535388/1‑150 (MQ=255)
                                                                                                                                    cacaGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGt                   >  1:250315/1‑150 (MQ=25)
                                                                                                                                     acaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTCGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTc                  >  2:1513418/1‑150 (MQ=255)
                                                                                                                                     acaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTc                  >  1:1053723/1‑150 (MQ=255)
                                                                                                                                     acaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTc                  >  1:320362/1‑150 (MQ=255)
                                                                                                                                     acaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTc                  >  2:409579/1‑150 (MQ=255)
                                                                                                                                     acaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTc                  >  2:4837830/1‑150 (MQ=255)
                                                                                                                                     acaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTc                  >  1:5224292/1‑150 (MQ=255)
                                                                                                                                      caGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCt                 >  1:6085506/1‑150 (MQ=255)
                                                                                                                                       aGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTg                >  2:3891106/1‑150 (MQ=255)
                                                                                                                                         gtccAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCAc                                                                >  2:6624427/4‑100 (MQ=18)
                                                                                                                                         gtccAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCAc                                                                <  1:6624427/97‑1 (MQ=18)
                                                                                                                                         gtccAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg                     >  1:2850746/4‑143 (MQ=37)
                                                                                                                                         gtccAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg                     <  2:2850746/140‑1 (MQ=37)
                                                                                                                                         gtccAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTggg              >  2:2552981/4‑150 (MQ=37)
                                                                                                                                          ttCAGCACAAGTCTGAGGGTTGTGCCGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGt             >  2:2394954/1‑150 (MQ=255)
                                                                                                                                           tCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTc            >  1:3292239/1‑150 (MQ=255)
                                                                                                                                           ccAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTc            >  2:4517603/2‑150 (MQ=38)
                                                                                                                                            cAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCAGGCGCCGCGGCGGGACCCGCGAAGAAGGCCACGCGAGCCCAAAGCTGATAGCTACGCGGCGCCCCCAACCCCCGGGTCTGGGTCa           >  1:1002722/1‑150 (MQ=38)
                                                                                                                                              ccacaCGTCTGCGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGt         <  2:5168536/149‑1 (MQ=255)
                                                                                                                                              ccacaCGTCTGAGGGTTGTGCTGGACCTGGGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGt         <  2:4396900/149‑1 (MQ=255)
                                                                                                                                               cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTg        >  2:4718674/1‑150 (MQ=255)
                                                                                                                                                acaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTgg       >  1:5474571/1‑150 (MQ=255)
                                                                                                                                                acaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTgg       >  2:5712091/1‑150 (MQ=255)
                                                                                                                                                     tCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCacac  >  1:6819480/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACAC  >  NC_002947/2963643‑2963941

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: