Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,792,240 Δ1 bp 15.6% intergenic (‑220/+27) 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,2400T.15.6% 352.7 / 67.3 128intergenic (‑220/+27)ydiV/nlpCanti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (57/45);  new base . (9/11);  total (69/59)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.34e-01

Only 100 of 128 total aligned reads displayed.

CCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCT  >  NC_000913/1792143‑1792337
                                                                                                 |                                                                                                  
ccTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat                                                                                                  <  1:1517138/100‑1 (MQ=255)
ccTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat                                                                                                  >  2:3255438/1‑100 (MQ=255)
 cTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat‑c                                                                                                >  2:1669047/1‑100 (MQ=255)
  tttttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGCCGATGTat‑ca                                                                                               >  1:2445762/1‑100 (MQ=255)
    tttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat‑catc                                                                                             <  2:2868191/100‑1 (MQ=255)
        tAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGGCTGCAAAATACTACTGTATTCAGGATTAAATGTACGTCGACGTAT‑CATCGCCg                                                                                         <  1:1310971/100‑1 (MQ=255)
         aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGt                                                                                        >  2:1889634/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                                                                                       <  1:890987/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTc                                                                                       <  2:1218644/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTc                                                                                       <  2:587097/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTc                                                                                       <  2:844668/100‑1 (MQ=39)
            ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTCCa                                                                                     <  2:2694168/100‑1 (MQ=38)
            ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATG‑ATACATCGCCGTCCa                                                                                     <  1:3073009/100‑1 (MQ=38)
                ccTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCAtttt                                                                                 >  2:3680533/1‑100 (MQ=255)
                    aTTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGa                                                                             <  1:350293/100‑1 (MQ=255)
                     ttATAAACGAATAATTCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGAc                                                                            <  2:3672959/100‑1 (MQ=255)
                     ttATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGAc                                                                            <  2:1615453/100‑1 (MQ=255)
                     ttATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGAc                                                                            <  2:391720/100‑1 (MQ=255)
                     ttATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGAc                                                                            >  2:2152024/1‑100 (MQ=255)
                      tataAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACa                                                                           >  1:555778/1‑100 (MQ=255)
                      tataAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTAGACa                                                                           >  1:2429485/1‑100 (MQ=255)
                         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:2715760/100‑1 (MQ=255)
                              aataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGAt                                                                   <  1:3585792/100‑1 (MQ=255)
                               ataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGAt                                                                   <  1:3346232/99‑1 (MQ=255)
                                  atCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGt                                                               >  2:1056834/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‑CATCGCGGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGc                                                      >  2:2619231/1‑100 (MQ=255)
                                           cGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTAGACATAATAGATTAGTCTCTCCTGc                                                      >  1:660782/1‑100 (MQ=255)
                                            ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGTAAGGCGATGTAT‑CATCGCCGTACATTTTTCGAGATACTAGATTAGTGTTGcc                                                         >  1:3213949/1‑96 (MQ=255)
                                             gAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCa                                                    >  1:3098303/1‑100 (MQ=255)
                                             gAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCa                                                    >  1:1334147/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)
                                                       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                                        >  2:1053010/1‑100 (MQ=255)
                                                         gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                        >  1:2409729/1‑100 (MQ=255)
                                                         gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                        >  1:2968577/1‑100 (MQ=255)
                                                          cAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGc                                       <  2:3288957/100‑1 (MQ=255)
                                                          cAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGc                                       >  2:2141230/1‑100 (MQ=255)
                                                           taaaTACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGcc                                      <  2:564924/99‑1 (MQ=255)
                                                              aTACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAAt                                   >  2:2389787/1‑100 (MQ=255)
                                                                  tactGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGAc                               <  2:1767668/100‑1 (MQ=255)
                                                                   actGTATTCATTAACAAATGGACGGCGATGTAT‑CATCGCTGTCCCTTTTTCGACATACTAGGTTCGTGCGGCCTGCCGGAAATTATTCTGCAGATAGACa                              >  1:2062370/1‑100 (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)
                                                                         ttCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCa                        <  1:3255416/100‑1 (MQ=255)
                                                                           ccgaaGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTACATTCGTCTCGTCTGCCAGAAATTTTTCTGACTATTGACATTATCaa                       <  1:690839/95‑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)
                                                                              gAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGa                   >  2:2260619/1‑100 (MQ=255)
                                                                              gAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGa                   >  2:1336386/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)
                                                                                 aaaaTGGACGGCGATG‑ATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGa                >  2:35450/2‑100 (MQ=38)
                                                                                  aaaTGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGaa               >  1:2802096/1‑100 (MQ=255)
                                                                                  aaaTGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGaa               >  1:2874417/1‑100 (MQ=255)
                                                                                  aaaTGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAATTTTCTGCCAATAGACATTATCAAGTGAGGaa               >  2:3369624/1‑100 (MQ=255)
                                                                                   aaTGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATTGACATTATCAAGTGAGGAAc              >  1:995515/1‑100 (MQ=255)
                                                                                   aaTGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAc              <  1:1378658/100‑1 (MQ=255)
                                                                                     tGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAgatcgga            >  1:1210707/1‑94 (MQ=255)
                                                                                      ngACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa           >  1:1828843/2‑100 (MQ=255)
                                                                                      ggACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa           <  1:1710286/100‑1 (MQ=255)
                                                                                      ggACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa           >  2:687025/1‑100 (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:3001484/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)
                                                                                         cGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCa        <  1:3680505/100‑1 (MQ=255)
                                                                                               tGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCt  <  2:2874436/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
CCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAT‑CATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCT  >  NC_000913/1792143‑1792337

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