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.

CCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGC  >  NC_000913/1792143‑1792340
                                                                                                    |                                                                                                  
ccTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat                                                                                                     >  2:3255438/1‑100 (MQ=255)
 cTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTatc                                                                                                   >  2:1669047/1‑100 (MQ=255)
  tttttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGCCGATGTatca                                                                                                  >  1:2445762/1‑100 (MQ=255)
   ttttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTatcat                                                                                                 >  1:2495032/1‑100 (MQ=255)
    tttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTatcatc                                                                                                <  2:2868191/100‑1 (MQ=255)
     ttttAGATTATTCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCg                                                                                               <  1:649401/100‑1 (MQ=255)
     ttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCg                                                                                               <  2:2786930/100‑1 (MQ=255)
        tAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGGCTGCAAAATACTACTGTATTCAGGATTAAATGTACGTCGACGTATCATCGCCg                                                                                            <  1:1310971/100‑1 (MQ=255)
         aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGt                                                                                           >  2:1889634/1‑100 (MQ=255)
         aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGt                                                                                           >  2:3255075/1‑100 (MQ=255)
         aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGt                                                                                           <  1:947569/100‑1 (MQ=39)
         aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGt                                                                                           <  2:3570963/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTc                                                                                          <  2:2561662/100‑1 (MQ=255)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                          <  1:1071614/100‑1 (MQ=255)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                          <  1:953752/100‑1 (MQ=255)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                          <  2:3154183/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                          <  2:844668/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                          <  2:587097/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                          <  1:890999/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                          <  1:890987/100‑1 (MQ=39)
           aTTATCCTTATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTcc                                                                                         <  1:1371615/100‑1 (MQ=38)
            ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTCCa                                                                                        <  2:2694168/100‑1 (MQ=38)
            ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTCCa                                                                                        <  1:3073009/100‑1 (MQ=38)
            ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATCGACGGCGATGTATCATCGCGGTc                                                                                          >  1:2408426/1‑98 (MQ=255)
                ccTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAtttt                                                                                    >  1:2948500/1‑100 (MQ=255)
                ccTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAtttt                                                                                    >  2:3680533/1‑100 (MQ=255)
                     ttATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGAc                                                                               <  2:1615453/100‑1 (MQ=255)
                     ttATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGAc                                                                               <  2:391720/100‑1 (MQ=255)
                        taAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTCCAtttactcctgaa                                                                            >  1:3572134/1‑91 (MQ=25)
                         aaaCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTCCAtttactcctgaat                                                                           <  2:3572162/100‑11 (MQ=21)
                          aaCGAATAATCCATTTCCGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAATTGATCGCCGATGTATCATCGCTGACCATTTTTAGACATACt                                                                          >  1:2859554/1‑100 (MQ=255)
                          aaCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACt                                                                          <  1:1780116/100‑1 (MQ=255)
                            cGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACGTACTAg                                                                        <  2:1594091/100‑1 (MQ=255)
                            cGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAg                                                                        <  2:3584229/100‑1 (MQ=255)
                              aataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGAt                                                                      <  1:3585792/100‑1 (MQ=255)
                              aataatCCATTAACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGAt                                                                      >  1:1958549/1‑100 (MQ=255)
                               ataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGAt                                                                      <  1:3346232/99‑1 (MQ=255)
                                    ccATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGtct                                                                >  2:3486147/1‑100 (MQ=255)
                                    ccATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGtct                                                                >  1:2984741/1‑100 (MQ=255)
                                        ttACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGcc                                                            <  1:2244383/100‑1 (MQ=255)
                                          aCGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTg                                                          >  2:574643/1‑100 (MQ=255)
                                           cGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTAGACATAATAGATTAGTCTCTCCTGc                                                         >  1:660782/1‑100 (MQ=255)
                                            ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGTAAGGCGATGTATCATCGCCGTACATTTTTCGAGATACTAGATTAGTGTTGcc                                                            >  1:3213949/1‑96 (MQ=255)
                                            ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGcc                                                        <  2:2478363/100‑1 (MQ=255)
                                             gAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCa                                                       >  1:3098303/1‑100 (MQ=255)
                                             gAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCa                                                       >  2:2573614/1‑100 (MQ=255)
                                               aTTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGa                                                     <  2:84508/100‑1 (MQ=255)
                                                 ttttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGaaa                                                   <  2:2300165/100‑1 (MQ=255)
                                                   ttGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAtt                                                 <  1:1901443/100‑1 (MQ=255)
                                                      tCTGCAAGATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAttt                                                >  2:1954484/1‑98 (MQ=255)
                                                      tCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAttttt                                              <  2:1018940/100‑1 (MQ=255)
                                                       cTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                             >  1:2052342/1‑100 (MQ=255)
                                                        tGCAAAATACTCCTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGCCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTc                                             <  2:1677751/99‑1 (MQ=255)
                                                        tGCAAAATACTACTGTATTCAGGAGTAAATGGACTGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCt                                            >  1:2461919/1‑100 (MQ=255)
                                                        tGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCt                                            >  2:101436/1‑100 (MQ=255)
                                                        tGCAAAATACTACTGTATTCAGGAGTAAATGGAAGGCGATGTATCATCGACGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAACTTTTTCt                                            >  1:2984923/1‑100 (MQ=255)
                                                         gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                           >  1:2409729/1‑100 (MQ=255)
                                                         gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                           >  2:1053010/1‑100 (MQ=255)
                                                          cAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGc                                          >  2:2141230/1‑100 (MQ=255)
                                                                  tactGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGAc                                  <  2:1767668/100‑1 (MQ=255)
                                                                   actGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATCCTAGATTCGTCTCGCCTGCAAGAAATTTTTCTGCCAATAGACa                                 >  1:570802/1‑100 (MQ=255)
                                                                   actGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGAc                                  >  1:1073869/1‑99 (MQ=255)
                                                                   actGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACa                                 >  2:913127/1‑100 (MQ=255)
                                                                   actGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACa                                 <  1:311685/100‑1 (MQ=255)
                                                                   actGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACa                                 <  1:1056828/100‑1 (MQ=255)
                                                                      gTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTa                              >  2:928193/1‑100 (MQ=255)
                                                                         ttCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTTGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCa                           <  2:2636856/100‑1 (MQ=255)
                                                                           cAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAg                         <  2:707330/100‑1 (MQ=255)
                                                                             tccgatcttGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                       <  2:1210713/92‑1 (MQ=255)
                                                                             ggAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTCTCAAGTg                       >  1:1090647/1‑100 (MQ=255)
                                                                             ggAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                       >  1:3268712/1‑100 (MQ=255)
                                                                             ggAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                       <  1:1538671/100‑1 (MQ=255)
                                                                             ggAGTAAATGGACGGCGATGTATCATCGCCGTCCATATTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                       >  2:2635050/1‑100 (MQ=255)
                                                                              gAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGa                      >  2:1336386/1‑100 (MQ=255)
                                                                              gAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGa                      >  2:2260619/1‑100 (MQ=255)
                                                                              cgaaaaaTGGACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACCTTATCAAGTGa                      >  1:1203772/5‑100 (MQ=37)
                                                                               aGTAAATTGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAg                     <  2:158970/100‑1 (MQ=255)
                                                                                gTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAgg                    <  2:1313562/100‑1 (MQ=255)
                                                                                gTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAgg                    <  2:3336403/100‑1 (MQ=255)
                                                                                  aaaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGaa                  >  1:2874417/1‑100 (MQ=255)
                                                                                  aaaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAATTTTCTGCCAATAGACATTATCAAGTGAGGaa                  >  2:3369624/1‑100 (MQ=255)
                                                                                   aaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAc                 <  1:1378658/100‑1 (MQ=255)
                                                                                     tGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAgatcgga               >  1:1210707/1‑94 (MQ=255)
                                                                                     tGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGc               >  2:1126331/1‑100 (MQ=255)
                                                                                      ngACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa              >  1:1828843/2‑100 (MQ=255)
                                                                                      ggACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa              <  1:1710286/100‑1 (MQ=255)
                                                                                      ggACGGCGATGATACATCGCCGTCCATTTTTTTTCTTTTTATATTCGTCTCGCCTGCCAGAAATTTTTCTCCCTTTAGACATTATCAAGTTACGAACGCa              >  1:3409014/1‑100 (MQ=16)
                                                                                      ggACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa              >  2:105738/1‑100 (MQ=39)
                                                                                       gACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATGAAGTGAGGAACGCAt             >  2:3667181/1‑100 (MQ=255)
                                                                                       gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt             >  1:1449578/1‑100 (MQ=39)
                                                                                       gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt             >  2:905340/1‑100 (MQ=39)
                                                                                       gACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATGTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAt             >  2:2626317/1‑100 (MQ=38)
                                                                                        aCGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATc            <  2:1573513/100‑1 (MQ=255)
                                                                                        aCGGCGATGATACATCGCCGTCCATTTTTCGACATACTGGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGAGAGGAACGCATc            >  2:2182506/1‑100 (MQ=38)
                                                                                         cGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCa           <  1:3680505/100‑1 (MQ=255)
                                                                                            cGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTc        >  2:622742/1‑100 (MQ=255)
                                                                                               tGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCt     >  2:423600/1‑100 (MQ=255)
                                                                                               tGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCt     <  2:2874436/100‑1 (MQ=255)
                                                                                                  atcatcGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGc  <  2:1305746/100‑1 (MQ=255)
                                                                                                   tcatcGTCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGc  <  1:3022780/99‑1 (MQ=255)
                                                                                                    |                                                                                                  
CCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGC  >  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: