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.

CCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCT  >  NC_000913/1792143‑1792337
                                                                                                 |                                                                                                  
ccTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat                                                                                                  <  1:1517138/100‑1 (MQ=255)
ccTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTat                                                                                                  >  2:3255438/1‑100 (MQ=255)
 cTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTatc                                                                                                >  2:1669047/1‑100 (MQ=255)
  tttttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGCCGATGTatca                                                                                               >  1:2445762/1‑100 (MQ=255)
    tttttAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTatcatc                                                                                             <  2:2868191/100‑1 (MQ=255)
        tAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGGCTGCAAAATACTACTGTATTCAGGATTAAATGTACGTCGACGTATCATCGCCg                                                                                         <  1:1310971/100‑1 (MQ=255)
         aGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGt                                                                                        >  2:1889634/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                                                                                       <  1:890987/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                       <  2:1218644/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                       <  2:587097/100‑1 (MQ=39)
          gATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                       <  2:844668/100‑1 (MQ=39)
            ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTCCa                                                                                     <  2:2694168/100‑1 (MQ=38)
            ttATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTCCa                                                                                     <  1:3073009/100‑1 (MQ=38)
                ccTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCAtttt                                                                                 >  2:3680533/1‑100 (MQ=255)
                    aTTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGa                                                                             <  1:350293/100‑1 (MQ=255)
                     ttATAAACGAATAATTCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGAc                                                                            <  2:3672959/100‑1 (MQ=255)
                     ttATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGAc                                                                            <  2:1615453/100‑1 (MQ=255)
                     ttATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGAc                                                                            <  2:391720/100‑1 (MQ=255)
                     ttATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGAc                                                                            >  2:2152024/1‑100 (MQ=255)
                      tataAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACa                                                                           >  1:555778/1‑100 (MQ=255)
                      tataAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTAGACa                                                                           >  1:2429485/1‑100 (MQ=255)
                         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:2715760/100‑1 (MQ=255)
                              aataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGAt                                                                   <  1:3585792/100‑1 (MQ=255)
                               ataatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGAt                                                                   <  1:3346232/99‑1 (MQ=255)
                                  atCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGt                                                               >  2:1056834/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)
                                           cGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCGGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGc                                                      >  2:2619231/1‑100 (MQ=255)
                                           cGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTAGACATAATAGATTAGTCTCTCCTGc                                                      >  1:660782/1‑100 (MQ=255)
                                            ggAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGTAAGGCGATGTATCATCGCCGTACATTTTTCGAGATACTAGATTAGTGTTGcc                                                         >  1:3213949/1‑96 (MQ=255)
                                             gAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCa                                                    >  1:3098303/1‑100 (MQ=255)
                                             gAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCa                                                    >  1:1334147/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)
                                                       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                                        >  2:1053010/1‑100 (MQ=255)
                                                         gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                        >  1:2409729/1‑100 (MQ=255)
                                                         gCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTg                                        >  1:2968577/1‑100 (MQ=255)
                                                          cAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGc                                       <  2:3288957/100‑1 (MQ=255)
                                                          cAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGc                                       >  2:2141230/1‑100 (MQ=255)
                                                           taaaTACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGcc                                      <  2:564924/99‑1 (MQ=255)
                                                              aTACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAAt                                   >  2:2389787/1‑100 (MQ=255)
                                                                  tactGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGAc                               <  2:1767668/100‑1 (MQ=255)
                                                                   actGTATTCATTAACAAATGGACGGCGATGTATCATCGCTGTCCCTTTTTCGACATACTAGGTTCGTGCGGCCTGCCGGAAATTATTCTGCAGATAGACa                              >  1:2062370/1‑100 (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)
                                                                         ttCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCa                        <  1:3255416/100‑1 (MQ=255)
                                                                           ccgaaGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTACATTCGTCTCGTCTGCCAGAAATTTTTCTGACTATTGACATTATCaa                       <  1:690839/95‑1 (MQ=255)
                                                                           cAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAg                      <  2:707330/100‑1 (MQ=255)
                                                                             tccgatcttGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTg                    <  2:1210713/92‑1 (MQ=255)
                                                                             ggAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTCTCAAGTg                    >  1:1090647/1‑100 (MQ=255)
                                                                              gAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGa                   >  2:2260619/1‑100 (MQ=255)
                                                                              gAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGa                   >  2:1336386/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)
                                                                                 aaaaTGGACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGa                >  2:35450/2‑100 (MQ=38)
                                                                                  aaaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGaa               >  1:2802096/1‑100 (MQ=255)
                                                                                  aaaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGaa               >  1:2874417/1‑100 (MQ=255)
                                                                                  aaaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAAATTTTCTGCCAATAGACATTATCAAGTGAGGaa               >  2:3369624/1‑100 (MQ=255)
                                                                                   aaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATTGACATTATCAAGTGAGGAAc              >  1:995515/1‑100 (MQ=255)
                                                                                   aaTGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAAc              <  1:1378658/100‑1 (MQ=255)
                                                                                     tGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAgatcgga            >  1:1210707/1‑94 (MQ=255)
                                                                                      ngACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa           >  1:1828843/2‑100 (MQ=255)
                                                                                      ggACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa           <  1:1710286/100‑1 (MQ=255)
                                                                                      ggACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCa           >  2:687025/1‑100 (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:3001484/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)
                                                                                         cGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCa        <  1:3680505/100‑1 (MQ=255)
                                                                                               tGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCt  <  2:2874436/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
CCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCT  >  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: