Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,960,755 T→C 16.1% intergenic (+56/+86) 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,960,7550TC16.1% 267.5 / 33.0 110intergenic (+56/+86)PP_2589/PP_2590aldehyde dehydrogenase family protein/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (44/46);  new base C (8/9);  total (52/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

GCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTC  >  NC_002947/2960608‑2960904
                                                                                                                                                   |                                                                                                                                                     
gCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCa                                                                                                                                                     <  2:1458158/150‑1 (MQ=255)
 cAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGGTCCAg                                                                                                                                                    <  2:482097/150‑1 (MQ=255)
   aGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGca                                                                                                                                                  <  1:4631970/150‑1 (MQ=255)
      aGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACaa                                                                                                                                               <  2:1709831/150‑1 (MQ=255)
      aGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACaa                                                                                                                                               <  2:4453711/150‑1 (MQ=255)
      aGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACaa                                                                                                                                               <  2:4453582/150‑1 (MQ=255)
      aGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACaa                                                                                                                                               <  1:2863875/150‑1 (MQ=255)
       gTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAg                                                                                                                                              >  2:841880/1‑150 (MQ=255)
            ggTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGa                                                                                                                                         <  1:592143/150‑1 (MQ=255)
                tGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAg                                                                                                                                                    >  1:4018997/1‑135 (MQ=255)
                tGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAg                                                                                                                                                    <  2:4018997/135‑1 (MQ=255)
                  gCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtg                                                                                                                                   >  2:3256228/1‑150 (MQ=255)
                  gCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtg                                                                                                                                   >  2:3256295/1‑150 (MQ=255)
                  gCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtg                                                                                                                                   <  1:1154765/150‑1 (MQ=255)
                  gCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtg                                                                                                                                   <  2:885182/150‑1 (MQ=255)
                   ccGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgt                                                                                                                                  <  2:3662274/150‑1 (MQ=255)
                     gcgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAACCCTCAGACTTGTGc                                                                                                                                <  1:3290692/150‑1 (MQ=39)
                     gcgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGc                                                                                                                                <  2:2893495/150‑1 (MQ=38)
                     gcgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGGTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGc                                                                                                                                <  1:827764/150‑1 (MQ=255)
                      cgATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCt                                                                                                                               <  1:2134254/150‑1 (MQ=255)
                         tCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGa                                                                                                                            <  1:3699319/150‑1 (MQ=255)
                         tCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGa                                                                                                                            >  1:1961618/1‑150 (MQ=255)
                         tCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTgaa                                                                                                                            <  2:4130833/150‑3 (MQ=38)
                          cTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGAc                                                                                                                           <  1:4960190/150‑1 (MQ=255)
                          cTGTCACTGAATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCAGAAGTCTGAGGGTTGTGCTGGAc                                                                                                                           <  1:3749168/150‑1 (MQ=255)
                              cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                       <  2:1342235/150‑1 (MQ=255)
                              cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                       <  1:5008051/150‑1 (MQ=255)
                              cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAACCCTCAGACTTGTGCTGAACCtgt                                                                                                                       <  1:502808/150‑1 (MQ=38)
                               ccTGCCTTCGTTCGACACGTCGGCGCAACTCCAGCCCCCCTGGTACCAGGGGCGCGGCTGATGTATGGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGGTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgtg                                                                                                                      <  2:4130155/149‑1 (MQ=255)
                               aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGa                                                                                                                                         <  2:1398488/131‑1 (MQ=255)
                               aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGa                                                                                                                                         <  2:210820/131‑1 (MQ=255)
                               aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGa                                                                                                                                         >  1:1398488/1‑131 (MQ=255)
                               aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGa                                                                                                                                         >  1:210820/1‑131 (MQ=255)
                               aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgtg                                                                                                                      <  2:387144/150‑1 (MQ=255)
                               aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgtg                                                                                                                      <  2:2711973/150‑1 (MQ=255)
                               aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCggtg                                                                                                                      >  1:1838249/1‑150 (MQ=255)
                                cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                     <  1:4356147/150‑1 (MQ=255)
                                cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                     <  2:4975158/150‑1 (MQ=255)
                                     ttcgttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAg                                                                                                                                        >  2:1036038/1‑126 (MQ=255)
                                     ttcgttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAg                                                                                                                                        <  1:1036038/126‑1 (MQ=255)
                                     ttcgttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAg                                                                                                                                        <  1:1036172/126‑1 (MQ=255)
                                     ttcgttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCg                                                                                                                >  1:2728998/1‑150 (MQ=255)
                                        gttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCGGAGGGTTGTGCTGGACCTGTGGGAGCGGgc                                                                                                             >  1:4835858/1‑150 (MQ=255)
                                        gttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATGGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCCCAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGgc                                                                                                             <  2:2100146/150‑1 (MQ=255)
                                              cAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGccc                                                                                                       <  2:1277638/150‑1 (MQ=18)
                                               aaGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCcg                                                                                                      >  1:1876145/1‑150 (MQ=255)
                                                  tATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGa                                                                                                   >  1:2730993/1‑150 (MQ=255)
                                                    ttCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGc                                                                                                                                >  1:460607/1‑119 (MQ=37)
                                                    ttCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGc                                                                                                                                <  2:460607/119‑1 (MQ=37)
                                                        cAACTCAAGACCACCTGGTACCCGTTGCGCGGCTGATGTATTTATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGc                                                                                             <  2:1200023/150‑1 (MQ=255)
                                                            tCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGCAGGGGCCGTa                                                                                         >  1:1013814/1‑150 (MQ=255)
                                                              aaGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGc                                                                                                                                <  2:4372442/109‑1 (MQ=255)
                                                              aaGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGc                                                                                                                                >  1:4372442/1‑109 (MQ=255)
                                                               cgACCACCGGGTACCAGGTGCGCGGCTGAGGTATGGATCTTTGACCGCGTCGTCGGCTGCGCGGGCAGGCCCGCTCCCCCCGGTTCCGCACACGGCTGAGGGTGGGGCTGGCCCTGGGGGCGCGGGCGCGCCCGCGAAGAGGCCGTATCt                                                                                      <  2:2989272/149‑1 (MQ=255)
                                                               acaccaccTGTTACCTGTTGCGCGGCTGATTTTTTGATCTTTGACCGCGTCGTCTGGTCCGCGGGCATGCCCGCTCCCACCGGTACCGCACAAGTCTGCGGGTTGTGGTGGACCTGTGGGAGGGGGCGCGCCCGCGAAGAGGCCGTATCt                                                                                      <  2:1876145/148‑1 (MQ=25)
                                                                gACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTg                                                                                     >  2:1424502/1‑150 (MQ=255)
                                                                gACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTgaa                                                                                                                            >  2:1064735/1‑109 (MQ=21)
                                                                gACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTgaa                                                                                                                            <  1:1064735/111‑3 (MQ=21)
                                                                      cTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGcgccg                                                                               >  2:4978900/1‑150 (MQ=255)
                                                                               gTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCAGGGGGAGCGGGCGCGCCCGCGCACAGGCCGTATCATGCGCCCCGGCgggcc                                                                      >  1:1342116/1‑148 (MQ=38)
                                                                                    gcgGCTGATTTATTGATCATTGACGGCGTCTTCTGCTTCGGGGGCATGCCCGCTCCCACCGGTTCCGCACACGTCTGTGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGt                                                                 <  2:1631552/150‑1 (MQ=38)
                                                                                    gcgGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACCGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGt                                                                 <  2:4376350/150‑1 (MQ=255)
                                                                                    gcgGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGt                                                                 >  1:2610356/1‑150 (MQ=255)
                                                                                          gATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCGGTGGGAGCGGGCGCGCCCGCGAAGAGGCCCTATCTGGCGCCGCGGCGGGACCCCGTTAGAAg                                                           >  1:2331325/1‑150 (MQ=255)
                                                                                              tATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACCGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGgc                                                                                                             <  1:2389369/96‑1 (MQ=255)
                                                                                              tATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCa                                                       >  1:4832553/1‑150 (MQ=255)
                                                                                              tATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCGGTGGGAGCGGgc                                                                                                             >  2:2389369/1‑96 (MQ=255)
                                                                                               aTTGATCTTTGCCGGGGGGGTGTGCTGCGCGGCAAGGCCGGCCCCCACGGGTGCCCCACACGTGGGGGGGTGGTGGGGGCCGTGGGGGTGGGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCAc                                                      <  1:1600843/150‑1 (MQ=255)
                                                                                                  gATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGtt                                                                                                                                                       >  2:1403253/1‑50 (MQ=255)
                                                                                                  gATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGtt                                                                                                                                                       <  1:1403253/50‑1 (MQ=255)
                                                                                                   aTCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACCGGGTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGCCCCGGTTAGAAGTCCACGCTa                                                  >  2:3220792/1‑150 (MQ=255)
                                                                                                    tCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAg                                                 >  1:2134895/1‑150 (MQ=255)
                                                                                                          cccGCGTGGTCTGGTTGGCGGGCATGCCCGCACCCACCGGTCCAGCACATGTCTGTGGGTTGTGCTGGACGTGTGGGCGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGaaa                                           <  2:1838249/149‑1 (MQ=255)
                                                                                                          acggcgTCGTCTGCTTCGCGGGCATGCCCGCTCCCCCAGGTTCAGCACATGTCTGAGGGTTGTGCTGGACCTGTGGGCGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGaaa                                           <  2:1040684/147‑1 (MQ=255)
                                                                                                          aCCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGaaa                                           >  2:3829997/1‑150 (MQ=255)
                                                                                                             gcgTCGTCTGCTTCGCGGGCAGGCCCGCTCCCCCCGGTTCAGCACAAGTCTGTGGGTTGTGCTGGACCTGTGGGAGGGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCt                                        <  2:1013814/150‑1 (MQ=255)
                                                                                                               gtcgtcTGCTTCGGGGGCTGGCCCGCCCCCACAGGTTCAGCACAAGTCGGTGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGa                                      <  2:2730993/150‑1 (MQ=39)
                                                                                                                  gtcGGCTTCGCGGGCTGGCCGGCGCCCACTGGTTCAGCACACGTCTGTGGGTGGTGGGGGCCCTGGGGGGGGGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAg                                   <  1:370967/150‑1 (MQ=255)
                                                                                                                   tcTGCTTCGCGGGCATGCCCGCTCCCCCAGGTACAGCACACGGTTTTGGTTTGTGCTGGACCTGTGGGAGGGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAgg                                  <  1:3524573/150‑1 (MQ=38)
                                                                                                                   tcTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCCCGAAGAGGCCGTATATGGCGCCGCGGCGGGACCCGGGTAGAAGTCCAAGCTAGCCGAAAGCTGGTAgg                                  >  1:487560/1‑150 (MQ=38)
                                                                                                                     tGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGGTGTGCTGGACCGGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTc                                >  2:713003/1‑150 (MQ=255)
                                                                                                                       cTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCcg                              <  2:4966896/150‑1 (MQ=255)
                                                                                                                        ttCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCcgc                             >  1:3649571/1‑150 (MQ=255)
                                                                                                                         tCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTGTCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCcgcg                            >  2:1242262/1‑150 (MQ=255)
                                                                                                                          cgcgGGCATGCCCGCCCCCACCGGTTCCGCACACGTCTGAGGGTGGTGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCgcgg                           <  2:5017004/150‑1 (MQ=39)
                                                                                                                           gcgGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGATGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGAGAGGTCCgcggc                          >  2:3136634/1‑150 (MQ=255)
                                                                                                                                 aTGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGccccc                    >  2:3085930/1‑150 (MQ=255)
                                                                                                                                    cccGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGCAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAAc                 >  1:200093/1‑150 (MQ=25)
                                                                                                                                       gCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAAcccc              >  1:1562297/1‑150 (MQ=25)
                                                                                                                                         tCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGaa                                                            <  1:729677/102‑1 (MQ=14)
                                                                                                                                         tCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGaa                                                            >  2:729677/1‑102 (MQ=14)
                                                                                                                                         tCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAg            >  1:2030480/1‑150 (MQ=25)
                                                                                                                                           ccACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAggg          >  1:3921645/1‑150 (MQ=37)
                                                                                                                                             acaGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTc        >  2:1134865/1‑150 (MQ=255)
                                                                                                                                                ggTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTgg     >  1:1298763/1‑150 (MQ=255)
                                                                                                                                                  ttCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGt   >  1:1329578/1‑150 (MQ=255)
                                                                                                                                                  ttCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGt   <  2:588841/150‑1 (MQ=255)
                                                                                                                                                   tCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTc  >  2:3017001/1‑150 (MQ=255)
                                                                                                                                                   tCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGCCCACGCTAGCCGAAAGCTGATa                                    <  2:4544649/116‑1 (MQ=255)
                                                                                                                                                   tCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGCCCACGCTAGCCGAAAGCTGATa                                    >  1:4544649/1‑116 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTC  >  NC_002947/2960608‑2960904

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: