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

CCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCT  >  NC_000913/1792143‑1792341
                                                                                                    |                                                                                                   
ccTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat                                                                                                      >  2:2363944/1‑100 (MQ=255)
   ttttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTatcat                                                                                                  <  2:1713054/100‑1 (MQ=255)
      tttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGc                                                                                               <  2:29180/100‑1 (MQ=255)
       ttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGcc                                                                                              <  2:2529955/100‑1 (MQ=255)
        tAGATTATCCTGATTGTAAACGAATAATGCAGTTACGGAATTTTTGTCTGCAAGGTACGACTGTATTCAGGAGTAAATGGACGGCTGTGTATCATTGc                                                                                               >  1:894286/1‑98 (MQ=255)
        tAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCg                                                                                             <  2:952551/100‑1 (MQ=255)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTc                                                                                           >  2:2486506/1‑100 (MQ=255)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                           <  2:1688256/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                           <  1:93624/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                           <  2:125086/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                           <  2:506962/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                           <  1:1658661/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                           <  2:743844/100‑1 (MQ=39)
           aTTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTcc                                                                                          <  2:734492/100‑1 (MQ=255)
           aTTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTcc                                                                                          <  1:1445520/100‑1 (MQ=39)
            ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTCCa                                                                                         <  1:1313993/100‑1 (MQ=38)
             tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAt                                                                                        <  2:210316/100‑1 (MQ=255)
             tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAt                                                                                        <  1:201812/100‑1 (MQ=255)
             tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAt                                                                                        >  2:620525/1‑100 (MQ=255)
              aTCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTCCAtt                                                                                       <  2:302076/100‑1 (MQ=38)
                  tGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTc                                                                                   <  2:1629672/100‑1 (MQ=255)
                     ttATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTATTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGAc                                                                                <  1:18549/100‑1 (MQ=255)
                      tataAACGAATAATCCATTTACGGGATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACa                                                                               <  1:912203/100‑1 (MQ=255)
                      tataAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTCCAtttactcctg                                                                               <  2:890956/100‑8 (MQ=25)
                        taAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGCAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATa                                                                             <  1:516534/100‑1 (MQ=255)
                            cGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAg                                                                         <  2:1231198/100‑1 (MQ=255)
                              aataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGAt                                                                       <  1:1165612/100‑1 (MQ=255)
                                 aatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCg                                                                    <  2:2273904/100‑1 (MQ=255)
                                 aatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCg                                                                    <  1:1060400/100‑1 (MQ=255)
                                   tCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGtc                                                                  >  1:1543822/1‑100 (MQ=255)
                                      aTTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCg                                                               <  2:507772/100‑1 (MQ=255)
                                        ttACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGcc                                                             <  1:499737/100‑1 (MQ=255)
                                         tACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCt                                                            >  1:960364/1‑100 (MQ=255)
                                          aCGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTg                                                           >  2:2077149/1‑100 (MQ=255)
                                          aCGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTg                                                           >  1:699902/1‑100 (MQ=255)
                                           cGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGc                                                          >  2:1788078/1‑100 (MQ=255)
                                           cGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGAGATGTATCATCACCGTCCATTTTTCGCCATACTAGATTCGTCTCGCCTGc                                                          >  2:1164467/1‑100 (MQ=255)
                                            ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGcc                                                         >  1:1294149/1‑100 (MQ=255)
                                             gAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCa                                                        >  1:1917842/1‑100 (MQ=255)
                                             gAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCa                                                        <  2:736460/100‑1 (MQ=255)
                                                  tttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAt                                                   >  2:641112/1‑100 (MQ=255)
                                                  tttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAt                                                   <  1:746593/100‑1 (MQ=255)
                                                      tCTGCAAAATACTACTTTATTCAGTAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAttttt                                               <  2:1507815/100‑1 (MQ=255)
                                                       cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                              <  1:2077135/100‑1 (MQ=255)
                                                       cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                              <  1:1275494/100‑1 (MQ=255)
                                                       cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                              <  2:1788317/100‑1 (MQ=255)
                                                       cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                              >  1:59079/1‑100 (MQ=255)
                                                       cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                              >  1:281805/1‑100 (MQ=255)
                                                        tGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCt                                             >  1:1140232/1‑100 (MQ=255)
                                                         gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                            >  2:975499/1‑100 (MQ=255)
                                                         gCAAAATACTACTGTATTCAGGAGTAAATGGACGGAGATGTATCATAGACGTCCATTTTTCGACATACTAGATTCGTCTAGACTGCCAGAAATTTTTCTg                                            >  1:1311920/1‑100 (MQ=255)
                                                          cAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGc                                           <  1:200602/100‑1 (MQ=255)
                                                          cAAAATAATACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGc                                           <  1:1095718/100‑1 (MQ=255)
                                                            aaaTACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCACCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTCCCAGAAATTTTTCTGCCa                                         >  1:561623/1‑100 (MQ=255)
                                                                actactGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCTAGAAATTTTTCTGCCAATAg                                     <  1:1586494/100‑1 (MQ=255)
                                                                actactGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAg                                     <  2:536707/100‑1 (MQ=255)
                                                                actactGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAg                                     <  2:1195073/100‑1 (MQ=255)
                                                                actactGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAg                                     <  1:451605/100‑1 (MQ=255)
                                                                   actGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATATACa                                  >  1:1926039/1‑100 (MQ=255)
                                                                   actGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACa                                  <  1:333873/100‑1 (MQ=255)
                                                                      gTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTAGCCTGCCAGAAATTTTTCTGCCAATAGACATTa                               >  1:2580546/1‑100 (MQ=255)
                                                                        aTTCAGGAGTAAAGGGACGGCCATGTATCATCTCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATc                             <  2:356420/100‑1 (MQ=255)
                                                                         ttCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCa                            <  1:620522/100‑1 (MQ=255)
                                                                           cAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAg                          >  2:1319575/1‑100 (MQ=255)
                                                                           cAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAg                          <  1:1319568/100‑1 (MQ=255)
                                                                            aGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGt                         <  1:1104322/100‑1 (MQ=255)
                                                                             ggAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCATCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                        <  1:1981352/100‑1 (MQ=255)
                                                                               aGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAg                      >  2:2076458/1‑100 (MQ=255)
                                                                                gTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGACTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                        >  1:2323295/1‑97 (MQ=255)
                                                                                  aaaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGaa                   >  2:1948647/1‑100 (MQ=255)
                                                                                   aaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAc                  <  2:784953/100‑1 (MQ=255)
                                                                                    aTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACg                 >  2:887065/1‑100 (MQ=255)
                                                                                    aTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACg                 <  2:1543831/100‑1 (MQ=255)
                                                                                      ggACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa               <  1:921261/100‑1 (MQ=255)
                                                                                      ggACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa               >  1:168502/1‑100 (MQ=255)
                                                                                       gACGGCGATGTATCATCGCCGTCCATTTTTTTTTTTATTTGATTCGTCTCGCCTGCAAGCAATTTTTCTACCACTAGACATTTTCAAGTGGGGAACGCAt              >  1:68167/1‑100 (MQ=255)
                                                                                       gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTTTTTTTTTTTTTATTATCAAGTGAGGAACGCAt              >  1:1714276/1‑100 (MQ=19)
                                                                                       gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  2:770618/1‑100 (MQ=39)
                                                                                       gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  1:1443931/1‑100 (MQ=39)
                                                                                       gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt              >  2:1889028/1‑100 (MQ=39)
                                                                                         cGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCa            <  2:1427363/100‑1 (MQ=255)
                                                                                         cGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCa            >  2:394640/1‑100 (MQ=255)
                                                                                         cGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCa            >  1:913928/1‑100 (MQ=255)
                                                                                           gCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGCGAGGAACGcaccac           >  1:1266952/1‑99 (MQ=255)
                                                                                               tGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCt      <  1:1345217/100‑1 (MQ=255)
                                                                                                gTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCtt     <  1:280420/100‑1 (MQ=255)
                                                                                                  atcatcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGc   >  2:636891/1‑100 (MQ=255)
                                                                                                   tcatcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCt  <  1:1874536/100‑1 (MQ=255)
                                                                                                   tcatcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCt  <  1:2234422/100‑1 (MQ=255)
                                                                                                   tcatcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCt  >  2:1202891/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
CCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCT  >  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: