Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 850,995 T→G 6.3% 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,9950TG 6.3% 266.5 / 8.6 95intergenic (+30/+19)ompX/opgEouter membrane protein X/OPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base T (41/47);  new base G (3/3);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

GTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCT  >  NC_000913/850847‑851141
                                                                                                                                                    |                                                                                                                                                  
gTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCctct                                                                                                                                                    <  3:265385/149‑1 (MQ=255)
    gggTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCggg                                                                                                                                                <  5:476532/149‑1 (MQ=255)
    gggTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCggg                                                                                                                                                <  6:130185/149‑1 (MQ=255)
    gggTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCACggg                                                                                                                                                >  5:354340/1‑149 (MQ=255)
      gTCTGCAGTTCAACGCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGc                                                                                                                                              <  8:15667/149‑1 (MQ=255)
      gTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGc                                                                                                                                              >  2:85742/1‑149 (MQ=255)
            aGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAttt                                                                                                                                        <  4:170988/149‑1 (MQ=255)
              ttaaaCGCGATGGAAAACGTTGCTCTGGACTTCTCTTTCGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCAGGATTGCCGGTGTTGGTTACCTCTTCTAATCTCTTTGGTGTTATAAAAAATCCGCCCCGAGAGGCGGAttttt                                                                                                                                      <  4:239270/146‑1 (MQ=255)
              ttCAACCCGATGGAAATCGTTGCTCTGGAGTTCTCTTACGAGCAGAGCCGTATTCGTTGCGTTGACGTAGGCACCTGTATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAttttt                                                                                                                                      <  8:154530/149‑1 (MQ=255)
              ttCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAttttt                                                                                                                                      <  2:334912/149‑1 (MQ=255)
              ttCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAttttt                                                                                                                                      <  6:364997/149‑1 (MQ=255)
              ttCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTATGGTGATATAAAAAATCCGCCTCTCGGGGCGGAttttt                                                                                                                                      <  5:421631/149‑1 (MQ=255)
              ttCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGCATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAttttt                                                                                                                                      <  8:419342/149‑1 (MQ=255)
                      gATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGttttt                                                                                                                                >  1:113723/1‑147 (MQ=255)
                      gATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTaa                                                                                                                              <  2:113723/149‑1 (MQ=255)
                          gAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAACCACTTTGGTGATATAAAAAATACGCCTCTCGGGGCGGATTTTTTTTTTTAAGGtt                                                                                                                          >  8:507279/1‑149 (MQ=255)
                           aaaaCGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTaa                                                                                                                              >  3:440765/1‑144 (MQ=255)
                           aaaaCGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTaa                                                                                                                              <  4:440765/144‑1 (MQ=255)
                           aaaaCGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGttt                                                                                                                         >  1:364414/1‑149 (MQ=255)
                            aaaCGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTc                                                                                                                        >  8:72947/1‑149 (MQ=255)
                                     tctGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGaa                                                                                                                >  1:227280/1‑148 (MQ=255)
                                       tGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGaaaa                                                                                                              >  4:399413/1‑148 (MQ=255)
                                         gACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAatat                                                                                                           <  8:406374/149‑1 (MQ=255)
                                           cTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCg                                                                                                         <  1:283753/149‑1 (MQ=255)
                                              ctctTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttttt                                                                                                                                     >  7:251903/1‑117 (MQ=39)
                                              ctctTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttttt                                                                                                                                     <  8:251903/118‑2 (MQ=39)
                                                 ttACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTcc                                                                                                    >  8:126159/1‑148 (MQ=255)
                                                 ttACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTcc                                                                                                    <  1:36496/148‑1 (MQ=255)
                                                  tACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttttt                                                                                                                                      <  7:105327/113‑1 (MQ=39)
                                                  tACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttttt                                                                                                                                      >  8:105327/1‑113 (MQ=39)
                                                   aCGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAttttt                                                                                                                                      >  7:6443/1‑112 (MQ=255)
                                                   aCGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAttttt                                                                                                                                      <  8:6443/112‑1 (MQ=255)
                                                     gAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAat                                                                                                             >  7:80582/1‑135 (MQ=255)
                                                     gAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAat                                                                                                             <  8:80582/135‑1 (MQ=255)
                                                       gCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATgg                                                                                             <  7:71052/149‑1 (MQ=255)
                                                              cGTATTCGTAGCGTTGACGTAGGCACCTGGAGTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAAt                                                                                      >  5:460322/1‑149 (MQ=255)
                                                                tATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGGAACCGCTTCTAATCTCTTTGTTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATcc                                                                                    <  6:321508/149‑1 (MQ=255)
                                                                    cGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCAcc                                                                                 <  7:126159/148‑1 (MQ=255)
                                                                     gTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCtt                                                                               >  3:215471/1‑149 (MQ=255)
                                                                     gTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCtt                                                                               >  6:115025/1‑149 (MQ=255)
                                                                     gTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCtt                                                                               >  6:113377/1‑149 (MQ=255)
                                                                     gTAGCGTTGACGTAGGCACCTGGATTGCCGCTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAATGTTTCGGGTAGAAAATATCGGTTCCGAGATGGTTGTAATCCTCCtt                                                                               >  3:100435/1‑149 (MQ=255)
                                                                      tAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTc                                                                              >  7:313641/1‑149 (MQ=255)
                                                                        gCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCt                                                                             <  2:364414/148‑1 (MQ=255)
                                                                         cgatGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATAGGTTCCGAGATGGTTGTAATCCACCTTCTg                                                                            >  8:467588/4‑148 (MQ=255)
                                                                         cGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATTGGTTCCGAGATGGTTGTAATCCACCTTCTg                                                                            <  2:227280/148‑1 (MQ=255)
                                                                         cGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGt                                                                           <  7:221976/149‑1 (MQ=255)
                                                                              aCGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCCCCTTCTGTAACTg                                                                      >  5:153819/1‑149 (MQ=255)
                                                                                 tAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTgg                                                                     >  4:145364/1‑147 (MQ=255)
                                                                                 tAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGaa                                                                   <  5:115025/149‑1 (MQ=255)
                                                                                 tAGGCACCTGGATTGACGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTgg                                                                     >  8:148746/1‑147 (MQ=255)
                                                                                  aGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAg                                                                  >  7:198981/1‑149 (MQ=255)
                                                                                    gCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGTGTCGAAAATATCGGTTCCGAGATGGTTTTAATCCACCTTCTGTAACTGGAAGtt                                                                >  1:462442/1‑149 (MQ=255)
                                                                                      aCCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTgg                                                              >  2:172446/1‑149 (MQ=255)
                                                                                      aCCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTg                                                               >  2:409905/1‑148 (MQ=255)
                                                                                       ccTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGt                                                             <  4:215471/149‑1 (MQ=255)
                                                                                        cTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTg                                                               <  1:409905/146‑1 (MQ=255)
                                                                                           gATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGa                                                           <  8:313641/147‑1 (MQ=255)
                                                                                           gATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGAt                                                          <  5:144081/148‑1 (MQ=255)
                                                                                            aTTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGt                                                        <  3:399413/149‑1 (MQ=255)
                                                                                            aTTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGt                                                        <  1:85742/149‑1 (MQ=255)
                                                                                                ccGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGAt                                                          >  4:135545/1‑143 (MQ=255)
                                                                                                ccGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGAt                                                          <  3:135545/143‑1 (MQ=255)
                                                                                                 cGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATc                                                   >  3:219817/1‑149 (MQ=255)
                                                                                                      ttGGTTAGCGCTTCTAACCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGcc                                               <  7:507279/148‑1 (MQ=255)
                                                                                                          tgaCCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCt                                          <  1:365174/147‑1 (MQ=255)
                                                                                                            aCCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATc                                                                                                          <  2:241693/83‑1 (MQ=255)
                                                                                                            aCCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCtt                                         <  3:145364/148‑1 (MQ=255)
                                                                                                              cGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTttct                                      >  4:39382/1‑149 (MQ=255)
                                                                                                                    tAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTttctt                                     >  8:5107/1‑144 (MQ=255)
                                                                                                                    tAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTttctt                                     <  7:5107/144‑1 (MQ=255)
                                                                                                                    tAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTtcttct                                   >  8:2581/1‑146 (MQ=255)
                                                                                                                      aTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCa                                 <  7:2581/146‑1 (MQ=255)
                                                                                                                      aTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGa                               >  7:254504/1‑148 (MQ=255)
                                                                                                                      aTCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGa                               <  4:219817/148‑1 (MQ=255)
                                                                                                                         aCTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGGAAGCCACCTTCTGTAACTGGAAGGTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGatat                            >  5:399855/1‑148 (MQ=255)
                                                                                                                           tttGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATaaa                         >  5:52135/1‑149 (MQ=255)
                                                                                                                             tGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAAc                        <  7:72947/148‑1 (MQ=255)
                                                                                                                             tGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAAc                        <  8:198981/148‑1 (MQ=255)
                                                                                                                               gTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACggg                     >  2:99192/1‑149 (MQ=255)
                                                                                                                                  atatAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGGCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTAt                  <  3:39382/149‑1 (MQ=255)
                                                                                                                                    ataAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGATAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATTTTAACGGGTtttt                >  3:304646/1‑149 (MQ=255)
                                                                                                                                     taAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTg               <  1:172446/149‑1 (MQ=255)
                                                                                                                                     acaaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCACAGATATAAACGGGTAttt                >  3:266610/3‑148 (MQ=255)
                                                                                                                                      aaaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGaa                                                                   <  8:281291/96‑1 (MQ=255)
                                                                                                                                      aaaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGaa                                                                   >  7:281291/1‑96 (MQ=255)
                                                                                                                                      aaaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTttcttc                                    >  2:181757/1‑127 (MQ=255)
                                                                                                                                      aaaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTttcttc                                    <  1:181757/127‑1 (MQ=255)
                                                                                                                                        aaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATg            >  1:393107/1‑149 (MQ=255)
                                                                                                                                             ccGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGAt             >  8:52360/1‑143 (MQ=255)
                                                                                                                                             ccGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGAt             <  7:52360/143‑1 (MQ=255)
                                                                                                                                              cGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAAtt      <  6:460322/149‑1 (MQ=255)
                                                                                                                                              cGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTACGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAAtt      >  8:418794/1‑149 (MQ=255)
                                                                                                                                               gCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCt                                          >  1:202817/1‑112 (MQ=255)
                                                                                                                                               gCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCt                                          <  2:202817/112‑1 (MQ=255)
                                                                                                                                                ccTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTc    <  5:113377/149‑1 (MQ=255)
                                                                                                                                                    tcGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTTATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCt  <  7:148746/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCT  >  NC_000913/850847‑851141

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