Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,792,242:1 +A 16.4% 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.A16.4% 398.6 / 86.2 122intergenic (‑222/+25)ydiV/nlpCanti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (56/46);  new base A (9/11);  total (65/57)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.14e-01

Only 100 of 131 total aligned reads displayed.

CCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGC  >  NC_000913/1792143‑1792340
                                                                                                    |                                                                                                  
ccTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat                                                                                                     >  2:3255438/1‑100 (MQ=255)
 cTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat‑c                                                                                                   >  2:1669047/1‑100 (MQ=255)
  tttttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGCCGATGTat‑ca                                                                                                  >  1:2445762/1‑100 (MQ=255)
   ttttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat‑cat                                                                                                 >  1:2495032/1‑100 (MQ=255)
    tttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat‑catc                                                                                                <  2:2868191/100‑1 (MQ=255)
     ttttAGATTATTCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCg                                                                                               <  1:649401/100‑1 (MQ=255)
     ttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCg                                                                                               <  2:2786930/100‑1 (MQ=255)
        tAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGGCTGCAAAATACTACTGTATTCAGGATTAAATGTACGTCGACGTAT‑CATCGCCg                                                                                            <  1:1310971/100‑1 (MQ=255)
         aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGt                                                                                           >  2:1889634/1‑100 (MQ=255)
         aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGt                                                                                           >  2:3255075/1‑100 (MQ=255)
         aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATA‑CATCGCCGt                                                                                           <  1:947569/100‑1 (MQ=39)
         aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGt                                                                                           <  2:3570963/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTc                                                                                          <  2:2561662/100‑1 (MQ=255)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATA‑CATCGCCGTc                                                                                          <  1:1071614/100‑1 (MQ=255)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATA‑CATCGCCGTc                                                                                          <  1:953752/100‑1 (MQ=255)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTc                                                                                          <  2:3154183/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTc                                                                                          <  2:844668/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTc                                                                                          <  2:587097/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTc                                                                                          <  1:890999/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTc                                                                                          <  1:890987/100‑1 (MQ=39)
           aTTATCCTTATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTcc                                                                                         <  1:1371615/100‑1 (MQ=38)
            ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTCCa                                                                                        <  2:2694168/100‑1 (MQ=38)
            ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTCCa                                                                                        <  1:3073009/100‑1 (MQ=38)
            ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATCGACGGCGATGTAT‑CATCGCGGTc                                                                                          >  1:2408426/1‑98 (MQ=255)
                ccTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCAtttt                                                                                    >  1:2948500/1‑100 (MQ=255)
                ccTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCAtttt                                                                                    >  2:3680533/1‑100 (MQ=255)
                     ttATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGAc                                                                               <  2:1615453/100‑1 (MQ=255)
                     ttATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGAc                                                                               <  2:391720/100‑1 (MQ=255)
                        taAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTCCAtttactcctgaa                                                                            >  1:3572134/1‑91 (MQ=25)
                         aaaCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTCCAtttactcctgaat                                                                           <  2:3572162/100‑11 (MQ=21)
                          aaCGAATAATCCATTTCCGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAATTGATCGCCGATGTAT‑CATCGCTGACCATTTTTAGACATACt                                                                          >  1:2859554/1‑100 (MQ=255)
                          aaCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACt                                                                          <  1:1780116/100‑1 (MQ=255)
                            cGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACGTACTAg                                                                        <  2:1594091/100‑1 (MQ=255)
                            cGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAg                                                                        <  2:3584229/100‑1 (MQ=255)
                              aataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGAt                                                                      <  1:3585792/100‑1 (MQ=255)
                              aataatCCATTAACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGAt                                                                      >  1:1958549/1‑100 (MQ=255)
                               ataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGAt                                                                      <  1:3346232/99‑1 (MQ=255)
                                    ccATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGtct                                                                >  2:3486147/1‑100 (MQ=255)
                                    ccATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGtct                                                                >  1:2984741/1‑100 (MQ=255)
                                        ttACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGcc                                                            <  1:2244383/100‑1 (MQ=255)
                                          aCGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTg                                                          >  2:574643/1‑100 (MQ=255)
                                           cGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTAGACATAATAGATTAGTCTCTCCTGc                                                         >  1:660782/1‑100 (MQ=255)
                                            ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGTAAGGCGATGTAT‑CATCGCCGTACATTTTTCGAGATACTAGATTAGTGTTGcc                                                            >  1:3213949/1‑96 (MQ=255)
                                            ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGcc                                                        <  2:2478363/100‑1 (MQ=255)
                                             gAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCa                                                       >  1:3098303/1‑100 (MQ=255)
                                             gAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCa                                                       >  2:2573614/1‑100 (MQ=255)
                                               aTTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGa                                                     <  2:84508/100‑1 (MQ=255)
                                                 ttttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGaaa                                                   <  2:2300165/100‑1 (MQ=255)
                                                   ttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAtt                                                 <  1:1901443/100‑1 (MQ=255)
                                                      tCTGCAAGATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAttt                                                >  2:1954484/1‑98 (MQ=255)
                                                      tCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAttttt                                              <  2:1018940/100‑1 (MQ=255)
                                                       cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                             >  1:2052342/1‑100 (MQ=255)
                                                        tGCAAAATACTCCTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGCCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                             <  2:1677751/99‑1 (MQ=255)
                                                        tGCAAAATACTACTGTATTCAGGAGTAAATGGACTGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCt                                            >  1:2461919/1‑100 (MQ=255)
                                                        tGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCt                                            >  2:101436/1‑100 (MQ=255)
                                                        tGCAAAATACTACTGTATTCAGGAGTAAATGGAAGGCGATGTAT‑CATCGACGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAACTTTTTCt                                            >  1:2984923/1‑100 (MQ=255)
                                                         gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                           >  1:2409729/1‑100 (MQ=255)
                                                         gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                           >  2:1053010/1‑100 (MQ=255)
                                                          cAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGc                                          >  2:2141230/1‑100 (MQ=255)
                                                                  tactGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGAc                                  <  2:1767668/100‑1 (MQ=255)
                                                                   actGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATCCTAGATTCGTCTCGCCTGCAAGAAATTTTTCTGCCAATAGACa                                 >  1:570802/1‑100 (MQ=255)
                                                                   actGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGAc                                  >  1:1073869/1‑99 (MQ=255)
                                                                   actGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACa                                 >  2:913127/1‑100 (MQ=255)
                                                                   actGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACa                                 <  1:311685/100‑1 (MQ=255)
                                                                   actGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACa                                 <  1:1056828/100‑1 (MQ=255)
                                                                      gTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTa                              >  2:928193/1‑100 (MQ=255)
                                                                         ttCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTTGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCa                           <  2:2636856/100‑1 (MQ=255)
                                                                           cAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAg                         <  2:707330/100‑1 (MQ=255)
                                                                             tccgatcttGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                       <  2:1210713/92‑1 (MQ=255)
                                                                             ggAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTCTCAAGTg                       >  1:1090647/1‑100 (MQ=255)
                                                                             ggAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                       >  1:3268712/1‑100 (MQ=255)
                                                                             ggAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                       <  1:1538671/100‑1 (MQ=255)
                                                                             ggAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATATTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                       >  2:2635050/1‑100 (MQ=255)
                                                                              gAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGa                      >  2:1336386/1‑100 (MQ=255)
                                                                              gAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGa                      >  2:2260619/1‑100 (MQ=255)
                                                                              cgaaaaaTGGACGGCGATGATA‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACCTTATCAAGTGa                      >  1:1203772/5‑100 (MQ=37)
                                                                               aGTAAATTGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAg                     <  2:158970/100‑1 (MQ=255)
                                                                                gTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAgg                    <  2:1313562/100‑1 (MQ=255)
                                                                                gTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAgg                    <  2:3336403/100‑1 (MQ=255)
                                                                                  aaaTGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGaa                  >  1:2874417/1‑100 (MQ=255)
                                                                                  aaaTGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAATTTTCTGCCAATAGACATTATCAAGTGAGGaa                  >  2:3369624/1‑100 (MQ=255)
                                                                                   aaTGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAc                 <  1:1378658/100‑1 (MQ=255)
                                                                                     tGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAgatcgga               >  1:1210707/1‑94 (MQ=255)
                                                                                     tGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGc               >  2:1126331/1‑100 (MQ=255)
                                                                                      ngACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa              >  1:1828843/2‑100 (MQ=255)
                                                                                      ggACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa              <  1:1710286/100‑1 (MQ=255)
                                                                                      ggACGGCGATG‑ATACATCGCCGTCCATTTTTTTTCTTTTTATATTCGTCTCGCCTGCCAGAAATTTTTCTCCCTTTAGACATTATCAAGTTACGAACGCa              >  1:3409014/1‑100 (MQ=16)
                                                                                      ggACGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa              >  2:105738/1‑100 (MQ=39)
                                                                                       gACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATGAAGTGAGGAACGCAt             >  2:3667181/1‑100 (MQ=255)
                                                                                       gACGGCGATGATA‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt             >  1:1449578/1‑100 (MQ=39)
                                                                                       gACGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt             >  2:905340/1‑100 (MQ=39)
                                                                                       gACGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATGTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt             >  2:2626317/1‑100 (MQ=38)
                                                                                        aCGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATc            <  2:1573513/100‑1 (MQ=255)
                                                                                        aCGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTGGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGAGAGGAACGCATc            >  2:2182506/1‑100 (MQ=38)
                                                                                         cGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCa           <  1:3680505/100‑1 (MQ=255)
                                                                                            cGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTc        >  2:622742/1‑100 (MQ=255)
                                                                                               tGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCt     >  2:423600/1‑100 (MQ=255)
                                                                                               tGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCt     <  2:2874436/100‑1 (MQ=255)
                                                                                                  at‑catcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGc  <  2:1305746/100‑1 (MQ=255)
                                                                                                   t‑catcGTCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGc  <  1:3022780/99‑1 (MQ=255)
                                                                                                    |                                                                                                  
CCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGC  >  NC_000913/1792143‑1792340

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