Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,314,509 T→C V341A (GTC→GCC)  rcsD → phosphotransfer intermediate protein in two‑component regulatory system with RcsBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,314,5090TC97.7% 288.1 / ‑3.0 87V341A (GTC→GCC) rcsDphosphotransfer intermediate protein in two‑component regulatory system with RcsBC
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (36/49);  total (37/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

ATTCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGTCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACGATC  >  NC_000913/2314412‑2314604
                                                                                                  |                                                                                               
aTTCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGccc                                                                                               <  2:958627/100‑1 (MQ=255)
aTTCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGccc                                                                                               <  2:1436682/100‑1 (MQ=255)
 ttCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGcccc                                                                                              >  1:1674337/1‑100 (MQ=255)
  tCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCcccc                                                                                             >  2:2116459/1‑100 (MQ=255)
  tCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCcccc                                                                                             <  1:2060076/100‑1 (MQ=255)
  tCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCcccc                                                                                             >  1:947010/1‑100 (MQ=255)
  tCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCcccc                                                                                             <  1:1576145/100‑1 (MQ=255)
  tCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCcccc                                                                                             <  2:661959/100‑1 (MQ=255)
   ctgctgCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCa                                                                                            <  2:1608991/100‑1 (MQ=255)
     gctgcCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGc                                                                                          <  2:1703386/100‑1 (MQ=255)
     gctgcCACTGCTGCTGAACATCGGTTTGCTGGCGCGGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGc                                                                                          <  2:2560726/100‑1 (MQ=255)
         ccACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCAcgg                                                                                      >  2:2002599/1‑100 (MQ=255)
          cACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCAcggc                                                                                     <  1:1236739/100‑1 (MQ=255)
           aCTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCGCCAGCCGCAGTACAGAAAACGCCCCCAGCAcggcg                                                                                    <  2:1222876/100‑1 (MQ=255)
            ctgctgctgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTCTACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCAcggcgg                                                                                   <  2:1842108/100‑1 (MQ=255)
            ctgctgctgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCAcggcgg                                                                                   <  2:548337/100‑1 (MQ=255)
            ctgctgctgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCAcggcgg                                                                                   <  1:1371310/100‑1 (MQ=255)
            ctgctgctgAACATCGGTTTGCTGGCGCTGGCGTTATTAGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCAcggcgg                                                                                   <  2:1189586/100‑1 (MQ=255)
             tgctgctgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGt                                                                                  <  2:931490/100‑1 (MQ=255)
              gctgctgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTc                                                                                 <  2:1289793/100‑1 (MQ=255)
                tgctgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCaa                                                                               >  1:1745257/1‑100 (MQ=255)
                  ctgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATa                                                                             <  2:1917332/100‑1 (MQ=255)
                   tgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATaa                                                                            >  1:1801697/1‑100 (MQ=255)
                    gAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAAc                                                                           >  1:671331/1‑100 (MQ=255)
                    gAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAAc                                                                           <  1:519149/100‑1 (MQ=255)
                    gAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAAc                                                                           >  2:1533502/1‑100 (MQ=255)
                      aCATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGa                                                                         >  1:2270555/1‑100 (MQ=255)
                         tCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAAtt                                                                      >  2:305470/1‑100 (MQ=255)
                          cGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTg                                                                     >  2:255476/1‑100 (MQ=255)
                          cGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTg                                                                     >  1:805230/1‑100 (MQ=255)
                               tGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCAt                                                                >  1:363122/1‑100 (MQ=255)
                                 ctggcgctggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCAttt                                                              <  2:2014172/100‑1 (MQ=255)
                                 ctggcgctggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCAttt                                                              <  1:582803/100‑1 (MQ=255)
                                 ctggcgctggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCAttt                                                              >  1:2207784/1‑100 (MQ=255)
                                   ggcgctggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTa                                                            >  2:920811/1‑100 (MQ=255)
                                    gcgctggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTAc                                                           <  1:1003041/100‑1 (MQ=255)
                                       ctggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACggg                                                        <  2:1772309/100‑1 (MQ=255)
                                       ctggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATAGCGCATTTTACggg                                                        <  1:1543715/100‑1 (MQ=255)
                                        tggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACTGGc                                                       <  1:2029416/100‑1 (MQ=255)
                                        tggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGc                                                       <  1:2303513/100‑1 (MQ=255)
                                         ggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGca                                                      >  1:2642653/1‑100 (MQ=255)
                                           cgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGcaa                                                     >  2:1870136/1‑100 (MQ=255)
                                           cgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGcaat                                                    >  1:1406929/1‑100 (MQ=255)
                                            gTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGcaatc                                                   <  1:1546102/100‑1 (MQ=255)
                                            gTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGcaatc                                                   >  2:398631/1‑100 (MQ=255)
                                             ttATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGcaatca                                                  >  1:1805108/1‑100 (MQ=255)
                                              tATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGcaatcaa                                                 <  2:1402660/100‑1 (MQ=255)
                                                    gCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAaga                                           <  1:1423090/100‑1 (MQ=255)
                                                    gCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAaga                                           <  1:1133208/100‑1 (MQ=255)
                                                    gCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAaga                                           <  1:1798033/100‑1 (MQ=255)
                                                     cTATACCACATGCCGCCATTTCTCCAGCCGGAGTACAGAAAACGCTCCCAGCACGGCGGGAAACAACGAATTGCGCATTTTACGGGCAATCAAGGGagag                                          <  2:2441077/100‑1 (MQ=255)
                                                         aCCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAg                                      <  1:511879/100‑1 (MQ=255)
                                                          ccACATGCCGCCATTTCTACAGCCGCAGTACAGAAAACGCCCCCAGTACGGCGCACAATAACGAATTGCACCTTGTACGGGCAATCAATGAAGAGATCGt                                     >  2:890748/1‑100 (MQ=255)
                                                           cacaTTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGtc                                    <  1:433706/100‑1 (MQ=255)
                                                            acaTTCCGCCAGTTCTCCAGCCGCAGTACAGAAAACGCCCGCAGCACGGCGGTCAATAACGAATTGCGCATTCTACGGGCAATCAAAGAAGAGATAGtct                                   <  1:2327404/100‑1 (MQ=255)
                                                               ttCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCAc                                <  1:2436600/100‑1 (MQ=255)
                                                               ctCCGCCATTTCGCCAGCCGCAGTACAGACAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCAc                                <  2:1528257/99‑1 (MQ=255)
                                                                tCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCAct                               <  1:2069824/100‑1 (MQ=255)
                                                                tCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCAct                               >  1:694340/1‑100 (MQ=255)
                                                                  cGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCActgc                             <  1:1031845/100‑1 (MQ=255)
                                                                  cGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCActgc                             >  1:1265057/1‑100 (MQ=255)
                                                                  cGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCActgc                             >  1:1575231/1‑100 (MQ=255)
                                                                   gCCATTTCTCCAGCCGCAGTACAGAAAACGTCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCActgct                            <  1:1287479/100‑1 (MQ=255)
                                                                   gCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCAAGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGCGATAGTCTCActgct                            <  2:516470/100‑1 (MQ=255)
                                                                    ccATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCActgctg                           >  1:25306/1‑100 (MQ=255)
                                                                     cATTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACtgctgc                          <  2:126573/100‑1 (MQ=255)
                                                                      aTTTCTCCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAGTTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGcc                         >  2:2322301/1‑100 (MQ=255)
                                                                          ctcCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTc                     >  1:1080237/1‑100 (MQ=255)
                                                                           tcCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCg                    <  1:311225/100‑1 (MQ=255)
                                                                           tcCAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCg                    >  1:85906/1‑100 (MQ=255)
                                                                            ccAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGGGATAGTCTCACTGCTGCCGCTCgg                   <  1:1346282/100‑1 (MQ=255)
                                                                            ccAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCgg                   <  2:1069090/100‑1 (MQ=255)
                                                                            ccAGCCGCAGTACAGAAAACGCCCCCAGCACGGCAGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCgg                   <  2:2018143/100‑1 (MQ=255)
                                                                             cAGCCGCAGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGc                  >  1:1594790/1‑100 (MQ=255)
                                                                                  gccgTACAGAAAACGCCCCCCGCACGGCTGTCCATAACGAAGTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCctgct             <  2:1157883/97‑1 (MQ=255)
                                                                                    aGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTgg           >  2:2418686/1‑100 (MQ=255)
                                                                                    aGTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTgg           <  1:329647/100‑1 (MQ=255)
                                                                                     gTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGt          >  2:501318/1‑100 (MQ=255)
                                                                                     gTACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGt          <  2:796305/100‑1 (MQ=255)
                                                                                      tACAGAAAACGCCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGt          <  1:1895192/100‑1 (MQ=255)
                                                                                      tACAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGtt         <  2:941308/100‑1 (MQ=255)
                                                                                       aCAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTc        <  2:1872247/100‑1 (MQ=255)
                                                                                        cAGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCa       >  2:2302408/1‑100 (MQ=255)
                                                                                         aGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCAc      >  2:807744/1‑100 (MQ=255)
                                                                                         aGAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCAc      >  1:22914/1‑100 (MQ=255)
                                                                                          gAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACg     >  1:2137752/1‑100 (MQ=255)
                                                                                          gAAAACGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACg     <  2:291152/100‑1 (MQ=255)
                                                                                            aaaCGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACGAt   >  2:805179/1‑100 (MQ=255)
                                                                                            aaaCGCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACGAt   >  1:2262845/1‑100 (MQ=255)
                                                                                             aaCGTCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACGATc  >  2:1782745/1‑100 (MQ=255)
                                                                                                  |                                                                                               
ATTCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGTCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACGATC  >  NC_000913/2314412‑2314604

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