Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,319,581 A→G 5.4% V99A (GTG→GCG) ‡ rcsC ← sensory histidine kinase RcsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,319,5810AG 5.4% 221.6 / 2.7 74V99A (GTG→GCG) ‡rcsCsensory histidine kinase RcsC
Reads supporting (aligned to +/- strand):  ref base A (33/37);  new base G (2/2);  total (35/39)
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.40e-01

CGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTG  >  NC_000913/2319432‑2319729
                                                                                                                                                     |                                                                                                                                                    
cGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGca                                                                                                                                                      <  2:910609/150‑1 (MQ=255)
     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                                                                                                                  <  2:1403863/150‑1 (MQ=255)
                                    aaaaaTTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGca                                                                                                                  <  1:1451836/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                                                                                             >  1:1034628/1‑150 (MQ=255)
                                                         gcATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGc                                                                                             <  2:388123/150‑1 (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                                               >  2:150589/1‑150 (MQ=255)
                                                                                                       cATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCg                                               >  1:463581/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)
                                                                                                                                                     |                                                                                                                                                    
CGATTAAAAATACCCGGTTTAGATCGTAAGCCGCAGAAAAATTATCGCGCCAGTAGCGCATAAACCACGCCAATGACTCCAGAGAACCTCGCCAGGTGTTACTCATTGCGGAACAATCGGAGTCGGCAAACAGCGGTTCAAACGCAGGCACATCCGCCTGCGTTTCTCGTCCACGCGGGGAAAGCACACCGTTTTCTGCCGATAAGCGATTTTCGGCGATGTACTTCAGCTCTTTCATCACATCAGAGGTGCGTTGAATAAAGCGCTGAGCCTGATCGGAACTCAGATTAAATTCCTG  >  NC_000913/2319432‑2319729

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: