Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 850,996 C→A 16.9% intergenic (+31/+18) ompX → / ← opgE outer membrane protein X/OPG biosynthetic transmembrane phosphoethanolamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913850,9960CA16.9% 404.0 / 64.1 152intergenic (+31/+18)ompX/opgEouter membrane protein X/OPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base C (67/59);  new base A (9/17);  total (76/76)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.91e-01

Only 100 of 160 total aligned reads displayed.

TACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCG  >  NC_000913/850897‑851093
                                                                                                   |                                                                                                 
tACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCtctc                                                                                                   >  2:1664312/1‑100 (MQ=255)
tACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCtctc                                                                                                   >  1:926629/1‑100 (MQ=255)
    aGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCgggg                                                                                               >  1:204603/1‑100 (MQ=255)
      cAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCg                                                                                             <  2:2247556/100‑1 (MQ=255)
          gCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAtt                                                                                         <  2:2981124/100‑1 (MQ=255)
            catATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCTCTTTGGTGATATAAAAAGACCGCCCCTAGAGGCGGAtttt                                                                                       <  1:1193730/98‑1 (MQ=255)
            cGTATTCGTAGTGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttt                                                                                       <  1:478691/100‑1 (MQ=38)
            cGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttt                                                                                       <  2:3597255/100‑1 (MQ=38)
             gTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttt                                                                                       <  2:1710451/99‑1 (MQ=38)
              tATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTg                                                                                     <  2:2183950/100‑1 (MQ=255)
              tATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTg                                                                                     <  2:1319071/100‑1 (MQ=255)
              tATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTg                                                                                     <  2:2176079/100‑1 (MQ=255)
              tATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttttt                                                                                     <  1:1917009/100‑2 (MQ=38)
              tATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttttt                                                                                     <  1:3165290/100‑2 (MQ=38)
              tATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttttt                                                                                     <  1:2636931/100‑2 (MQ=38)
              tATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttttt                                                                                     <  2:1238639/100‑2 (MQ=38)
               aTTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttttta                                                                                    <  2:3552258/100‑3 (MQ=38)
                ttCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGtt                                                                                   >  2:2389786/1‑100 (MQ=255)
                ttCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttttta                                                                                    <  1:826632/99‑3 (MQ=38)
                 tCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttttttata                                                                                  <  1:1986454/100‑5 (MQ=37)
                 tCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttttttata                                                                                  <  2:52879/100‑5 (MQ=37)
                   gTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGttttt                                                                                >  1:2788273/1‑100 (MQ=255)
                   gTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGttttt                                                                                >  1:369443/1‑100 (MQ=255)
                     aGCGTTGACGTAGGCACCTGGATTGCCTGTGTTGGTTGCCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttttttatatc                                                                                <  1:1943063/98‑7 (MQ=25)
                      gCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAg                                                                             >  2:739387/1‑100 (MQ=255)
                      gCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttttttatatcacc                                                                             <  2:2381265/100‑10 (MQ=25)
                       cGTTGGCGGAGGCACCTGGATTGCCGTTGGTGGTTGCCTCTTCTATTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAg                                                                             <  2:159186/99‑1 (MQ=255)
                            aCGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCg                                                                       <  2:1127161/100‑1 (MQ=255)
                             cGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCgg                                                                      >  2:538379/1‑100 (MQ=255)
                                aGGCACCTGGATTGCCGGGGTTGGTTACCGCTTCTAATCGCATTGGTAATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTGTAAgat                                                                           >  1:3154/1‑90 (MQ=255)
                                 ggCACCTGGGTTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCg                                                                  <  1:1000607/100‑1 (MQ=255)
                                 ggCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCg                                                                  >  1:3364660/1‑100 (MQ=255)
                                 ggCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCg                                                                  <  2:1736161/100‑1 (MQ=255)
                                  gCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGa                                                                 >  2:3520300/1‑100 (MQ=255)
                                  gCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTAGGTGATATAAAAAATCCGCCTCTCGGGGCGGCTTTTTGTTTTTAAGGTTTCGGGTCGa                                                                 >  2:133552/1‑100 (MQ=255)
                                   cACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCGGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGaa                                                                >  1:1442149/1‑100 (MQ=255)
                                    aCCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGACATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGaaa                                                               <  2:662515/100‑1 (MQ=255)
                                      cTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAat                                                             <  2:3274134/100‑1 (MQ=255)
                                        ggATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAata                                                            <  1:1189972/99‑1 (MQ=255)
                                         gATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTTATATAAAAAATCCGCCTCTCGGAGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATc                                                          <  1:960822/100‑1 (MQ=255)
                                          aTTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCg                                                         >  2:3815230/1‑100 (MQ=255)
                                           ttGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCgg                                                        <  2:1574395/100‑1 (MQ=255)
                                           ttGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCgg                                                        <  2:441665/100‑1 (MQ=255)
                                             gCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGtt                                                      >  2:3648105/1‑100 (MQ=255)
                                                ggTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCg                                                   <  1:1725941/100‑1 (MQ=255)
                                                ggTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCg                                                   <  1:2153449/100‑1 (MQ=255)
                                                 gtgtTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCga                                                  <  2:3325697/100‑1 (MQ=255)
                                                   gtTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCATCTCGGGGCGGATTTTTGTTTTTGAGGTTTCGGGTGGCAAATATCGGTTc                                                     >  1:3759823/1‑95 (MQ=255)
                                                   gtTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCAGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCgaga                                                >  1:727682/1‑100 (MQ=255)
                                                    ttGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGAt                                               <  1:3496242/100‑1 (MQ=255)
                                                    ttGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCAGAGAt                                               >  1:360613/1‑100 (MQ=255)
                                                       gTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGt                                            <  1:3499293/100‑1 (MQ=255)
                                                        ttACCGCTTCTAATCATTTTGTTGATATAAAAAATCCGACTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGtt                                           <  1:3651950/100‑1 (MQ=255)
                                                         tACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTg                                          <  2:3261478/100‑1 (MQ=255)
                                                           ccGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTa                                        >  1:1749342/1‑100 (MQ=255)
                                                            cGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAACATATCGGTTCCGAGATGGTTGTaa                                       >  1:553400/1‑100 (MQ=255)
                                                             gCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAAt                                      <  2:391016/100‑1 (MQ=255)
                                                                 cTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCAc                                  >  2:234322/1‑100 (MQ=255)
                                                                  tAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCAcc                                 >  2:140602/1‑100 (MQ=255)
                                                                    aTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCtt                               >  2:2448677/1‑100 (MQ=255)
                                                                     tCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTc                              >  2:1163468/1‑100 (MQ=255)
                                                                      cACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCt                             >  1:67859/1‑100 (MQ=255)
                                                                        cTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTTAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGt                           <  1:515570/100‑1 (MQ=255)
                                                                        cTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGt                           >  1:280545/1‑100 (MQ=255)
                                                                         tttGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTa                          <  2:3013380/100‑1 (MQ=255)
                                                                           tGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAAc                        >  2:2042047/1‑100 (MQ=255)
                                                                           tGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAAc                        <  1:1620217/100‑1 (MQ=255)
                                                                             gTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTg                      <  2:1669719/100‑1 (MQ=255)
                                                                             gTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTg                      <  1:327985/100‑1 (MQ=255)
                                                                             gTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTg                      <  1:2136872/100‑1 (MQ=255)
                                                                              tGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAGGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTgg                     <  1:3781259/100‑1 (MQ=255)
                                                                              tGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTgg                     >  2:1532715/1‑100 (MQ=255)
                                                                               gATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGa                    <  1:2469754/100‑1 (MQ=255)
                                                                                aaaacaaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGaa                   >  2:128125/6‑100 (MQ=37)
                                                                                aaaacaaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGaa                   >  1:3395961/6‑100 (MQ=37)
                                                                                 tataAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAg                  <  2:3207665/100‑1 (MQ=255)
                                                                                 tataAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAg                  >  1:2170325/1‑100 (MQ=255)
                                                                                  ataAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGt                 >  2:2968345/1‑100 (MQ=255)
                                                                                  aacaaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGt                 >  1:2882733/4‑100 (MQ=37)
                                                                                   taAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGt                 >  1:1041679/1‑99 (MQ=255)
                                                                                   acaaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTACTCCACCTTCTGTAACTGGAAGtt                >  1:3041362/3‑100 (MQ=37)
                                                                                   acaaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGtt                >  2:1826005/3‑100 (MQ=38)
                                                                                    caaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTg               >  2:3025277/2‑100 (MQ=38)
                                                                                    aaaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTgg                     >  1:1572220/1‑94 (MQ=255)
                                                                                    aaaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTg               >  2:671827/1‑100 (MQ=255)
                                                                                    aaaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCAGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAg                  >  1:3028457/1‑97 (MQ=255)
                                                                                    aaaaaaTCCGCCTCTCGGGGCGGATTTTTGTTATTAAGGTTTCGGGTCGAAAATATCGCTTCCTAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTg               >  2:2672225/1‑100 (MQ=255)
                                                                                     aaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTgg              <  2:1181510/100‑1 (MQ=255)
                                                                                     aaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTgg              <  2:1058382/100‑1 (MQ=255)
                                                                                      aaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGt             <  2:3732336/100‑1 (MQ=255)
                                                                                       aaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTg            <  2:3630718/100‑1 (MQ=255)
                                                                                           ccGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGt        <  1:1290555/100‑1 (MQ=255)
                                                                                            cGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGTTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGta       <  2:3069356/100‑1 (MQ=255)
                                                                                            cGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGta       <  1:1214634/100‑1 (MQ=255)
                                                                                            cGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGta       >  1:2166052/1‑100 (MQ=255)
                                                                                             gCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGtat      >  2:1646304/1‑100 (MQ=255)
                                                                                              ccTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGtata     <  2:3779550/100‑1 (MQ=255)
                                                                                               ctctCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTatat    <  2:198753/100‑1 (MQ=255)
                                                                                                tctcGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTatat    >  1:1386611/1‑99 (MQ=255)
                                                                                                 ctcGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCg  >  2:3234697/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
TACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCG  >  NC_000913/850897‑851093

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