Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,318,922 C→T 5.2% A319T (GCA→ACA)  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,318,9220CT 5.2% 246.4 / 2.7 77A319T (GCA→ACA) rcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base C (35/37);  new base T (2/2);  total (38/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 = 7.81e-01

GCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCA  >  NC_000913/2318778‑2319069
                                                                                                                                                |                                                                                                                                                   
gcgcGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGtt                                                                                                                                                  >  7:167732/1‑148 (MQ=255)
gcgcGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGtt                                                                                                                                                  <  1:340370/148‑1 (MQ=255)
gcgcGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCCCTTTCCGCCGGAATGATAATACGTCGCTCTTACATCAGTGCGAGAGTAAACAATGCATCTCCGGCAAGCACAATCAGCAAAATTGAGtt                                                                                                                                                  >  2:257618/1‑148 (MQ=255)
  gcgGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAg                                                                                                                                               <  2:58240/149‑1 (MQ=255)
   cgGAGGCGACAATCTTGCGATTGAACTGCGCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAAGTGCGTTAAg                                                                                                                                               >  6:394976/1‑148 (MQ=255)
      aGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATc                                                                                                                                           <  3:90655/149‑1 (MQ=255)
       ggCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCTTTAAGGATc                                                                                                                                           >  7:349935/1‑148 (MQ=255)
          gACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACa                                                                                                                                       >  5:302773/1‑149 (MQ=255)
             aaTCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGAACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGc                                                                                                                                    >  1:65403/1‑149 (MQ=255)
             aaTCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGAACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGc                                                                                                                                    >  1:149029/1‑149 (MQ=255)
                   gCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAAt                                                                                                                                <  8:167732/147‑1 (MQ=255)
                   gCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAAt                                                                                                                                <  6:302773/147‑1 (MQ=255)
                   gCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTAGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAAt                                                                                                                                <  1:244818/147‑1 (MQ=255)
                   gCGATTGAACTGCTCATGTTCTTCCAGTCGAAGGTCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGaa                                                                                                                                 <  6:247252/146‑1 (MQ=255)
                     gATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAg                                                                                                                                               <  5:121510/130‑1 (MQ=255)
                     gATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAg                                                                                                                                               >  6:121510/1‑130 (MQ=255)
                             tGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCAc                                                                                                                    >  5:263680/1‑149 (MQ=255)
                                tCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCt                                                                                                                  >  5:380200/1‑148 (MQ=255)
                                  attttcttcCAGTAGCAGGGCGTCTATTTCCGCCGGAATGTAAATACGTCGCTAGTACATCCGTGCGAGAGTAAACAATGCAGCTCCTGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTAt                                                                                                               <  5:335896/146‑1 (MQ=255)
                                  aTGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTa                                                                                                                <  5:44683/148‑1 (MQ=255)
                                  aTGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTAt                                                                                                               >  7:82439/1‑149 (MQ=255)
                                    gTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGc                                                                                                                                                     >  5:161816/1‑109 (MQ=255)
                                    gTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGc                                                                                                                                                     <  6:161816/109‑1 (MQ=255)
                                        ttcCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAg                                                                                                                                               >  2:411104/1‑111 (MQ=255)
                                          ccAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATg                                                                                                                               >  5:315940/1‑125 (MQ=255)
                                          ccAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATg                                                                                                                               <  6:315940/125‑1 (MQ=255)
                                          ccAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAAcc                                                                                                          <  6:263680/146‑1 (MQ=255)
                                          ccAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCgg                                                                                                        <  5:394976/148‑1 (MQ=255)
                                          ccAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGc                                                                                                       >  5:397243/1‑149 (MQ=255)
                                           cAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGc                                                                                                       <  8:82439/148‑1 (MQ=255)
                                            aGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCAc                                                                                                     <  7:369461/149‑1 (MQ=255)
                                                gCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGTCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGaa                                                                                                 <  8:349935/149‑1 (MQ=255)
                                                        tCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGAt                                                                                         >  8:121882/1‑149 (MQ=255)
                                                          gCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGa                                                                                          <  4:109767/146‑1 (MQ=255)
                                                          gCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGa                                                                                          >  3:109767/1‑146 (MQ=255)
                                                          gCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGc                                                                                       >  2:286584/1‑149 (MQ=255)
                                                                  ccgGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTACGGAAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGAt                                                                               <  4:9915/149‑1 (MQ=255)
                                                                     gAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCTATGAg                                                                            >  6:395818/1‑149 (MQ=255)
                                                                         gAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTg                                                                        >  4:59824/1‑149 (MQ=255)
                                                                          aaaaTACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGACAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCg                                                                                 <  2:193521/139‑1 (MQ=255)
                                                                          aaaaTACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGACAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCg                                                                                 >  1:193521/1‑139 (MQ=255)
                                                                            aaTACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCa                                                                     <  6:397243/149‑1 (MQ=255)
                                                                                    gCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGTGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTttc                                                             >  8:309224/1‑149 (MQ=255)
                                                                                    gCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGTGTTAACGATCAACATGCCAATGCGTTCCAGCACCTAATAAACCGTCACCGAATACACGATGCTTAGCGATGAGGCTGGCAGAtattt                                                              >  8:303536/1‑148 (MQ=255)
                                                                                    gCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGAttttt                                                              <  4:37819/148‑1 (MQ=255)
                                                                                    gCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTttc                                                             <  3:59824/149‑1 (MQ=255)
                                                                                     cTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTttct                                                            >  2:133478/1‑149 (MQ=255)
                                                                                        gTACATCCGTGCGAGAGTAAACAATGCATCTCCGGCAAGCACATTCAGCAAAATTGCGTTAATGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTttcttc                                                          <  7:364676/148‑1 (MQ=255)
                                                                                        gTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGc                                                                                                                                    >  8:304490/1‑74 (MQ=255)
                                                                                        gTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGc                                                                                                                                    <  7:304490/74‑1 (MQ=255)
                                                                                        gTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTttctt                                                           <  6:49214/147‑1 (MQ=255)
                                                                                         tccATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGACTGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGCGGGTGGCAGATTTTTCTTCa                                                         <  6:329733/146‑1 (MQ=255)
                                                                                         tACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAATTTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAg                                                        <  1:257618/149‑1 (MQ=255)
                                                                                         tACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGGGGCAGATTTTTCTTCAg                                                        >  5:329733/1‑149 (MQ=255)
                                                                                               cGTGCGATAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCa                                                   <  2:65403/148‑1 (MQ=255)
                                                                                               cGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCa                                                   <  2:149029/148‑1 (MQ=255)
                                                                                                      gagTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCg                                            >  5:30998/1‑148 (MQ=255)
                                                                                                        gTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGAt                                                                               >  5:136131/1‑111 (MQ=255)
                                                                                                        gTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGAt                                                                               <  6:136131/111‑1 (MQ=255)
                                                                                                         tAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATAAACCGGCACCGAATACACGATACTTAGCGATGAGGGTTGCAGATTTTTCTTCAGCACCCGATACCGGaa                                         >  7:360219/1‑148 (MQ=255)
                                                                                                           aaCAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATc                                                                                                              <  4:65915/77‑1 (MQ=255)
                                                                                                           aaCAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATc                                                                                                              >  3:65915/1‑77 (MQ=255)
                                                                                                            aCAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCt                                     <  7:121882/149‑1 (MQ=255)
                                                                                                            aCAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCt                                     >  1:175275/1‑149 (MQ=255)
                                                                                                              aaTGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTc                                   >  8:66144/1‑149 (MQ=255)
                                                                                                              aaTGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATCCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTc                                   >  3:74658/1‑149 (MQ=255)
                                                                                                               aTGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTcc                                  <  1:286584/149‑1 (MQ=255)
                                                                                                                 gCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCa                                                         >  6:155410/1‑124 (MQ=255)
                                                                                                                 gCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATGTTTCTTCa                                                         <  5:155410/124‑1 (MQ=255)
                                                                                                                             aGCACATTCAGCAAAATTGCGTAAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCa                    <  6:380200/149‑1 (MQ=255)
                                                                                                                              gCACATTCAGCAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTGATCAACCGGCACCGAATACACGATGGTTAGCGATCAGGGTGGCAGATTTTGCTTAAGCACCAGCTCCCGGAACCCGTCCGTAAAGCCAAACCAg                   >  6:115823/1‑148 (MQ=255)
                                                                                                                                acaTTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGa                 <  6:30998/149‑1 (MQ=255)
                                                                                                                                 caTTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAg                >  7:16616/1‑149 (MQ=255)
                                                                                                                                  aTTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGt                                >  5:83007/1‑132 (MQ=255)
                                                                                                                                  aTTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGt                                <  6:83007/132‑1 (MQ=255)
                                                                                                                                  aTTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGc               >  3:21693/1‑149 (MQ=255)
                                                                                                                                      aGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTc           >  7:403495/1‑149 (MQ=255)
                                                                                                                                       gCAAAATTGTGTTAATGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTcc          <  7:303536/149‑1 (MQ=255)
                                                                                                                                       gCAAAATTGTGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTcc          <  7:309224/149‑1 (MQ=255)
                                                                                                                                                cGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCa  <  7:66144/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCGCGGAGGCGACAATCTTGCGATTGAACTGCTCATGTTCTTCCAGTCGCAGGGCGTCGCTTTCCGCCGGAATGAAAATACGTCGCTCGTACATCCGTGCGAGAGTAAACAATGCAGCTCCGGCAAGCACATTCAGCAAAATTGCGTTAAGGATCAACATGCGAATGCGTTCCAGCACCTTATCAACCGGCACCGAATACACGATGCTTAGCGATGAGGGTGGCAGATTTTTCTTCAGCACCAGCTCCCGGAACCCTTCCGTATAGCCAAACCAGGAGCGTTCCTGCATCCA  >  NC_000913/2318778‑2319069

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