Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,319,581 A→G 8.8% V99A (GTG→GCG) ‡ rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,319,5810AG 8.8% 244.7 / 13.0 91V99A (GTG→GCG) ‡rcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base A (36/47);  new base G (4/4);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

CCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGA  >  NC_000913/2319444‑2319730
                                                                                                                                         |                                                                                                                                                     
cccGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAACCCACGCCAATGACTCAAGAGAACCTCGCCAGGTGTTAATCATTGCGGACCAATCGGAGTCGGCAAACAGCGGTTCAAACGCCGCCACATCCGCCTGCg                                                                                                                                           >  2:1074837/1‑150 (MQ=255)
 ccGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGt                                                                                                                                          >  1:385085/1‑150 (MQ=255)
  cgttttAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGtt                                                                                                                                         <  1:461983/147‑1 (MQ=255)
  cGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGtt                                                                                                                                         <  2:308120/150‑1 (MQ=255)
  cGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGtt                                                                                                                                         <  2:993625/150‑1 (MQ=255)
  cGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGtt                                                                                                                                         <  1:1181968/150‑1 (MQ=255)
  cGGATTAGATCGTAAGACGCAGAAAAATTATCGCGCCATTAGCGCATAAACCAAGCCAATGACTCCAGAGAACCTCGCAAGGTGTTACTCATTGCGGAACAATCGGAGTCGTCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGtt                                                                                                                                         <  1:93383/150‑1 (MQ=255)
   ggTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGttt                                                                                                                                        <  2:1114001/150‑1 (MQ=255)
   ggTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATACGCCTGCGttt                                                                                                                                        >  1:701352/1‑150 (MQ=255)
   ggTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTACAAACGCAGGCACATCCGCCTGCGttt                                                                                                                                        >  1:887134/1‑150 (MQ=255)
   ggTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGGCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGttt                                                                                                                                        <  2:1108263/150‑1 (MQ=255)
     tttAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTtct                                                                                                                                      >  2:1228289/1‑150 (MQ=255)
       tAGATCGTAAGCCGGAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCg                                                                                                                                    <  2:949308/150‑1 (MQ=255)
       tAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCg                                                                                                                                    <  2:1716622/150‑1 (MQ=255)
           tCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCTCCAGGTGTTACTCATTGCGGATCAATCGGAGTCTGCAAACAGCGGTTCAAACGCAGGCACCTCCGCCTGCGTTTCTCGTCCa                                                                                                                                <  1:884615/150‑1 (MQ=255)
               aaGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAcgcg                                                                                                                            >  1:676383/1‑150 (MQ=255)
                 gCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCggg                                                                                                                          <  1:1389738/150‑1 (MQ=255)
                 gCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCggg                                                                                                                          >  2:1354839/1‑150 (MQ=255)
                  ccGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCgggg                                                                                                                         <  1:402385/150‑1 (MQ=255)
                   cGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGGGTTACTAATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTCCGTTTCTCGTCCACGCGGGGa                                                                                                                        >  1:770648/1‑150 (MQ=255)
                       tcaaaaTTATCGCGCCAGTAGCGCATAAACCCCGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGGTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGc                                                                                                                    <  1:1078555/148‑1 (MQ=255)
                         aaaaTTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGcac                                                                                                                  >  1:23429/1‑150 (MQ=255)
                         aaaaTTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGcac                                                                                                                  >  2:1238503/1‑150 (MQ=255)
                           aaTTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCacac                                                                                                                >  1:757701/1‑150 (MQ=255)
                           aaTTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCacac                                                                                                                >  1:1670136/1‑150 (MQ=255)
                             ttATCGCGCCATTATCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTTCTCATTTCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACAcca                                                                                                              <  2:1024746/150‑2 (MQ=255)
                                 cgcgCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGGGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGtttt                                                                                                          <  1:974942/150‑1 (MQ=255)
                                  gcgcCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTc                                                                                                         <  2:614489/150‑1 (MQ=255)
                                  gcgcCAGTAGAGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTc                                                                                                         <  2:391183/150‑1 (MQ=255)
                                   cgcCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCt                                                                                                        <  2:1683267/150‑1 (MQ=255)
                                   cgcCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCt                                                                                                        <  1:636691/150‑1 (MQ=255)
                                     ccAGTAGCGCATAAACCACGCCAATGACTCCAGAGTACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAACACACCGTTTTCTGc                                                                                                      <  2:85972/150‑1 (MQ=255)
                                      cAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGcc                                                                                                     <  1:1354839/150‑1 (MQ=255)
                                           gcgcATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATaa                                                                                                >  2:522950/1‑150 (MQ=255)
                                            cgcATAACCCACGCCAATGCCTCCAGCGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCATACCGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAACGCACCACGTTTTCTGCCGATAAg                                                                                               >  2:78588/1‑150 (MQ=255)
                                             gcATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGc                                                                                              <  1:647966/150‑1 (MQ=255)
                                              cATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCg                                                                                             >  1:1680651/1‑150 (MQ=255)
                                                tAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGAt                                                                                           <  1:458489/150‑1 (MQ=255)
                                                       cGCCAATTACTCCATAGAACCTCGCCAGGTGTTACTCATTGCGTAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGGGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGc                                                                                    <  1:1074837/150‑1 (MQ=255)
                                                       cGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGc                                                                                    <  2:527727/150‑1 (MQ=255)
                                                         ccAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGAAAATCGCTTATCGGcag                                                                                  <  1:265952/150‑3 (MQ=255)
                                                         ccAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCAGCCGATAAGCGATTTTCGGCGa                                                                                  >  2:393154/1‑150 (MQ=255)
                                                           aaTGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATg                                                                                >  1:1094887/1‑150 (MQ=255)
                                                              gACTCCAGATAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGAAAATCGCTTATCGGcagaaaac                                                                             <  2:1506348/150‑8 (MQ=255)
                                                                 tCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTc                                                                          <  2:1670136/150‑1 (MQ=255)
                                                                 tCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTc                                                                          <  2:757701/150‑1 (MQ=255)
                                                                    agagAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGc                                                                       <  2:716416/150‑1 (MQ=255)
                                                                          ccTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCTTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCttta                                                                 >  1:458297/1‑149 (MQ=255)
                                                                           cTCGCCAGGTGTTACTCATTGCCCAACAATCGGCGTCGGCAACCCGCGGCCCACACGCAGGCACATCCGCCAGCGTTTCTCGTCCACGCGGCGAAAGCACACCGTTATCAGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTtca                                                                <  2:1680651/150‑1 (MQ=255)
                                                                              gCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTtcatca                                                             <  1:751249/150‑1 (MQ=255)
                                                                              gCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTtcatca                                                             <  2:701352/150‑1 (MQ=255)
                                                                              gCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGAATTTCGGCGATGTACTTCAGCTCTTtcatca                                                             >  1:1056183/1‑150 (MQ=255)
                                                                               ccAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTcatcac                                                            >  1:603794/1‑150 (MQ=255)
                                                                                  ggTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTcatcacatc                                                         >  1:45340/1‑150 (MQ=255)
                                                                                   gtgtTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTcatcacatca                                                        <  1:1190481/150‑1 (MQ=255)
                                                                                   gtgtTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTATCGTCCACGCGGGGAAAGCAAACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACATCAGCTCTTTcaacacatca                                                        >  2:1353916/1‑150 (MQ=255)
                                                                                     gtTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCaga                                                      <  2:547788/150‑1 (MQ=255)
                                                                                         cTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTg                                                  <  2:385085/150‑1 (MQ=255)
                                                                                             ttGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGtt                                              <  2:23429/150‑1 (MQ=255)
                                                                                              tGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTg                                             <  1:623144/150‑1 (MQ=255)
                                                                                                 ggAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAAt                                          >  2:1458639/1‑150 (MQ=255)
                                                                                                    aCAATCGGAGTCGGCAAACAGCGGTTCAAACGAAGGCACATCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAg                                      <  1:375192/150‑1 (MQ=255)
                                                                                                     cAATCGTAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGTGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAg                                      <  2:352044/150‑1 (MQ=255)
                                                                                                             gTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCCCTGAGc                              >  2:74092/1‑150 (MQ=255)
                                                                                                              tCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGcc                             <  1:1123814/150‑1 (MQ=255)
                                                                                                              gcGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGcc                             <  2:676383/149‑1 (MQ=255)
                                                                                                               cGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCt                            <  2:968862/150‑1 (MQ=255)
                                                                                                                ggCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCATGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTg                           >  1:554869/1‑150 (MQ=255)
                                                                                                                 gCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACGCCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCTCTGAGCCTGa                          >  2:1480156/1‑150 (MQ=255)
                                                                                                                   aaaCAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATc                        <  1:1168458/150‑1 (MQ=255)
                                                                                                                     aCAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCgg                      >  2:693428/1‑150 (MQ=255)
                                                                                                                      cAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACAGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCAAATAAGAGTTGCGTTGAATAAAGAGCTGAGCCTGATCGGa                     >  2:878053/1‑150 (MQ=255)
                                                                                                                           gTTCAAACGCAGGCACATCCGCCTGCGTTTGTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCa                <  2:1161813/150‑1 (MQ=255)
                                                                                                                           gTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCa                >  1:1524544/1‑150 (MQ=255)
                                                                                                                           gTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCa                >  2:1016208/1‑150 (MQ=255)
                                                                                                                           gTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCa                >  2:1455692/1‑150 (MQ=255)
                                                                                                                            gagaaaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGAGTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATAGGATCTCAg               >  2:1153775/4‑150 (MQ=255)
                                                                                                                             tCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGa              >  2:328783/1‑150 (MQ=255)
                                                                                                                             agaaaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGa              >  2:1194326/3‑150 (MQ=255)
                                                                                                                              ctaaCTCAGGCACATCCTCCTGCGTTTCTCGGCCACGCTGGGAACGCACACCGTTTTCAGCCGATTAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTGGAATAAAGCGCTGAGCCTGATCGGTACTCAGAt             <  1:850495/148‑1 (MQ=255)
                                                                                                                              cAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGAt             >  2:1118755/1‑150 (MQ=255)
                                                                                                                                aaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTa           >  2:216204/1‑150 (MQ=255)
                                                                                                                                aaCGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTa           <  1:1228289/150‑1 (MQ=255)
                                                                                                                                  cGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGAGCGGAACTCAGATTtaa         >  2:827694/1‑150 (MQ=255)
                                                                                                                                  cGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTaaa         >  1:618840/1‑150 (MQ=255)
                                                                                                                                  cGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTaaa         >  1:117807/1‑150 (MQ=255)
                                                                                                                                   gCAGGCACATCCGCCTGCGTTTCTCGTCCTCGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAAt        <  1:1156795/150‑1 (MQ=255)
                                                                                                                                   gCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAAt        >  1:480315/1‑150 (MQ=255)
                                                                                                                                     aGGCACATCCGCATGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTc      <  2:1261336/150‑1 (MQ=255)
                                                                                                                                      ggCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCTGATTAAATTcc     <  1:516518/150‑1 (MQ=255)
                                                                                                                                       gCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCt    <  2:887134/150‑1 (MQ=255)
                                                                                                                                         acaTCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGa  <  2:1524544/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                     
CCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGA  >  NC_000913/2319444‑2319730

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