Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 850,995 T→G 15.1% intergenic (+30/+19) 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,9950TG15.1% 356.4 / 54.4 149intergenic (+30/+19)ompX/opgEouter membrane protein X/OPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base T (66/59);  new base G (8/15);  total (74/75)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

Only 100 of 154 total aligned reads displayed.

TACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGG  >  NC_000913/850897‑851094
                                                                                                  |                                                                                                   
tACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCtctc                                                                                                    <  1:3027105/100‑1 (MQ=255)
tACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCtctc                                                                                                    >  2:2274320/1‑100 (MQ=255)
    aGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCgggg                                                                                                >  1:3656490/1‑100 (MQ=255)
    aGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCgggg                                                                                                >  1:3704776/1‑100 (MQ=255)
          gCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAtt                                                                                          >  2:3956455/1‑100 (MQ=255)
           ccGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttn                                                                                         <  2:2134793/100‑2 (MQ=255)
           ccGTAGTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCGAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttt                                                                                         <  2:599512/100‑1 (MQ=38)
            cGTATTCGTAGCGTTGACGTAGGCAGCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAGATCCGCCCCGAGAGGCGGAtt                                                                                          <  2:3358245/98‑1 (MQ=255)
            cGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAtttt                                                                                        >  1:545934/1‑100 (MQ=255)
             gTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAttttt                                                                                       <  1:1726567/100‑1 (MQ=255)
              tATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttttt                                                                                      <  2:3417820/100‑2 (MQ=38)
              tATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttttt                                                                                      <  1:2702947/100‑2 (MQ=38)
                ttCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttttttat                                                                                    <  2:3652017/100‑4 (MQ=37)
                ttCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttttttat                                                                                    <  1:3962964/100‑4 (MQ=37)
                 tCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGTGCGGATTTTTGttt                                                                                   <  2:1477072/100‑1 (MQ=255)
                 tCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGttt                                                                                   <  2:4099901/100‑1 (MQ=255)
                 tCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttttttata                                                                                   <  1:1548061/100‑5 (MQ=37)
                  cGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGtttt                                                                                  <  1:2638621/100‑1 (MQ=255)
                  cGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGtttt                                                                                  <  1:74203/100‑1 (MQ=255)
                  cGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGtttt                                                                                  <  1:3974861/100‑1 (MQ=255)
                  cGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGttt                                                                                   <  2:1093359/99‑1 (MQ=255)
                  cGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTGCTAATCACTTTGTTTATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGtttt                                                                                  <  2:2054886/100‑1 (MQ=255)
                   gTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttttttatatc                                                                                 <  1:916260/100‑7 (MQ=37)
                    tAGCTTTGACGTAGGCACCTGGATTTCCGATATTGGTTACCGCGTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttttttatat                                                                                  <  2:240465/98‑6 (MQ=21)
                     aGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttttttatatcac                                                                               <  1:1318247/100‑9 (MQ=25)
                      gCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAg                                                                              <  2:1870751/100‑1 (MQ=255)
                      gCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTGGGTGATATAAAAAATACGCCTCTCGGGGCGGATTTTTGTTTTGaa                                                                               >  2:3065790/1‑99 (MQ=255)
                      gCGTTGACGTAGGCACCTGGATGGCCGGTGTTGGTTACCGCTGCTAATCGCTTTGGTGATATAAAAAAACCGCCGCTCGGGGCGGATTTTTGTTTTTa                                                                                <  2:2680013/98‑1 (MQ=255)
                       cGTTGACGTCGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTACTCACTTTGCCGCTATCAACAATCAACCTCTCGGGCCAGCtttttcatc                                                                                   >  2:3800347/1‑90 (MQ=255)
                       cGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAgg                                                                             >  2:3307766/1‑100 (MQ=255)
                       cGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAgg                                                                             <  1:93241/100‑1 (MQ=255)
                       cGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAgg                                                                             <  1:1082493/100‑1 (MQ=255)
                           gACGTAGGCACCTGGATTGCCGGTGCTGGTTACCGCTTCTAATCACTTTGGTGATATACAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGAttt                                                                          >  2:1810668/1‑99 (MQ=255)
                              gTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAGAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCggg                                                                      >  1:487433/1‑100 (MQ=255)
                              gTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCggg                                                                      >  1:2831929/1‑100 (MQ=255)
                               tAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGt                                                                     >  1:152615/1‑100 (MQ=255)
                                     ccTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGaaaa                                                               >  2:3195539/1‑100 (MQ=255)
                                        ggATTGCCGGTGTTGGATACCGCTTCTAATAACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAatat                                                            >  2:19384/1‑100 (MQ=255)
                                         gATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATc                                                           <  2:3142459/100‑1 (MQ=255)
                                         gATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATc                                                           <  2:389199/100‑1 (MQ=255)
                                         gATTGCCGGTGTTGGTTACCGCTTCTAATCACTGTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTGTCGGGTCGAAAATATc                                                           >  1:186675/1‑100 (MQ=255)
                                          aTTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCg                                                          <  1:2742328/100‑1 (MQ=255)
                                          aTTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCg                                                          <  2:163229/100‑1 (MQ=255)
                                             gCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGGTTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGtt                                                       >  1:2039848/1‑100 (MQ=255)
                                             gCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGtt                                                       >  1:2294554/1‑100 (MQ=255)
                                             gCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGtt                                                       >  2:3791099/1‑100 (MQ=255)
                                                ggTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCg                                                    <  2:1476609/100‑1 (MQ=255)
                                                  tgtTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCgag                                                  <  1:2885958/100‑1 (MQ=255)
                                                    ttGGTTACCGCTTCTAATCACTTTGGTGTTATAAAAAATCCGCCTCTCGGGGCGGATTGTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGAt                                                <  1:221504/100‑1 (MQ=255)
                                                     tGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCGGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATg                                               >  2:1791842/1‑100 (MQ=255)
                                                     tGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATg                                               <  1:2261745/100‑1 (MQ=255)
                                                     tGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATg                                               >  1:3213604/1‑100 (MQ=255)
                                                      ggTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATgg                                              >  1:1724966/1‑100 (MQ=255)
                                                      ggTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATgg                                              >  1:2472202/1‑100 (MQ=255)
                                                      agtaactgCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGATAATATCGGtt                                                       >  2:2430465/8‑91 (MQ=255)
                                                        ttACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGtt                                            >  2:1646969/1‑100 (MQ=255)
                                                         tACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAgttggttg                                           >  2:1329827/1‑100 (MQ=255)
                                                         tACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGtt                                            >  2:3869203/1‑99 (MQ=255)
                                                         tACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTg                                           <  1:156491/100‑1 (MQ=255)
                                                            cGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTaa                                        >  1:632964/1‑100 (MQ=255)
                                                             gCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGACAATATCGGTTCCGAGATGGTTGTaa                                        >  2:3461743/1‑99 (MQ=255)
                                                               ttCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATcc                                     <  1:2978513/100‑1 (MQ=255)
                                                               ttCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATcc                                     >  2:1292492/1‑100 (MQ=255)
                                                               ttCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATcc                                     <  1:1486129/100‑1 (MQ=255)
                                                                tCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCa                                    <  2:1900081/100‑1 (MQ=255)
                                                                  tAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCAcc                                  >  2:3116084/1‑100 (MQ=255)
                                                                  tAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCAcc                                  >  1:753443/1‑100 (MQ=255)
                                                                     tCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTc                               <  2:3663544/100‑1 (MQ=255)
                                                                       aCTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTg                             <  2:283943/100‑1 (MQ=255)
                                                                        cTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGt                            >  1:187291/1‑100 (MQ=255)
                                                                        cTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGt                            >  1:4137303/1‑100 (MQ=255)
                                                                         tttGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTa                           <  1:1074101/100‑1 (MQ=255)
                                                                          ttGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTaa                          >  2:102008/1‑100 (MQ=255)
                                                                          ttGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAAGTTTCGGGTCGAAAATATCGGCTCCGAGATGGTTGTAATCCACCTTCTGTa                           >  2:2105263/1‑99 (MQ=255)
                                                                           tGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAAc                         >  1:2645188/1‑100 (MQ=255)
                                                                            cttaaaaacaaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACt                        >  1:1172907/10‑100 (MQ=255)
                                                                            cttaaaaacaaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAACATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACt                        >  2:1934167/10‑100 (MQ=25)
                                                                                  ataAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATAGGGTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGt                  >  2:553617/1‑100 (MQ=255)
                                                                                   taAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGtt                 >  1:3801289/1‑100 (MQ=255)
                                                                                   acaaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGACAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGtt                 >  2:2747549/3‑100 (MQ=38)
                                                                                   acaaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGtt                 >  1:3399683/3‑100 (MQ=38)
                                                                                    caaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGCTGGTTGTAATCCACCTTCTGTAACTGGAAGTTg                >  1:3448923/2‑100 (MQ=38)
                                                                                     aaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTgg               <  1:46221/100‑1 (MQ=255)
                                                                                     aaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTgg               <  1:3081094/100‑1 (MQ=255)
                                                                                     aaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTgg               >  1:2948069/1‑100 (MQ=255)
                                                                                     aaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTgg               <  1:2299712/100‑1 (MQ=255)
                                                                                     aaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTgg               >  1:528822/1‑100 (MQ=38)
                                                                                     aaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTgg               >  1:4097801/1‑100 (MQ=38)
                                                                                      aaaaTCCTCCTCTCGGGGTGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGt              <  2:2834879/100‑1 (MQ=255)
                                                                                      aaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGt              >  1:3105637/1‑100 (MQ=255)
                                                                                       aaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTg             >  1:3426592/1‑100 (MQ=255)
                                                                                         aTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGAt           >  1:2362408/1‑100 (MQ=255)
                                                                                          tCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATg          <  1:3568139/100‑1 (MQ=255)
                                                                                            cGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGta        <  2:964088/100‑1 (MQ=255)
                                                                                            cGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGta        <  2:3085136/100‑1 (MQ=255)
                                                                                             gCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGtat       >  2:273507/1‑100 (MQ=255)
                                                                                             gCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAAACCACCTTCTGTAACTGGAAGTTGGTGATGta        <  2:1784025/99‑1 (MQ=255)
                                                                                                tctcGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATc    <  1:3157946/100‑1 (MQ=255)
                                                                                                  tcGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCgg  <  1:3111049/100‑1 (MQ=255)
                                                                                                  tcGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCgg  <  1:2970933/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGG  >  NC_000913/850897‑851094

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