Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,792,242:1 +A 17.5% intergenic (‑222/+25) ydiV ← / ← nlpC anti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,792,2421.A17.5% 241.1 / 58.0 80intergenic (‑222/+25)ydiV/nlpCanti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (29/37);  new base A (8/6);  total (37/43)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

TTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCT  >  NC_000913/1792147‑1792341
                                                                                                |                                                                                                   
tttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat‑catc                                                                                                 >  2:1727789/1‑100 (MQ=255)
tttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat‑catc                                                                                                 >  2:548826/1‑100 (MQ=255)
 ttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTGTGTATGCAAAATACTACGGGAGTGAGGTGTATATGGACGGCGATGTAT‑CATCg                                                                                                <  2:1365356/100‑1 (MQ=255)
  tttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGc                                                                                               >  1:1695813/1‑100 (MQ=255)
  tttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGc                                                                                               >  1:264292/1‑100 (MQ=255)
   ttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGcc                                                                                              <  1:1074427/100‑1 (MQ=255)
   ttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGcc                                                                                              >  1:2255049/1‑100 (MQ=255)
    tAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCg                                                                                             <  2:1932698/100‑1 (MQ=255)
    tAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCg                                                                                             >  2:502121/1‑100 (MQ=255)
     aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGt                                                                                            <  2:330929/100‑1 (MQ=39)
     aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGt                                                                                            <  1:1757939/100‑1 (MQ=39)
      gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATA‑CATCGCCGTc                                                                                           <  1:31147/100‑1 (MQ=255)
      gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTc                                                                                           <  1:1585326/100‑1 (MQ=39)
       aTTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTcc                                                                                          <  2:1863235/100‑1 (MQ=39)
       aTTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTcc                                                                                          <  1:1101720/100‑1 (MQ=39)
        ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCACCGTCCa                                                                                         <  1:1715673/100‑1 (MQ=38)
         tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCAt                                                                                        <  1:295541/100‑1 (MQ=255)
         tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCAt                                                                                        <  2:1408100/100‑1 (MQ=255)
         tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCAt                                                                                        >  2:1724742/1‑100 (MQ=255)
             cTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCAttttt                                                                                    >  2:2165128/1‑100 (MQ=255)
             cTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCAttttt                                                                                    <  2:1571455/100‑1 (MQ=255)
              tGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTc                                                                                   <  1:1913327/100‑1 (MQ=255)
                 ttATAAACGAATAGTCCATTTACGGAATTTTTGTCTGCAAAAGACTACTGTATTCAGTAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGAc                                                                                <  2:2294958/100‑1 (MQ=255)
                   ataAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACAt                                                                              >  1:1732980/1‑100 (MQ=255)
                    taAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATa                                                                             >  2:2379422/1‑100 (MQ=255)
                     aaaCGAATAATCCATTTACGGACTTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGGCGGCGATGTAT‑CATCCCCGTCCATTTTTCGACATAc                                                                            >  2:617149/1‑100 (MQ=255)
                     aaaCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATAc                                                                            <  2:2134762/100‑1 (MQ=255)
                      aaCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACt                                                                           >  1:1821036/1‑100 (MQ=255)
                      aaCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGAAGATGTAT‑TATCGACGTCCATTTTTCGACATACt                                                                           >  2:2181058/1‑100 (MQ=255)
                       aCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTa                                                                          <  2:2319108/100‑1 (MQ=255)
                          aataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGAt                                                                       <  1:25534/100‑1 (MQ=255)
                                 cATTTACGGAATTTTTGTTTGTAAAATACTACTGTATTCAGGAGTAAATGTACGTCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGtctc                                                                >  2:2280470/1‑100 (MQ=255)
                                    ttACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGcc                                                             <  2:2001188/100‑1 (MQ=255)
                                     tACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCt                                                            >  2:112233/1‑100 (MQ=255)
                                        ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGcc                                                         <  2:1353109/100‑1 (MQ=255)
                                               ttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAt                                                   <  2:1257640/99‑1 (MQ=255)
                                                 gTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAtttt                                                <  1:1817040/100‑1 (MQ=255)
                                                  tCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAttttt                                               <  2:1236438/100‑1 (MQ=255)
                                                   cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                              <  2:1360811/100‑1 (MQ=255)
                                                     gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                            <  1:2269723/100‑1 (MQ=255)
                                                     gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                            <  2:658675/100‑1 (MQ=255)
                                                          aTACTCCTGTATTCAGGAGTAAATGGCAGGCTATGTAT‑CATCGTCGTCCATTTTTCGAGCTACGACATTCGTCCCGCTTGCCAGAAAATTTTCTGCCa                                         <  2:2372344/98‑1 (MQ=255)
                                                             ctactGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGa                                    >  1:2074422/1‑100 (MQ=255)
                                                               actGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACa                                  >  2:399689/1‑100 (MQ=255)
                                                               actGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACa                                  >  2:1221818/1‑100 (MQ=255)
                                                                ctGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACAt                                 >  1:382813/1‑100 (MQ=255)
                                                                ctGTATCCAGGAGTAAATGGATGGCGCTGTAT‑CAGCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCGGAAATTTTTCTGCCAATAg                                     <  2:1463531/96‑1 (MQ=255)
                                                                 tGTATTCAGGAGTAAATGGACAGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACAtt                                <  2:1719990/100‑1 (MQ=255)
                                                                  gTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTa                               >  1:900602/1‑100 (MQ=255)
                                                                   tATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTAt                              <  1:502123/100‑1 (MQ=255)
                                                                   tATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTAt                              <  2:13334/100‑1 (MQ=255)
                                                                    aTTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATc                             >  1:447685/1‑100 (MQ=255)
                                                                    aTTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATc                             <  1:2036171/100‑1 (MQ=255)
                                                                     ttCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCa                            >  1:339211/1‑100 (MQ=255)
                                                                     ttCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCa                            >  2:354302/1‑100 (MQ=255)
                                                                       cAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAg                          >  2:1710112/1‑100 (MQ=255)
                                                                       cAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAg                          >  1:810957/1‑100 (MQ=255)
                                                                        aGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGt                         <  1:2058393/100‑1 (MQ=255)
                                                                         ggAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                        <  1:730974/100‑1 (MQ=255)
                                                                         ggAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                        <  1:1033928/100‑1 (MQ=255)
                                                                           aGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAg                      >  2:1397472/1‑100 (MQ=255)
                                                                            gTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAgg                     <  2:1085900/100‑1 (MQ=255)
                                                                             tAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGa                    <  2:311976/100‑1 (MQ=255)
                                                                             tAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGa                    <  2:197534/100‑1 (MQ=255)
                                                                               aaTGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAc                  <  1:1068702/100‑1 (MQ=255)
                                                                               aaTGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAc                  >  2:1218487/1‑100 (MQ=255)
                                                                                 tGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGc                <  2:852339/100‑1 (MQ=255)
                                                                                  ggACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa               <  1:1201931/100‑1 (MQ=255)
                                                                                  ggACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa               <  1:194395/100‑1 (MQ=255)
                                                                                  ggACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATCTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa               >  1:177080/1‑100 (MQ=255)
                                                                                  ggACGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa               >  1:382075/1‑100 (MQ=39)
                                                                                   gACGGCGATGTAT‑CATCGCCGTCCATTTTTTTATATATTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  2:5389/1‑100 (MQ=255)
                                                                                   gACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  1:2372957/1‑100 (MQ=255)
                                                                                   gACGGCGATG‑ATACATCGCCGTCCATTTTTCTACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  1:715899/1‑100 (MQ=38)
                                                                                   gACGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  1:2236596/1‑100 (MQ=39)
                                                                                   gACGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  1:1222050/1‑100 (MQ=39)
                                                                                   gACGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  1:1848807/1‑100 (MQ=39)
                                                                                   gACGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  1:2000857/1‑100 (MQ=39)
                                                                                   gACGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCAAGAAATTTTTTTGAGAATAGACATTATCAAGTGAGGAACGAAt              >  2:1468068/1‑100 (MQ=37)
                                                                                    aCGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATc             <  2:2255029/100‑1 (MQ=255)
                                                                                    aCGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATc             >  1:33839/1‑100 (MQ=39)
                                                                                             tAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTg    >  2:2268486/1‑100 (MQ=255)
                                                                                              at‑catcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGc   <  1:112234/100‑1 (MQ=255)
                                                                                              at‑catcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACACTATCAAGTGAGGAACGCATCACTCCCTTGc   >  2:2294534/1‑100 (MQ=255)
                                                                                               t‑catcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCt  <  2:1732966/100‑1 (MQ=255)
                                                                                                |                                                                                                   
TTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCT  >  NC_000913/1792147‑1792341

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: ‑