Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,792,242:1 +A 15.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.A15.5% 263.1 / 50.6 84intergenic (‑222/+25)ydiV/nlpCanti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (30/41);  new base A (4/9);  total (34/50)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.84e-01

CCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCT  >  NC_000913/1792143‑1792341
                                                                                                    |                                                                                                   
ccTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat                                                                                                      >  2:2363944/1‑100 (MQ=255)
   ttttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat‑cat                                                                                                  <  2:1713054/100‑1 (MQ=255)
      tttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGc                                                                                               <  2:29180/100‑1 (MQ=255)
       ttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGcc                                                                                              <  2:2529955/100‑1 (MQ=255)
        tAGATTATCCTGATTGTAAACGAATAATGCAGTTACGGAATTTTTGTCTGCAAGGTACGACTGTATTCAGGAGTAAATGGACGGCTGTGTAT‑CATTGc                                                                                               >  1:894286/1‑98 (MQ=255)
        tAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCg                                                                                             <  2:952551/100‑1 (MQ=255)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTc                                                                                           >  2:2486506/1‑100 (MQ=255)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATA‑CATCGCCGTc                                                                                           <  2:1688256/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTc                                                                                           <  1:93624/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTc                                                                                           <  2:125086/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTc                                                                                           <  2:506962/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTc                                                                                           <  1:1658661/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTc                                                                                           <  2:743844/100‑1 (MQ=39)
           aTTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTcc                                                                                          <  2:734492/100‑1 (MQ=255)
           aTTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTcc                                                                                          <  1:1445520/100‑1 (MQ=39)
            ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTCCa                                                                                         <  1:1313993/100‑1 (MQ=38)
             tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCAt                                                                                        <  2:210316/100‑1 (MQ=255)
             tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCAt                                                                                        <  1:201812/100‑1 (MQ=255)
             tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCAt                                                                                        >  2:620525/1‑100 (MQ=255)
              aTCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTCCAtt                                                                                       <  2:302076/100‑1 (MQ=38)
                  tGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTc                                                                                   <  2:1629672/100‑1 (MQ=255)
                     ttATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTATTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGAc                                                                                <  1:18549/100‑1 (MQ=255)
                      tataAACGAATAATCCATTTACGGGATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACa                                                                               <  1:912203/100‑1 (MQ=255)
                      tataAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTCCAtttactcctg                                                                               <  2:890956/100‑8 (MQ=25)
                        taAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGCAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATa                                                                             <  1:516534/100‑1 (MQ=255)
                            cGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAg                                                                         <  2:1231198/100‑1 (MQ=255)
                              aataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGAt                                                                       <  1:1165612/100‑1 (MQ=255)
                                 aatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCg                                                                    <  2:2273904/100‑1 (MQ=255)
                                 aatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCg                                                                    <  1:1060400/100‑1 (MQ=255)
                                   tCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGtc                                                                  >  1:1543822/1‑100 (MQ=255)
                                      aTTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCg                                                               <  2:507772/100‑1 (MQ=255)
                                        ttACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGcc                                                             <  1:499737/100‑1 (MQ=255)
                                         tACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCt                                                            >  1:960364/1‑100 (MQ=255)
                                          aCGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTg                                                           >  2:2077149/1‑100 (MQ=255)
                                          aCGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTg                                                           >  1:699902/1‑100 (MQ=255)
                                           cGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGc                                                          >  2:1788078/1‑100 (MQ=255)
                                           cGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGAGATGTAT‑CATCACCGTCCATTTTTCGCCATACTAGATTCGTCTCGCCTGc                                                          >  2:1164467/1‑100 (MQ=255)
                                            ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGcc                                                         >  1:1294149/1‑100 (MQ=255)
                                             gAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCa                                                        >  1:1917842/1‑100 (MQ=255)
                                             gAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCa                                                        <  2:736460/100‑1 (MQ=255)
                                                  tttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAt                                                   >  2:641112/1‑100 (MQ=255)
                                                  tttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAt                                                   <  1:746593/100‑1 (MQ=255)
                                                      tCTGCAAAATACTACTTTATTCAGTAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAttttt                                               <  2:1507815/100‑1 (MQ=255)
                                                       cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                              <  1:2077135/100‑1 (MQ=255)
                                                       cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                              <  1:1275494/100‑1 (MQ=255)
                                                       cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                              <  2:1788317/100‑1 (MQ=255)
                                                       cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                              >  1:59079/1‑100 (MQ=255)
                                                       cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                              >  1:281805/1‑100 (MQ=255)
                                                        tGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCt                                             >  1:1140232/1‑100 (MQ=255)
                                                         gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                            >  2:975499/1‑100 (MQ=255)
                                                         gCAAAATACTACTGTATTCAGGAGTAAATGGACGGAGATGTAT‑CATAGACGTCCATTTTTCGACATACTAGATTCGTCTAGACTGCCAGAAATTTTTCTg                                            >  1:1311920/1‑100 (MQ=255)
                                                          cAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGc                                           <  1:200602/100‑1 (MQ=255)
                                                          cAAAATAATACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGc                                           <  1:1095718/100‑1 (MQ=255)
                                                            aaaTACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCACCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTCCCAGAAATTTTTCTGCCa                                         >  1:561623/1‑100 (MQ=255)
                                                                actactGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCTAGAAATTTTTCTGCCAATAg                                     <  1:1586494/100‑1 (MQ=255)
                                                                actactGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAg                                     <  2:536707/100‑1 (MQ=255)
                                                                actactGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAg                                     <  2:1195073/100‑1 (MQ=255)
                                                                actactGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAg                                     <  1:451605/100‑1 (MQ=255)
                                                                   actGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATATACa                                  >  1:1926039/1‑100 (MQ=255)
                                                                   actGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACa                                  <  1:333873/100‑1 (MQ=255)
                                                                      gTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTAGCCTGCCAGAAATTTTTCTGCCAATAGACATTa                               >  1:2580546/1‑100 (MQ=255)
                                                                        aTTCAGGAGTAAAGGGACGGCCATGTAT‑CATCTCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATc                             <  2:356420/100‑1 (MQ=255)
                                                                         ttCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCa                            <  1:620522/100‑1 (MQ=255)
                                                                           cAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAg                          >  2:1319575/1‑100 (MQ=255)
                                                                           cAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAg                          <  1:1319568/100‑1 (MQ=255)
                                                                            aGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGt                         <  1:1104322/100‑1 (MQ=255)
                                                                             ggAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCATCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                        <  1:1981352/100‑1 (MQ=255)
                                                                               aGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAg                      >  2:2076458/1‑100 (MQ=255)
                                                                                gTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGACTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                        >  1:2323295/1‑97 (MQ=255)
                                                                                  aaaTGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGaa                   >  2:1948647/1‑100 (MQ=255)
                                                                                   aaTGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAc                  <  2:784953/100‑1 (MQ=255)
                                                                                    aTGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACg                 >  2:887065/1‑100 (MQ=255)
                                                                                    aTGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACg                 <  2:1543831/100‑1 (MQ=255)
                                                                                      ggACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa               <  1:921261/100‑1 (MQ=255)
                                                                                      ggACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa               >  1:168502/1‑100 (MQ=255)
                                                                                       gACGGCGATGTAT‑CATCGCCGTCCATTTTTTTTTTTATTTGATTCGTCTCGCCTGCAAGCAATTTTTCTACCACTAGACATTTTCAAGTGGGGAACGCAt              >  1:68167/1‑100 (MQ=255)
                                                                                       gACGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTTTTTTTTTTTTTATTATCAAGTGAGGAACGCAt              >  1:1714276/1‑100 (MQ=19)
                                                                                       gACGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  2:770618/1‑100 (MQ=39)
                                                                                       gACGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  1:1443931/1‑100 (MQ=39)
                                                                                       gACGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  2:1889028/1‑100 (MQ=39)
                                                                                         cGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCa            <  2:1427363/100‑1 (MQ=255)
                                                                                         cGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCa            >  2:394640/1‑100 (MQ=255)
                                                                                         cGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCa            >  1:913928/1‑100 (MQ=255)
                                                                                           gCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGCGAGGAACGcaccac           >  1:1266952/1‑99 (MQ=255)
                                                                                               tGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCt      <  1:1345217/100‑1 (MQ=255)
                                                                                                gTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCtt     <  1:280420/100‑1 (MQ=255)
                                                                                                  at‑catcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGc   >  2:636891/1‑100 (MQ=255)
                                                                                                   t‑catcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCt  <  1:1874536/100‑1 (MQ=255)
                                                                                                   t‑catcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCt  <  1:2234422/100‑1 (MQ=255)
                                                                                                   t‑catcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCt  >  2:1202891/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCT  >  NC_000913/1792143‑1792341

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