Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,101,039 (C)7→6 coding (585/720 nt) yggN ← DUF2884 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,101,0330C.96.2% 335.1 / 7.7 81coding (591/720 nt)yggNDUF2884 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base . (41/36);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCC  >  NC_000913/3100939‑3101128
                                                                                              |                                                                                               
cTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGcccccc                                                                                            >  2:806808/1‑100 (MQ=255)
 ttCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGcccccca                                                                                           >  1:1232995/1‑99 (MQ=255)
 ttCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGcccccca                                                                                           >  2:2188600/1‑99 (MQ=255)
  tCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGccccccag                                                                                          >  2:57358/1‑98 (MQ=255)
   ccAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGccccccagg                                                                                         <  1:2119186/100‑4 (MQ=255)
   ccAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGccccccagg                                                                                         <  2:415781/100‑4 (MQ=255)
   ccAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGccccccagg                                                                                         <  2:781100/100‑4 (MQ=255)
       agTCACAACGCGGCTACAAACATCTTTGCCAAACTGGTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTc                                                                                    <  1:2025347/100‑1 (MQ=255)
        gTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTcc                                                                                   >  1:1107902/1‑100 (MQ=255)
         tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTccc                                                                                  >  2:2844128/1‑100 (MQ=255)
         tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTccc                                                                                  <  2:2091604/100‑1 (MQ=255)
          cacaACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCa                                                                                 >  2:1315841/1‑100 (MQ=255)
             aaCGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGca                                                                              <  1:1931624/100‑1 (MQ=255)
              aCGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGcac                                                                             <  2:2089242/100‑1 (MQ=255)
                gcgGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACAt                                                                           >  2:1143825/1‑100 (MQ=255)
                  ggCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTc                                                                         >  1:1839392/1‑100 (MQ=255)
                  ggCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTc                                                                         <  1:2150922/100‑1 (MQ=255)
                   gCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCt                                                                        <  1:2787673/100‑1 (MQ=255)
                   gCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCt                                                                        <  1:2770968/100‑1 (MQ=255)
                   gCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCCGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCt                                                                        >  2:364779/1‑100 (MQ=255)
                     tACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGt                                                                      >  1:2833045/1‑100 (MQ=255)
                     tACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGt                                                                      >  2:534826/1‑100 (MQ=255)
                     tACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGt                                                                      <  1:1103047/100‑1 (MQ=255)
                     tACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATCGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGt                                                                      <  2:1778070/100‑1 (MQ=255)
                      aCAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTa                                                                     <  1:1018093/100‑1 (MQ=255)
                       cAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTaa                                                                    >  1:2159674/1‑100 (MQ=255)
                       cAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTaa                                                                    >  1:101245/1‑100 (MQ=255)
                         aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATg                                                                  >  1:2153921/1‑100 (MQ=255)
                         aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATg                                                                  >  2:788734/1‑100 (MQ=255)
                         aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATg                                                                  <  2:1633791/100‑1 (MQ=255)
                          aCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATgg                                                                 <  1:2095779/100‑1 (MQ=255)
                           cATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTTTTTTTTTTTTTTTTTGAGGATTGCGGCCGCTCAGGCTTTCCAGCACATTCTGTAAt                                                                   >  2:1800208/1‑97 (MQ=255)
                           cATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATggg                                                                <  1:806809/100‑1 (MQ=255)
                           cATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATggg                                                                <  2:2059794/100‑1 (MQ=255)
                            aTCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGt                                                               <  1:1464233/100‑1 (MQ=255)
                             tCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGtt                                                              >  1:673790/1‑100 (MQ=255)
                               tttGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAc                                                            >  1:229939/1‑100 (MQ=255)
                                   ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgcc                                                        <  1:1826494/100‑1 (MQ=255)
                                    cAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                       <  1:2225803/100‑1 (MQ=255)
                                     aaaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                       >  1:1499101/1‑99 (MQ=255)
                                     aaaCTGCTGGAAATCTTTTTCCTGCTTTTTCAACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACcgccgc                                                      >  1:1413201/1‑100 (MQ=255)
                                        ctgctgGAAATCTTTTTCCTGCTTTTTTTTTTTGGTTTGGATTGAGGATTGCAGCCTCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCttt                                                   >  2:212892/1‑100 (MQ=255)
                                        ctgctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCttt                                                   <  1:1889217/100‑1 (MQ=255)
                                          gctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTc                                                 <  1:2318557/100‑1 (MQ=255)
                                           ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCa                                                >  2:2726383/1‑100 (MQ=255)
                                            tgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAg                                               <  1:378852/100‑1 (MQ=255)
                                             ggAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGc                                              >  1:396946/1‑100 (MQ=255)
                                             ggAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGc                                              <  2:276479/100‑1 (MQ=255)
                                             ggAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGc                                              <  1:171608/100‑1 (MQ=255)
                                             ggAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGc                                              >  2:2440598/1‑100 (MQ=255)
                                              gAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGc                                              >  1:1867353/1‑100 (MQ=255)
                                               aaaTCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCAc                                            >  1:2741967/1‑100 (MQ=255)
                                                   cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCt                                        <  1:1479979/100‑1 (MQ=255)
                                                   cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCt                                        >  2:2270062/1‑100 (MQ=255)
                                                       ttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcg                                    <  2:189055/100‑1 (MQ=255)
                                                       ttCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcg                                    >  1:176764/1‑100 (MQ=255)
                                                          cTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCgcgc                                 >  2:2488850/1‑100 (MQ=255)
                                                          cTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCgcgc                                 <  1:378548/100‑1 (MQ=255)
                                                            gCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGccc                               >  1:2104224/1‑100 (MQ=255)
                                                            gCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGccc                               >  2:1866888/1‑100 (MQ=255)
                                                            gCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGccc                               >  2:1867017/1‑100 (MQ=255)
                                                             cTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCa                              <  1:876183/100‑1 (MQ=255)
                                                               ttttCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAtt                            >  1:1143915/1‑100 (MQ=255)
                                                                tttCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAttt                           >  2:67743/1‑100 (MQ=255)
                                                                   tcACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAt                        <  2:1570607/99‑1 (MQ=255)
                                                                    cACTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCAta                       >  1:1570616/1‑99 (MQ=255)
                                                                     aCTCGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTa                      >  2:2365268/1‑100 (MQ=255)
                                                                        cGGTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                   <  2:2741946/100‑1 (MQ=255)
                                                                          gTTTGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCt                 >  2:2827801/1‑100 (MQ=255)
                                                                           tttGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCt                 >  2:1427116/1‑100 (MQ=255)
                                                                           tttGGATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                <  2:919213/100‑1 (MQ=255)
                                                                            ttGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                <  1:1867609/100‑1 (MQ=255)
                                                                             tGGATTGAGGATTGCAGTCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTc              <  2:2234282/100‑1 (MQ=255)
                                                                              ggATTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTcc             <  1:534828/100‑1 (MQ=255)
                                                                                aTTGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTg           <  1:1335597/100‑1 (MQ=255)
                                                                                  tgatggttGCAGCCTCCCAGGCTTCCCAGCACACTCTGTAATGGGTTTTCGCCGCTTTTCGTCACCGCTTTCGCGCCCATTTCATTGATGCTGTCCTGTa         <  2:1781840/94‑1 (MQ=255)
                                                                                  tGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa         <  2:361882/100‑1 (MQ=255)
                                                                                  tGAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTa         >  2:746729/1‑100 (MQ=255)
                                                                                   gAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTaa        >  1:2263750/1‑100 (MQ=255)
                                                                                   gAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTaa        >  2:780492/1‑100 (MQ=255)
                                                                                   gAGGATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTaa        >  2:1152045/1‑100 (MQ=255)
                                                                                     ggATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTaaaa      >  1:2837421/1‑100 (MQ=255)
                                                                                     ggATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTaaaa      >  2:1473720/1‑100 (MQ=255)
                                                                                     ggATTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTaaaa      >  2:1845036/1‑100 (MQ=255)
                                                                                       aTTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAAtt    <  1:1857786/100‑1 (MQ=255)
                                                                                       aTTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAAtt    >  2:11278/1‑100 (MQ=255)
                                                                                       aTTGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAAtt    <  2:2335563/100‑1 (MQ=255)
                                                                                         tGCAGCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTcc  >  1:205206/1‑100 (MQ=255)
                                                                                              |                                                                                               
CTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCC  >  NC_000913/3100939‑3101128

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