Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,958,833 T→C 9.7% 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,958,8330TC 9.7% 318.1 / 17.5 113intergenic (+56/+86)PP_2589/PP_2590aldehyde dehydrogenase family protein/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (52/50);  new base C (3/8);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 116 total aligned reads displayed.

GGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGG  >  NC_002947/2958684‑2958980
                                                                                                                                                     |                                                                                                                                                   
gggCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGtt                                                                                                                                                     <  1:2537852/150‑1 (MQ=255)
 ggCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGtt                                                                                                                                                     >  2:264904/1‑149 (MQ=255)
 ggCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGtt                                                                                                                                                     <  1:264904/149‑1 (MQ=255)
 ggCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTc                                                                                                                                                    >  2:1946098/1‑150 (MQ=255)
  gCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCa                                                                                                                                                   <  1:1946098/150‑1 (MQ=255)
   aaaGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAg                                                                                                                                                  <  2:2896661/149‑1 (MQ=255)
        aGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACaa                                                                                                                                             <  2:1485411/150‑1 (MQ=255)
            ggggTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCt                                                                                                                                         >  1:2495342/1‑150 (MQ=255)
            ggggTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCt                                                                                                                                         >  2:1971581/1‑150 (MQ=255)
              ggTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGa                                                                                                                                       <  2:2276853/150‑1 (MQ=255)
                   tgccactgacccagaccctggggttgggggcgccgcggacctatcagctttcggctagcgtggacttctaaccgggtcccgccgcggcgccagatacggcctcttcgcgggcgcGCCCGCTCCCACAGGTCCAGCACAAGTCTGAGGGtt                                                                                                                                  >  1:1015898‑M1/115‑150 (MQ=255)
                             tGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                               >  1:2857987/1‑143 (MQ=255)
                             tGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                               <  2:2857987/143‑1 (MQ=255)
                             tGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGAcc                                                                                                                        <  1:1971581/150‑1 (MQ=255)
                              gTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCt                                                                                                                       <  2:2859236/150‑1 (MQ=255)
                               tCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtg                                                                                                                      >  2:2346861/1‑150 (MQ=255)
                                cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCtgt                                                                                                                     <  1:2537923/150‑1 (MQ=37)
                                 aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCtgtg                                                                                                                    <  1:1429449/150‑1 (MQ=37)
                                  cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGa                                                                                                                          >  1:1767613/1‑143 (MQ=255)
                                  cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGa                                                                                                                          <  2:1767613/143‑1 (MQ=255)
                                  cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                   <  2:2495342/150‑1 (MQ=255)
                                   tGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                     >  1:2231953/1‑147 (MQ=255)
                                   tGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                     <  2:2231953/147‑1 (MQ=255)
                                   tGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                   >  1:1881966/1‑149 (MQ=255)
                                   tGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                   <  2:1881966/149‑1 (MQ=255)
                                        tcgttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCgg                                                                                                             <  1:2346861/150‑1 (MQ=255)
                                          gttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGcac                                                                                                                                               <  2:2444326/114‑1 (MQ=255)
                                          gttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGcac                                                                                                                                               >  1:2444326/1‑114 (MQ=255)
                                              gACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCGAAAAGCTGAGGGTTGTGCTGGACCTGTGAGATCggaagagc                                                                                                       >  1:2336506/1‑144 (MQ=255)
                                              gACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGc                                                                                                                              >  1:1962099/1‑127 (MQ=255)
                                              gACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGc                                                                                                                              <  2:1962099/127‑1 (MQ=255)
                                               aCAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGgcatgcc                                                                                                      <  1:1067678/150‑6 (MQ=18)
                                                cAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGgcatgccc                                                                                                     <  1:574355/150‑7 (MQ=17)
                                                cAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGgcatgccc                                                                                                     >  1:1811230/1‑144 (MQ=17)
                                                 aaGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCcg                                                                                                    >  2:2696592/1‑150 (MQ=255)
                                                     aTTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGaa                                                                                                <  2:1811230/150‑1 (MQ=17)
                                                      ttCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCGCGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAag                                                                                               <  1:2683002/150‑1 (MQ=17)
                                                         gCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                   <  1:794010/127‑1 (MQ=255)
                                                         gCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                   >  2:794010/1‑127 (MQ=255)
                                                         gCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAgg                                                                                            >  2:823404/1‑150 (MQ=255)
                                                            aCTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGa                                                                                                                                       <  2:644266/104‑1 (MQ=255)
                                                              tCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACTTGTGGGAGCGGGCGCGCCCGCGGAGAGGCCGTa                                                                                       >  2:2691476/1‑150 (MQ=255)
                                                                aaGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGtt                                                                                                                                  <  1:39651/105‑1 (MQ=255)
                                                                 aGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAGGTCTGAGGGTtgtg                                                                                                                               >  1:2040455/1‑107 (MQ=255)
                                                                 aGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTtgtg                                                                                                                               <  2:2040455/107‑1 (MQ=255)
                                                                 aGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCt                                                                                    >  1:208119/1‑150 (MQ=255)
                                                                  gACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCAGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGccc                                                                                                     >  2:2560953/1‑132 (MQ=255)
                                                                   accaccTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                   >  2:1211185/1‑117 (MQ=255)
                                                                   accaccTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                   <  1:1211185/117‑1 (MQ=255)
                                                                           gTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcgg                                                                          >  2:2761001/1‑150 (MQ=255)
                                                                           gTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCAATGGTTTCAGGGTAGGTCAGTTCCGGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcgg                                                                          >  1:361382/1‑150 (MQ=37)
                                                                             cgtggacttctaaccgggtcccgccgcggcgccagatacggcctcttcgcgggcgcGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGa                                                                                                                          <  1:2197053‑M1/44‑1 (MQ=255)
                                                                                    ttctaaccgggtcccgccgcggcgccagatacggcctcttcgcgggcgcGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGAc                                                                                                                         >  2:1563477‑M1/50‑94 (MQ=255)
                                                                                      gcgGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACGGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGt                                                               <  1:338863/150‑1 (MQ=255)
                                                                                       taaccgggtcccgccgcggcgccagatacggcctcttcgcgggcgcGCCCGCTCCCACAGGTTCAGCACAAGt                                                                                                                                           >  1:1359401‑M1/47‑73 (MQ=255)
                                                                                       taaccgggtcccgccgcggcgccagatacggcctcttcgcgggcgcGCCCGCTCCCACAGGTTCAGCACAAGt                                                                                                                                           <  2:1359401‑M1/27‑1 (MQ=255)
                                                                                        ggCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTa                                                             >  1:881711/1‑150 (MQ=255)
                                                                                        ggCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTa                                                             >  1:2878589/1‑150 (MQ=255)
                                                                                        ggCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGAGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTa                                                             >  2:404826/1‑150 (MQ=255)
                                                                                            gATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGc                                                                                                               >  2:2336152/1‑96 (MQ=255)
                                                                                            gATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGc                                                                                                               <  1:2336152/96‑1 (MQ=255)
                                                                                             aTGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAgg                                                                                                                                     >  2:2218811/1‑73 (MQ=255)
                                                                                             aTGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAgg                                                                                                                                     <  1:2218811/73‑1 (MQ=255)
                                                                                                tATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGgcgc                                                                               <  2:257115/124‑1 (MQ=255)
                                                                                                tATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGgcgc                                                                               >  1:257115/1‑124 (MQ=255)
                                                                                                tATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCa                                                     >  1:1395135/1‑150 (MQ=255)
                                                                                                 aTTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCAc                                                    <  2:208119/150‑1 (MQ=255)
                                                                                                    gATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCt                                                 >  2:380724/1‑150 (MQ=255)
                                                                                                       cactGACCGCGTCGTCTGCTTCGCGGGCACGCCCGCTCCCACAGGTTCAGGACAAGTTTCAGGGTTGTGCTGGACCTGCGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGc                                              <  1:2691476/147‑1 (MQ=38)
                                                                                                        aTTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGTTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGcc                                             >  2:2721811/1‑150 (MQ=255)
                                                                                                              cgcgTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                     <  1:2041341/72‑1 (MQ=255)
                                                                                                              cgcgTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                     >  2:2041341/1‑72 (MQ=255)
                                                                                                              cgcgTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGt                                                                                        <  1:2239662/101‑1 (MQ=255)
                                                                                                              cgcgTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGt                                                                                        >  2:2239662/1‑101 (MQ=255)
                                                                                                                 gtcgtcTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGa                                    >  1:396268/1‑150 (MQ=255)
                                                                                                                    ggcctcttcgcggtctcGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGc                                                                                                               >  2:1190069‑M1/18‑72 (MQ=255)
                                                                                                                    ggcctcttcgcgggcgcGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGc                                                                                                               <  1:1190069‑M1/55‑1 (MQ=255)
                                                                                                                       tGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTc                              >  2:1372669/1‑150 (MQ=255)
                                                                                                                       tGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTACACGCTAGCCGAAAGCTGATAGGTc                              >  2:2149349/1‑150 (MQ=255)
                                                                                                                       ctcttcgcgggcgcGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                     >  1:1493770‑M1/15‑63 (MQ=255)
                                                                                                                          ttCGCGGGCATGCCCGCGCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCcgc                           <  1:1203791/150‑1 (MQ=255)
                                                                                                                           tCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTTGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCcgcg                          >  2:1154378/1‑150 (MQ=255)
                                                                                                                                  cATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGcccc                   <  1:823404/150‑1 (MQ=255)
                                                                                                                                    tGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCa                 <  2:1395135/150‑1 (MQ=255)
                                                                                                                                       CCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGcc                                                                                          >  1:1209820‑M1/1‑74 (MQ=255)
                                                                                                                                       CCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCAc                                                    <  2:1827976‑M1/112‑1 (MQ=255)
                                                                                                                                       CCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCAc                                                    >  1:1827976‑M1/1‑112 (MQ=255)
                                                                                                                                        CGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAAccc             <  1:2420475‑M1/150‑1 (MQ=255)
                                                                                                                                            CCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTACACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg         >  1:2149349‑M1/1‑150 (MQ=255)
                                                                                                                                             CCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGaa                                                          >  1:2451864‑M1/1‑100 (MQ=255)
                                                                                                                                             CCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAggg        <  2:942185‑M1/150‑1 (MQ=255)
                                                                                                                                              CACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGt       <  1:668935‑M1/150‑1 (MQ=255)
                                                                                                                                               ACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGtt                                                              >  2:1116143‑M1/1‑94 (MQ=255)
                                                                                                                                               ACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGtt                                                              <  1:1116143‑M1/94‑1 (MQ=255)
                                                                                                                                               ACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTc      <  1:1353683‑M1/150‑1 (MQ=255)
                                                                                                                                                CAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCt     <  2:965137‑M1/150‑1 (MQ=255)
                                                                                                                                                 AGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTg    <  1:1439208‑M1/150‑1 (MQ=255)
                                                                                                                                                  GGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTgg   >  1:2696592‑M1/1‑150 (MQ=255)
                                                                                                                                                   GTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTggg  <  2:1333599‑M1/150‑1 (MQ=255)
                                                                                                                                                   GTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTggg  <  1:1882542‑M1/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGGCAAGCAGTCGGGGTTTGGCCGCGATCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGG  >  NC_002947/2958684‑2958980

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: