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

ATTCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑TCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACGATC  >  NC_000913/2314412‑2314604
                                                                                                  |                                                                                               
aTTCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑ccc                                                                                               <  2:958627/100‑1 (MQ=255)
aTTCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑ccc                                                                                               <  2:1436682/100‑1 (MQ=255)
 ttCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑cccc                                                                                              >  1:1674337/1‑100 (MQ=255)
  tCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑Ccccc                                                                                             >  2:2116459/1‑100 (MQ=255)
  tCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑Ccccc                                                                                             <  1:2060076/100‑1 (MQ=255)
  tCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑Ccccc                                                                                             >  1:947010/1‑100 (MQ=255)
  tCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑Ccccc                                                                                             <  1:1576145/100‑1 (MQ=255)
  tCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑Ccccc                                                                                             <  2:661959/100‑1 (MQ=255)
   ctgctgCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCa                                                                                            <  2:1608991/100‑1 (MQ=255)
     gctgcCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGc                                                                                          <  2:1703386/100‑1 (MQ=255)
     gctgcCACTGCTGCTGAACATCGGTTTGCTGGCGCGGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGc                                                                                          <  2:2560726/100‑1 (MQ=255)
         ccACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCAcgg                                                                                      >  2:2002599/1‑100 (MQ=255)
          cACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCAcggc                                                                                     <  1:1236739/100‑1 (MQ=255)
           aCTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCGCCAGCCGCAGTACAGAAAACG‑CCCCCAGCAcggcg                                                                                    <  2:1222876/100‑1 (MQ=255)
            ctgctgctgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTCTACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCAcggcgg                                                                                   <  2:1842108/100‑1 (MQ=255)
            ctgctgctgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCAcggcgg                                                                                   <  2:548337/100‑1 (MQ=255)
            ctgctgctgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCAcggcgg                                                                                   <  1:1371310/100‑1 (MQ=255)
            ctgctgctgAACATCGGTTTGCTGGCGCTGGCGTTATTAGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCAcggcgg                                                                                   <  2:1189586/100‑1 (MQ=255)
             tgctgctgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGt                                                                                  <  2:931490/100‑1 (MQ=255)
              gctgctgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTc                                                                                 <  2:1289793/100‑1 (MQ=255)
                tgctgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCaa                                                                               >  1:1745257/1‑100 (MQ=255)
                  ctgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATa                                                                             <  2:1917332/100‑1 (MQ=255)
                   tgAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATaa                                                                            >  1:1801697/1‑100 (MQ=255)
                    gAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAAc                                                                           >  1:671331/1‑100 (MQ=255)
                    gAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAAc                                                                           <  1:519149/100‑1 (MQ=255)
                    gAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAAc                                                                           >  2:1533502/1‑100 (MQ=255)
                      aCATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGa                                                                         >  1:2270555/1‑100 (MQ=255)
                         tCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAAtt                                                                      >  2:305470/1‑100 (MQ=255)
                          cGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTg                                                                     >  2:255476/1‑100 (MQ=255)
                          cGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTg                                                                     >  1:805230/1‑100 (MQ=255)
                               tGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCAt                                                                >  1:363122/1‑100 (MQ=255)
                                 ctggcgctggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCAttt                                                              <  2:2014172/100‑1 (MQ=255)
                                 ctggcgctggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCAttt                                                              <  1:582803/100‑1 (MQ=255)
                                 ctggcgctggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCAttt                                                              >  1:2207784/1‑100 (MQ=255)
                                   ggcgctggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTa                                                            >  2:920811/1‑100 (MQ=255)
                                    gcgctggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTAc                                                           <  1:1003041/100‑1 (MQ=255)
                                       ctggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACggg                                                        <  2:1772309/100‑1 (MQ=255)
                                       ctggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATAGCGCATTTTACggg                                                        <  1:1543715/100‑1 (MQ=255)
                                        tggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACTGGc                                                       <  1:2029416/100‑1 (MQ=255)
                                        tggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGc                                                       <  1:2303513/100‑1 (MQ=255)
                                         ggcgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGca                                                      >  1:2642653/1‑100 (MQ=255)
                                           cgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACGCCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGcaa                                                     >  2:1870136/1‑100 (MQ=255)
                                           cgTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGcaat                                                    >  1:1406929/1‑100 (MQ=255)
                                            gTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGcaatc                                                   <  1:1546102/100‑1 (MQ=255)
                                            gTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGcaatc                                                   >  2:398631/1‑100 (MQ=255)
                                             ttATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGcaatca                                                  >  1:1805108/1‑100 (MQ=255)
                                              tATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGcaatcaa                                                 <  2:1402660/100‑1 (MQ=255)
                                                    gCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAaga                                           <  1:1423090/100‑1 (MQ=255)
                                                    gCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAaga                                           <  1:1133208/100‑1 (MQ=255)
                                                    gCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAaga                                           <  1:1798033/100‑1 (MQ=255)
                                                     cTATACCACATGCCGCCATTTCTCCAGCCGGAGTACAGAAAACG‑CTCCCAGCACGGCGGGAAACAACGAATTGCGCATTTTACGGGCAATCAAGGGagag                                          <  2:2441077/100‑1 (MQ=255)
                                                         aCCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAg                                      <  1:511879/100‑1 (MQ=255)
                                                          ccACATGCCGCCATTTCTACAGCCGCAGTACAGAAAACG‑CCCCCAGTACGGCGCACAATAACGAATTGCACCTTGTACGGGCAATCAATGAAGAGATCGt                                     >  2:890748/1‑100 (MQ=255)
                                                           cacaTTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGtc                                    <  1:433706/100‑1 (MQ=255)
                                                            acaTTCCGCCAGTTCTCCAGCCGCAGTACAGAAAACG‑CCCGCAGCACGGCGGTCAATAACGAATTGCGCATTCTACGGGCAATCAAAGAAGAGATAGtct                                   <  1:2327404/100‑1 (MQ=255)
                                                               ttCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCAc                                <  1:2436600/100‑1 (MQ=255)
                                                               ctCCGCCATTTCGCCAGCCGCAGTACAGACAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCAc                                <  2:1528257/99‑1 (MQ=255)
                                                                tCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCAct                               <  1:2069824/100‑1 (MQ=255)
                                                                tCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCAct                               >  1:694340/1‑100 (MQ=255)
                                                                  cGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCActgc                             <  1:1031845/100‑1 (MQ=255)
                                                                  cGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCActgc                             >  1:1265057/1‑100 (MQ=255)
                                                                  cGCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCActgc                             >  1:1575231/1‑100 (MQ=255)
                                                                   gCCATTTCTCCAGCCGCAGTACAGAAAACG‑TCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCActgct                            <  1:1287479/100‑1 (MQ=255)
                                                                   gCCATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCAAGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGCGATAGTCTCActgct                            <  2:516470/100‑1 (MQ=255)
                                                                    ccATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCActgctg                           >  1:25306/1‑100 (MQ=255)
                                                                     cATTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACtgctgc                          <  2:126573/100‑1 (MQ=255)
                                                                      aTTTCTCCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAGTTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGcc                         >  2:2322301/1‑100 (MQ=255)
                                                                          ctcCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTc                     >  1:1080237/1‑100 (MQ=255)
                                                                           tcCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCg                    <  1:311225/100‑1 (MQ=255)
                                                                           tcCAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCg                    >  1:85906/1‑100 (MQ=255)
                                                                            ccAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGGGATAGTCTCACTGCTGCCGCTCgg                   <  1:1346282/100‑1 (MQ=255)
                                                                            ccAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCgg                   <  2:1069090/100‑1 (MQ=255)
                                                                            ccAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCAGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCgg                   <  2:2018143/100‑1 (MQ=255)
                                                                             cAGCCGCAGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGc                  >  1:1594790/1‑100 (MQ=255)
                                                                                  gccgTACAGAAAACG‑CCCCCCGCACGGCTGTCCATAACGAAGTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCctgct             <  2:1157883/97‑1 (MQ=255)
                                                                                    aGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTgg           >  2:2418686/1‑100 (MQ=255)
                                                                                    aGTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTgg           <  1:329647/100‑1 (MQ=255)
                                                                                     gTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGt          >  2:501318/1‑100 (MQ=255)
                                                                                     gTACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGt          <  2:796305/100‑1 (MQ=255)
                                                                                      tACAGAAAACGCCCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGt          <  1:1895192/100‑1 (MQ=255)
                                                                                      tACAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGtt         <  2:941308/100‑1 (MQ=255)
                                                                                       aCAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTc        <  2:1872247/100‑1 (MQ=255)
                                                                                        cAGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCa       >  2:2302408/1‑100 (MQ=255)
                                                                                         aGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCAc      >  2:807744/1‑100 (MQ=255)
                                                                                         aGAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCAc      >  1:22914/1‑100 (MQ=255)
                                                                                          gAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACg     >  1:2137752/1‑100 (MQ=255)
                                                                                          gAAAACG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACg     <  2:291152/100‑1 (MQ=255)
                                                                                            aaaCG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACGAt   >  2:805179/1‑100 (MQ=255)
                                                                                            aaaCG‑CCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACGAt   >  1:2262845/1‑100 (MQ=255)
                                                                                             aaCG‑TCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACGATc  >  2:1782745/1‑100 (MQ=255)
                                                                                                  |                                                                                               
ATTCTGCTGCCACTGCTGCTGAACATCGGTTTGCTGGCGCTGGCGTTATTTGGCTATACCACATTCCGCCATTTCTCCAGCCGCAGTACAGAAAACG‑TCCCCAGCACGGCGGTCAATAACGAATTGCGCATTTTACGGGCAATCAATGAAGAGATAGTCTCACTGCTGCCGCTCGGCCTGCTGGTTCACGATC  >  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: ‑