Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,960,781 G→A 16.1% intergenic (+82/+60) 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,7810GA16.1% 274.2 / 39.1 118intergenic (+82/+60)PP_2589/PP_2590aldehyde dehydrogenase family protein/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (67/31);  new base A (8/10);  total (75/43)
Fisher's exact test for biased strand distribution p-value = 6.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-01

Only 100 of 122 total aligned reads displayed.

TCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTAGTCCTTGTC  >  NC_002947/2960633‑2960930
                                                                                                                                                    |                                                                                                                                                     
tCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGa                                                                                                                                                      <  1:4039656/150‑1 (MQ=255)
tCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGa                                                                                                                                                      <  2:1349481/150‑1 (MQ=255)
 cTGTCACTGCATTCGTTCGACAAGTATTCGCTACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAAc                                                                                                                                                     <  2:900937/150‑1 (MQ=37)
 cTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTTTTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAAc                                                                                                                                                     <  2:901187/150‑1 (MQ=37)
 cTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGAc                                                                                                                                                     <  2:1573844/150‑1 (MQ=255)
 cTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAAc                                                                                                                                                     <  2:175092/150‑1 (MQ=37)
   gTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTTGCGGGCGCGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCt                                                                                                                                                   <  1:2056586/150‑1 (MQ=255)
   gTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCt                                                                                                                                                   >  2:2942741/1‑150 (MQ=255)
   gTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCt                                                                                                                                                   <  2:4686678/150‑1 (MQ=37)
     cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                                                 <  1:4068400/150‑1 (MQ=255)
     cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgt                                                                                                                                                 <  1:1792848/150‑1 (MQ=255)
     cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCtgt                                                                                                                                                 <  1:2591277/150‑1 (MQ=37)
     cACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTCCAGCACAACCCTCAGACTTGTGCTGAACCtgt                                                                                                                                                 <  1:2591202/150‑1 (MQ=37)
      aCTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCtgtg                                                                                                                                                >  2:4795552/1‑150 (MQ=255)
       cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                                               >  2:3926096/1‑150 (MQ=255)
       cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                                               >  2:3925974/1‑150 (MQ=255)
       cTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTgg                                                                                                                                               <  2:508164/150‑1 (MQ=255)
             tcgttcgACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGGTCAGCACAAGTCTGAGGGTTGTGCTGGACCGGTGGGAGCgg                                                                                                                                         >  1:2371872/1‑150 (MQ=255)
                   gACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGgcgcgc                                                                                                                                   >  1:1359039/1‑150 (MQ=255)
                   gACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGAGGGAGCGGgcgcgc                                                                                                                                   >  1:4025785/1‑150 (MQ=255)
                    aCAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGcc                                                                                                                                  >  1:4522198/1‑150 (MQ=255)
                    aCAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGcc                                                                                                                                  >  1:4522143/1‑150 (MQ=255)
                    aCAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGcc                                                                                                                                  >  1:4522118/1‑150 (MQ=255)
                    aCAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGcc                                                                                                                                  >  1:4522088/1‑150 (MQ=255)
                    aCAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGcc                                                                                                                                  >  2:4502236/1‑150 (MQ=255)
                     cAAGTATTCGCAACTCAAGCCCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGccc                                                                                                                                 <  2:4096843/150‑1 (MQ=255)
                      aaGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCcg                                                                                                                                >  1:607643/1‑150 (MQ=255)
                          tttCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTTATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGTTTCAGCACAAGTCTGTGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGaa                                                                                                                            <  1:5122271/149‑1 (MQ=255)
                          aTTCCCAACCCACGACCACCTGGTACCATGGGCGCGGCTGATTTATTGATCATTGCCCGCGTCGTCTGCTTCGCGGGCAGGCCCGCTCCCCCTGGTTCCGCACAAGACTGAGGGTTGTGCTGCACCTGTGGGAGCGGGCGCGCCCGCGaa                                                                                                                            <  1:2465139/150‑1 (MQ=37)
                           ttCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAag                                                                                                                           >  1:1658512/1‑150 (MQ=255)
                           ttCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAag                                                                                                                           >  2:2602825/1‑150 (MQ=255)
                              gCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAgg                                                                                                                        >  2:791999/1‑150 (MQ=255)
                              gCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCACGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACGTGTGGGAGCGGGCGCGCCCGCGAAGTgg                                                                                                                        >  2:1947244/1‑150 (MQ=255)
                               cAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGc                                                                                                                       <  2:1106484/150‑1 (MQ=255)
                                aaCTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGcc                                                                                                                      >  2:1297895/1‑150 (MQ=255)
                                   tCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCGGAGGGGTGTGCTGGCCCAGTGGGAGCGGGCGCGCCCGCGACGAGGCCGTa                                                                                                                   >  1:4733555/1‑150 (MQ=39)
                                          caccTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGCGGCCGTATCTGgcg                                                                                                            >  1:2227140/1‑150 (MQ=255)
                                                  aCCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGCAGAGGCCGTATCTGGCGCCgcggcg                                                                                                    >  1:3210829/1‑150 (MQ=255)
                                                   ccAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGCAGAGGCCGTATCTGGCGCCCcggcgg                                                                                                   >  1:2621092/1‑150 (MQ=255)
                                                      gTTGCGCGGCTTATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCCGCACAAGTCGGAGGGTTGTGCTGGACCGGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCgcggcgccag                                                                                                <  2:1498884/150‑5 (MQ=39)
                                                      gTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGAc                                                                                                <  2:3378205/150‑1 (MQ=255)
                                                          cggggCTGTTGTATTGATCTTTGACCGCGTCGTCTGCGCGGCGGGCTGGCCCGGACCCACAGGTTCCGCACATTTATGTGGGTTGGGCTGGACCTGTGGTAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCgg                                                                                            <  2:4733555/147‑1 (MQ=255)
                                                             gggtgttgttttgtttATTGACCGCGTCTTCTGCTTGGGGGGCTTGCCCGCTCCCACTGGTTCTGCACTTGTCTGTGGGTTGTGCTGGTCGTGTGGGTGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGTCCCGGTTa                                                                                         <  2:2979968/144‑1 (MQ=255)
                                                             ggCTTATGTTTTGATCTTTGCCGGCGTCGTCTGCTTCGCGGGCAGGCCCGCTCCCCCCGGTTCAGCACACGTCTGTGGGTTGTGGTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTa                                                                                         <  2:1909988/150‑1 (MQ=38)
                                                             ggCTGATTTTTTTATCATTGACCGGGTCGTCTGCTTGGCGGGCATGCCCGCCCCCACCGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTa                                                                                         <  2:1910022/150‑1 (MQ=39)
                                                             ggCTGATGATTTGATCATTGACCGCGTCGTCTGCTTCGGGGGCATGCCCGCTCCCACAGGTTCAGCAAATGTGTTAGGTTTGTGCTGGACCTTTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTa                                                                                         <  2:4511396/150‑1 (MQ=38)
                                                              tcttttttATTTATCTTTGACCGGGTCGTCTGCTTCGCGGGCTTGCCCGCTCCCACTGTTTCAGCACAAGTCTGTGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAg                                                                                        <  2:3919143/143‑1 (MQ=37)
                                                                    gTATTGATAATTGACCGCGTCGTGTGCTTCGCGGGCATGCCCGCTCCCCCCGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTcc                                                                                  <  2:3177951/150‑1 (MQ=255)
                                                                     tATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCa                                                                                 >  1:2595955/1‑150 (MQ=255)
                                                                      aTTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGCCCCGGGTTGGAGTCCAc                                                                                >  1:1380862/1‑150 (MQ=255)
                                                                       ttGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGGCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGGTGGACCTGTGGGAGCGGGCGCGGCCGCGTAGTGGCCGTATCTGGGGCCGGGGCGGGACCCGGTTAGAAGTACACg                                                                               >  1:1899359/1‑150 (MQ=255)
                                                                        tGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCCCTGGTTCCGCACAAGTCTGAGGGTTGGGCTGGACCTGTGGGCGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGc                                                                              <  2:2597127/150‑1 (MQ=255)
                                                                           tCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAg                                                                           >  2:2477294/1‑150 (MQ=255)
                                                                                  ccGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACCGGTTCAGCACACGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAg                                                                    <  2:2120342/150‑1 (MQ=255)
                                                                                   cgcgTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCCGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATATGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGc                                                                   >  1:1111512/1‑150 (MQ=255)
                                                                                   cgcgTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCGGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCCGTTAGCAGTCCACGCTAGCCGCAAGc                                                                   >  1:2256481/1‑150 (MQ=255)
                                                                                      gtcgtcTGCTTCGCGGGCTTGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGa                                                                <  2:246820/150‑1 (MQ=255)
                                                                                      gtcgtcTGCTTCGCGGGCATGCCCGCGCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCGGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGa                                                                <  2:4456679/150‑1 (MQ=255)
                                                                                      ggcgtcGGCTTCGCGGGCTGGCCCGCGCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGa                                                                <  1:2942741/148‑1 (MQ=255)
                                                                                          tcTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAgg                                                            >  2:1537814/1‑150 (MQ=255)
                                                                                           cTGCTTCGCGGGCATGCCCGCTGCCACAGGTTCAGCACAAGTCTGAGGGTGGGGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGt                                                           >  2:23896/1‑150 (MQ=255)
                                                                                           cTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGt                                                           >  2:2577954/1‑150 (MQ=255)
                                                                                                        aTGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGccccc                                              >  1:3689988/1‑150 (MQ=255)
                                                                                                             cGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAAccc                                         >  2:1712706/1‑150 (MQ=255)
                                                                                                                 cccACGGTTTCCGCACACGTGTGTGGGTTGTGCTGCCCGTGGGGGCGGGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAgg                                     <  1:3925974/150‑1 (MQ=255)
                                                                                                                  ccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAggg                                    >  1:2315693/1‑150 (MQ=255)
                                                                                                                  ccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAggg                                    >  2:1050637/1‑150 (MQ=255)
                                                                                                                  ccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAggg                                    >  2:4195280/1‑150 (MQ=255)
                                                                                                                  ccACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCGGATAGGTCCGCGGCGCCCCCAACCCCAggg                                    >  1:2315611/1‑150 (MQ=255)
                                                                                                                   cacaGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGAGAGGTCCGCGGCGCCCCCAACCCCAGGGt                                   >  2:2919204/1‑150 (MQ=25)
                                                                                                                    acaGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCCCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTACGCGGCGCCCCCAACCCCAGGGTc                                  >  1:1257107/1‑150 (MQ=37)
                                                                                                                      aGGTCCAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGCGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATATGTCCGCGGCGCCCCCAACCCCAGGGTCTg                                >  2:1789322/1‑150 (MQ=25)
                                                                                                                        gtccAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTggg                              >  2:4200592/4‑150 (MQ=37)
                                                                                                                        gTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTggg                              >  2:4048411/1‑150 (MQ=255)
                                                                                                                        gTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGGTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTggg                              >  1:59568/1‑150 (MQ=255)
                                                                                                                         tccAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGTCGCCCCCAACCCCAGGGTCTGGGt                             >  1:4182234/3‑150 (MQ=37)
                                                                                                                         tccAGCACAACCCTCAGACTTGTGCTGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGt                             >  1:3457255/3‑150 (MQ=25)
                                                                                                                          tCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTc                            >  1:1755478/1‑150 (MQ=255)
                                                                                                                           cAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCa                           >  1:240160/1‑150 (MQ=255)
                                                                                                                            aGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAg                          >  1:4704305/1‑150 (MQ=255)
                                                                                                                             gCACAAGTCTGAGGGTGGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGt                         <  2:4024535/150‑1 (MQ=255)
                                                                                                                              cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTg                        >  2:3686274/1‑150 (MQ=255)
                                                                                                                              cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTg                        >  2:3686232/1‑150 (MQ=255)
                                                                                                                              cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTg                        >  2:3686446/1‑150 (MQ=255)
                                                                                                                              cacaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTg                        >  2:1121917/1‑150 (MQ=255)
                                                                                                                               acaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGtggccaa                       >  2:4888030/1‑144 (MQ=255)
                                                                                                                               acaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTgg                       >  2:2145992/1‑150 (MQ=255)
                                                                                                                               acaAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTgg                       >  2:1542248/1‑150 (MQ=255)
                                                                                                                                    tCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCacac                  >  2:3234140/1‑150 (MQ=255)
                                                                                                                                      tgtggggTGGGCGGGCCCTGGGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACAccc                <  1:4795552/147‑1 (MQ=255)
                                                                                                                                      tGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGCCCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCACGGTCTGGGTCAGTGGCACAccc                >  1:4709173/1‑150 (MQ=255)
                                                                                                                                      tGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACAccc                >  1:1428047/1‑150 (MQ=255)
                                                                                                                                               tgCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTAGTcc       >  2:5065183/1‑150 (MQ=255)
                                                                                                                                                gCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTAGTCCt      >  1:3863208/1‑150 (MQ=255)
                                                                                                                                                 cTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTAGTCCtt     <  2:2371872/150‑1 (MQ=255)
                                                                                                                                                  tGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTAGTCCTTg    <  2:1359039/150‑1 (MQ=255)
                                                                                                                                                   ggACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAATTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTAGTCCTTGt   >  1:1340049/1‑150 (MQ=255)
                                                                                                                                                    gACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTAGTCCTTGTc  >  1:1771093/1‑150 (MQ=255)
                                                                                                                                                    gACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTAGTCCTTGTc  >  2:754653/1‑150 (MQ=255)
                                                                                                                                                    gACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTAGTCCTTGTc  >  2:4088093/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCTGTCACTGCATTCGTTCGACAAGTATTCGCAACTCAAGACCACCTGGTACCAGTTGCGCGGCTGATGTATTGATCATTGACCGCGTCGTCTGCTTCGCGGGCATGCCCGCTCCCACAGGTTCAGCACAAGTCTGAGGGTTGTGCTGGACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGTATCTGGCGCCGCGGCGGGACCCGGTTAGAAGTCCACGCTAGCCGAAAGCTGATAGGTCCGCGGCGCCCCCAACCCCAGGGTCTGGGTCAGTGGCACACCCCAGTAGTCCTTGTC  >  NC_002947/2960633‑2960930

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: