Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 624,733 T→A 5.8% intergenic (‑223/‑152) fepB ← / → entC ferric enterobactin ABC transporter periplasmic binding protein/isochorismate synthase EntC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913624,7330TA 5.8% 267.1 / 3.4 82intergenic (‑223/‑152)fepB/entCferric enterobactin ABC transporter periplasmic binding protein/isochorismate synthase EntC
Reads supporting (aligned to +/- strand):  ref base T (38/39);  new base A (3/2);  total (41/41)
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

CCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGAT  >  NC_000913/624585‑624881
                                                                                                                                                    |                                                                                                                                                    
ccccTTTGGGGAGAGGGTTAGAGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTtctc                                                                                                                                                     >  2:713115/1‑150 (MQ=255)
 cccTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCa                                                                                                                                                    <  1:234319/150‑1 (MQ=255)
       gtggATAGGGTTAGGGTGTTGGGAACAGCCTGGACTGGTGCGAACATTAACCCTCACGACAGCCCTCACGCTGGAATGGAGAGGGTGCAGAACTGCGCAGGACATCACATTGCGGTTATGCGAATCCATCAATAATGCTTCTCATTTTCa                                                                                                                                              <  2:1184642/148‑1 (MQ=255)
       ggggAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCa                                                                                                                                              >  1:1241416/1‑150 (MQ=255)
       ggggAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCa                                                                                                                                              >  2:1188662/1‑150 (MQ=255)
         ggAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAtt                                                                                                                                            >  2:238226/1‑150 (MQ=255)
         ggAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAtt                                                                                                                                            <  2:1195841/150‑1 (MQ=255)
          gagaGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTg                                                                                                                                           <  1:1418097/150‑1 (MQ=255)
            gagGGTTAGAGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTa                                                                                                                                         <  1:713115/150‑1 (MQ=255)
               ggTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCa                                                                                                                                              >  1:126163/1‑142 (MQ=255)
               ggTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCa                                                                                                                                              <  2:126163/142‑1 (MQ=255)
                gTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACca                                                                                                                                     <  1:1188662/150‑1 (MQ=255)
                    gggTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAAc                                                                                                                                 >  1:1047765/1‑150 (MQ=255)
                           gggAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGc                                                                                                                          <  2:1241416/150‑1 (MQ=255)
                                 agccagcACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAAcccc                                                                                                                    >  2:94674/1‑150 (MQ=255)
                                       cACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACTCGGGGTTg                                                                                                              <  2:1213429/150‑1 (MQ=255)
                                       cACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTg                                                                                                              >  2:90271/1‑150 (MQ=255)
                                        aCTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGc                                                                                                             <  1:473511/150‑1 (MQ=255)
                                         cTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACTCGGGGTTGCa                                                                                                            <  1:24299/150‑1 (MQ=255)
                                         cTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCa                                                                                                            >  2:1033172/1‑150 (MQ=255)
                                            gTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGAt                                                                                                         >  2:899662/1‑150 (MQ=255)
                                            gTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGAt                                                                                                         <  1:1023680/150‑1 (MQ=255)
                                              gCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTg                                                                                                       >  2:154963/1‑150 (MQ=255)
                                                gAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCg                                                                                                     <  1:217778/150‑1 (MQ=255)
                                                   cATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAAcccc                                                                                                                    <  2:1062645/132‑1 (MQ=255)
                                                   cATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAAcccc                                                                                                                    >  1:1062645/1‑132 (MQ=255)
                                                        aCCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCa                                                                                             >  1:301455/1‑150 (MQ=255)
                                                         cccTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCaa                                                                                            >  1:235667/1‑150 (MQ=255)
                                                          ccTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAag                                                                                           <  1:39264/150‑1 (MQ=255)
                                                            tCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAagag                                                                                         >  2:1404869/1‑150 (MQ=255)
                                                               ccccAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGAGATCACATTGCGCTTATGCGAATCCATAAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTg                                                                                      <  2:1322772/150‑1 (MQ=255)
                                                               ccccAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTg                                                                                      <  1:792275/150‑1 (MQ=255)
                                                                   cccccTCTCCCTTCCAGGGTGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCACATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACAt                                                                                  >  2:1388505/3‑150 (MQ=255)
                                                                   cccccTCTCCCTTCCAGGGTGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCACATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACAt                                                                                  >  2:1383591/3‑150 (MQ=255)
                                                                    ccccTCTCCCTTCCAGGGTGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATa                                                                                 >  1:680395/2‑150 (MQ=255)
                                                                         cACCCTGGAAGGGAGAGGGGGCAGAACGGCGCACGACATAACATTGCGCTTATGCGAATCCATCAATAATCCATCACATTTTCATTCTAACCACAACCAGATGCAACCCCGAGATCCAGATTGCGTTACCTCAAGAGTTGACATAGTgcg                                                                            >  1:614075/1‑150 (MQ=255)
                                                                           cccTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGcgcg                                                                          >  2:626711/1‑150 (MQ=255)
                                                                           cccTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGcgcg                                                                          <  1:566744/150‑1 (MQ=255)
                                                                           cccTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGcgcg                                                                          <  1:571405/150‑1 (MQ=255)
                                                                           cccTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGcgcg                                                                          >  2:881611/1‑150 (MQ=255)
                                                                            ccTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTc                                                                                              <  1:1245705/129‑1 (MQ=255)
                                                                            ccTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTc                                                                                              >  2:1245705/1‑129 (MQ=255)
                                                                                 aaGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCt                                                                    >  1:644682/1‑150 (MQ=255)
                                                                                    ggAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCtttt                                                                 >  2:814997/1‑150 (MQ=255)
                                                                                     gtgagGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTa                                                                <  2:310840/148‑1 (MQ=255)
                                                                                      tgagGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCACATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAg                                                               <  1:1388505/149‑1 (MQ=255)
                                                                                      tgagGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCACATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAg                                                               <  1:1383591/149‑1 (MQ=255)
                                                                                      agagGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAg                                                               >  2:228721/1‑150 (MQ=255)
                                                                                         gggggCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGtt                                                            >  2:839130/1‑150 (MQ=255)
                                                                                            ggCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGc                                                         <  1:39800/150‑1 (MQ=255)
                                                                                            ggCAGAACGGCGCAGGACAGCCCATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTACCACAACCAGATGCCACCCCGAGTTGCAGATTGCGTTACCTCAAGCGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCg                                                        >  2:20768/1‑150 (MQ=255)
                                                                                              cAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGa                                                       >  2:1060517/1‑150 (MQ=255)
                                                                                              cAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGa                                                       >  2:719347/1‑150 (MQ=255)
                                                                                                gAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGAcc                                                     >  1:927141/1‑150 (MQ=255)
                                                                                                  aCGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGa                                                   >  1:224180/1‑150 (MQ=255)
                                                                                                      cgcAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAat                                               <  1:238226/150‑1 (MQ=255)
                                                                                                        cAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAatat                                             >  1:364069/1‑150 (MQ=255)
                                                                                                         aGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGttt                                                                       <  2:884369/123‑1 (MQ=255)
                                                                                                         aGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGttt                                                                       >  1:884369/1‑123 (MQ=255)
                                                                                                          ggACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATaa                                           <  2:534283/150‑1 (MQ=255)
                                                                                                           gACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATaaa                                          <  1:719347/150‑1 (MQ=255)
                                                                                                           gACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATaaa                                          <  1:1023800/150‑1 (MQ=255)
                                                                                                               tCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGAt                                      <  1:90271/150‑1 (MQ=255)
                                                                                                                caaaTTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATa                                     >  2:232997/4‑150 (MQ=255)
                                                                                                                 acaTTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATaa                                    >  1:1330806/1‑150 (MQ=255)
                                                                                                                       cgcTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCatta                              >  1:924928/1‑150 (MQ=255)
                                                                                                                       cgcTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCatta                              >  2:166223/1‑150 (MQ=255)
                                                                                                                        gcTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCattat                             <  1:881611/150‑1 (MQ=255)
                                                                                                                         cTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCattatt                            <  2:934659/150‑1 (MQ=255)
                                                                                                                          ttATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCAttatta                           <  2:1047765/150‑1 (MQ=255)
                                                                                                                           tATGCGAATCCATCAATAATGATTCTCATTTTAATAGTAACCACAACCAGATGCAACCCAGAGTTGCAGATTGCGTTACCTAAAGAGTTGACAAAGTGCGCGTTTGCTTTTAGGTTAGAGACCGAAAATATAAATGATAATCATTATTaa                          >  1:544308/1‑150 (MQ=255)
                                                                                                                               cGAATCCATCAATAATGCTTCTCATTTTCATTGTAACTACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGcc                      >  1:172684/1‑150 (MQ=255)
                                                                                                                               cGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGca                      >  1:1431740/1‑149 (MQ=255)
                                                                                                                                gAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGGTGACATAGTGCGCGTTTGCTTTTAGGGTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCt                     <  2:655952/150‑1 (MQ=255)
                                                                                                                                    ccATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTAt                 <  2:1330806/150‑1 (MQ=255)
                                                                                                                                     cATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATc                <  1:814997/150‑1 (MQ=255)
                                                                                                                                     cATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATc                <  1:1095087/150‑1 (MQ=255)
                                                                                                                                      aTCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCa               >  1:1007605/1‑150 (MQ=255)
                                                                                                                                      aTCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCa               <  2:330004/150‑1 (MQ=255)
                                                                                                                                         aataatGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCAttt            >  2:518142/1‑150 (MQ=255)
                                                                                                                                          ataatGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCAtttt           <  1:839130/150‑1 (MQ=255)
                                                                                                                                               gCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTgga      <  1:228721/150‑1 (MQ=255)
                                                                                                                                                 ttCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTggagg    <  2:845581/150‑1 (MQ=255)
                                                                                                                                                   ctcATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGgat  >  2:250790/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGAT  >  NC_000913/624585‑624881

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