Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,319,585 C→G 5.3% D98H (GAT→CAT)  rcsC ← sensory histidine kinase RcsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,319,5850CG 5.3% 248.5 / 4.2 75D98H (GAT→CAT) rcsCsensory histidine kinase RcsC
Reads supporting (aligned to +/- strand):  ref base C (31/39);  new base G (2/2);  total (34/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 = 6.06e-01

AAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGAC  >  NC_000913/2319437‑2319731
                                                                                                                                                    |                                                                                                                                                  
aaaaaTACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATcc                                                                                                                                                   >  2:897591/1‑150 (MQ=255)
aaaaaTACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATcc                                                                                                                                                   >  2:466153/1‑150 (MQ=255)
   aaTACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGcc                                                                                                                                                <  2:1310640/150‑1 (MQ=255)
     tACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTg                                                                                                                                              <  1:557201/150‑1 (MQ=255)
     tACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGGGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTg                                                                                                                                              <  2:637677/150‑1 (MQ=255)
       cccGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCg                                                                                                                                            >  2:1053148/1‑150 (MQ=255)
        ccGGTTTAGATCGTAAGCCGCAGAAAAATTATCGAGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACGTCGCCAGGTTTTACTCATTGCGGAACAATCGGAGTCGGCAAAGAGCGGTTCAAACGCAGGCACATACGCCTGGGt                                                                                                                                           >  2:663764/1‑150 (MQ=255)
          ggTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCTAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGttt                                                                                                                                         <  2:1486826/150‑1 (MQ=255)
          ggTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCGGATGTGCCTGCGttt                                                                                                                                         <  1:1257219/150‑1 (MQ=255)
            tttAGATCGTAAGCCTCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTtct                                                                                                                                       <  1:1061955/150‑1 (MQ=255)
            tttAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTtct                                                                                                                                       <  2:1517237/150‑1 (MQ=255)
                  tCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCa                                                                                                                                 <  1:897591/150‑1 (MQ=255)
                    gTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCAcg                                                                                                                               <  1:273498/150‑1 (MQ=255)
                       aGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCgg                                                                                                                            >  2:274502/1‑150 (MQ=255)
                       aGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCgg                                                                                                                            >  1:1170236/1‑150 (MQ=255)
                         ccGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCgggg                                                                                                                          <  1:726138/150‑1 (MQ=255)
                          cGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCCCGCGGGGa                                                                                                                         >  2:1203018/1‑150 (MQ=255)
                              gAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGc                                                                                                                     <  2:254516/150‑1 (MQ=255)
                              gAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGc                                                                                                                     >  1:881654/1‑150 (MQ=255)
                              gAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGc                                                                                                                     >  2:874482/1‑150 (MQ=255)
                               aaaaaTTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGca                                                                                                                    <  1:1451836/150‑1 (MQ=255)
                               aaaaaTTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGca                                                                                                                    <  2:1403863/150‑1 (MQ=255)
                               aaaaaTTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGca                                                                                                                    <  2:603795/150‑1 (MQ=255)
                                  aaTTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCacac                                                                                                                 <  2:1170236/150‑1 (MQ=255)
                                       tCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGttt                                                                                                            <  2:266700/150‑1 (MQ=255)
                                         gcgcCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTc                                                                                                          >  1:1221663/1‑150 (MQ=255)
                                           gcCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTg                                                                                                        <  2:613737/150‑1 (MQ=255)
                                             cAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGcc                                                                                                      <  1:466153/150‑1 (MQ=255)
                                               gTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCCGGCACATCCGCCTGCGGTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGa                                                                                                    >  2:239088/1‑150 (MQ=255)
                                                tAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCTCCAGGTGTTACTCATTGCGGAACACTCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGAt                                                                                                   <  1:663764/150‑1 (MQ=255)
                                                 aGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATa                                                                                                  >  1:241470/1‑150 (MQ=255)
                                                   ctcATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCAGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAg                                                                                                <  2:1094601/148‑1 (MQ=255)
                                                    gcATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGc                                                                                               <  2:388123/150‑1 (MQ=255)
                                                    gcATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGc                                                                                               >  1:1034628/1‑150 (MQ=255)
                                                     cATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGc                                                                                                       >  2:30709/1‑141 (MQ=255)
                                                     cATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGc                                                                                                       <  1:30709/141‑1 (MQ=255)
                                                     cATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCg                                                                                              <  1:175140/150‑1 (MQ=255)
                                                         aaCCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCTAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCGCACCGTTTTCTGCCGATAAGCGAtgt                                                                                          >  2:614754/1‑148 (MQ=255)
                                                         aaCCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGAttt                                                                                          >  2:508952/1‑150 (MQ=255)
                                                              cGCCAATGACTCCAGAGAACCTCGCCAGGTGTCACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGc                                                                                     <  1:443815/150‑1 (MQ=255)
                                                                  aaTGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCTGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATg                                                                                 <  2:897520/150‑1 (MQ=255)
                                                                      aCTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACt                                                                             <  1:227210/150‑1 (MQ=255)
                                                                      aCTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACt                                                                             >  2:232999/1‑150 (MQ=255)
                                                                       cTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGAAAATCGCTTATCGGcagaaaac                                                                              <  2:1143641/148‑8 (MQ=255)
                                                                       cTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGAAAATCGCTTATCGGcagaaaac                                                                              >  1:1143641/1‑141 (MQ=255)
                                                                            gagaACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGct                                                                       <  1:547720/150‑1 (MQ=255)
                                                                             agaACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGctc                                                                      <  2:1131629/150‑1 (MQ=255)
                                                                                 ccTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTtc                                                                  <  1:1136498/150‑1 (MQ=255)
                                                                                  cTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTtca                                                                 <  1:117431/150‑1 (MQ=255)
                                                                                    cGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTtcatc                                                               >  2:885614/1‑150 (MQ=255)
                                                                                                  cATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCg                                                 >  1:463581/1‑150 (MQ=255)
                                                                                                  cATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCg                                                 >  2:150589/1‑150 (MQ=255)
                                                                                                  cATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCg                                                 >  2:786939/1‑150 (MQ=255)
                                                                                                  cATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCAACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCg                                                 >  2:41813/1‑150 (MQ=255)
                                                                                                       cGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGaa                                            >  1:1068399/1‑150 (MQ=255)
                                                                                                        ggAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAAt                                           <  1:1084528/150‑1 (MQ=255)
                                                                                                              aTCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAgcg                                     >  1:72758/1‑150 (MQ=255)
                                                                                                                cGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCt                                   >  2:1123802/1‑150 (MQ=255)
                                                                                                                 ggAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTg                                  <  1:786939/150‑1 (MQ=255)
                                                                                                                     tCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGcc                              <  2:1352605/150‑1 (MQ=255)
                                                                                                                      cGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCCCATCAGAGGTGCGTTGAATAAAGCGCTGCGCCt                             >  2:1019367/1‑150 (MQ=255)
                                                                                                                       ggCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTg                            >  1:195814/1‑150 (MQ=255)
                                                                                                                         cAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGAt                          >  2:640505/1‑150 (MQ=255)
                                                                                                                         aaaaCATCGGTTCAAACGCAGGAACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGATGTGCGTTGAATAAAGCGCTGAGCATGAt                          <  2:106096/149‑1 (MQ=255)
                                                                                                                           aaCAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCg                        >  1:58742/1‑150 (MQ=255)
                                                                                                                                gacgagaaaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACt                   >  1:1168635/7‑150 (MQ=255)
                                                                                                                                gacgagaaaCGCAGGCGGATGTGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACt                   >  1:876741/7‑150 (MQ=255)
                                                                                                                                 ggTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTc                  <  1:885614/150‑1 (MQ=255)
                                                                                                                                     cAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACCCCGTTTTCTGCCGATAAGCGATTTTCGGCGAGGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGCACTCAGAt              >  2:825886/1‑150 (MQ=255)
                                                                                                                                     cAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGAt              <  2:942697/150‑1 (MQ=255)
                                                                                                                                           cAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGa                                             >  2:268846/1‑113 (MQ=255)
                                                                                                                                           cAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGa                                             <  1:268846/113‑1 (MQ=255)
                                                                                                                                             ggCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTcc      >  2:1221155/1‑150 (MQ=255)
                                                                                                                                              gCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGTTTTTCGGCGATGTACTTCATGTTTTTTATTAAATCAGAGGTGCATTGAATTAAGTTCTGAGCCTGATCGGAATTCAGATTATATTctt     <  2:1221663/150‑3 (MQ=255)
                                                                                                                                               cacaTCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTg    >  1:1046472/1‑150 (MQ=255)
                                                                                                                                                 caTCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGAc  <  2:1503234/150‑1 (MQ=255)
                                                                                                                                                 caTCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGAc  <  1:1123802/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTGAC  >  NC_000913/2319437‑2319731

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: