Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 850,995 T→G 7.8% 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 7.8% 216.2 / 10.8 64intergenic (+30/+19)ompX/opgEouter membrane protein X/OPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base T (26/33);  new base G (2/3);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

GCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAA  >  NC_000913/850849‑851135
                                                                                                                                                  |                                                                                                                                            
gCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCg                                                                                                                                            <  4:258903/149‑1 (MQ=255)
     tCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCg                                                                                                                                       >  1:449609/1‑149 (MQ=255)
       tGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGa                                                                                                                                     <  1:294974/149‑1 (MQ=255)
        gCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAt                                                                                                                                    <  4:443933/149‑1 (MQ=255)
           gTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttt                                                                                                                                 <  3:338788/149‑1 (MQ=255)
            ttCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttttt                                                                                                                                <  4:356177/149‑1 (MQ=255)
                aCCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGttt                                                                                                                            >  1:109662/1‑149 (MQ=255)
                  ccGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGttttt                                                                                                                          >  4:128033/1‑149 (MQ=255)
                                  ctctGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCg                                                                                                                 >  1:381308/1‑142 (MQ=255)
                                  ctctGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCg                                                                                                                 <  2:381308/142‑1 (MQ=255)
                                     tGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGc                                                                                                                                        <  2:439908/116‑1 (MQ=255)
                                     tGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGc                                                                                                                                        >  1:439908/1‑116 (MQ=255)
                                         cTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCg                                                                                                   >  1:381300/1‑149 (MQ=255)
                                         cTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCg                                                                                                   <  2:109662/149‑1 (MQ=255)
                                          ttCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCgg                                                                                                  >  2:594469/1‑149 (MQ=255)
                                            ctctTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGtt                                                                                                >  1:586856/1‑149 (MQ=255)
                                             tctTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTc                                                                                               >  1:573839/1‑149 (MQ=255)
                                               ttACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCg                                                                                             >  2:546414/1‑149 (MQ=255)
                                               ttACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAt                                                                                                                                    >  4:469397/1‑110 (MQ=255)
                                               ttACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAt                                                                                                                                    >  2:311846/1‑110 (MQ=255)
                                                    aGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATg                                                                                        <  1:60266/149‑1 (MQ=255)
                                                     gCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATgg                                                                                       <  1:546414/149‑1 (MQ=255)
                                                     gCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTCCGGGTCGAAAATATCGGTTCCGAGATgg                                                                                       <  1:660640/149‑1 (MQ=255)
                                                       agagCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGtt                                                                                     <  3:128033/149‑1 (MQ=255)
                                                       agagCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGtt                                                                                     <  3:599199/149‑1 (MQ=255)
                                                        gagCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGGCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTcc                                                                                              >  2:556320/1‑139 (MQ=255)
                                                        gagCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTcc                                                                                              <  1:556320/139‑1 (MQ=255)
                                                           ccGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTaa                                                                                 >  1:508257/1‑149 (MQ=255)
                                                              tATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATcc                                                                              <  3:557635/149‑1 (MQ=255)
                                                               aTTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAttttt                                                                                                                                >  3:193616/1‑98 (MQ=255)
                                                               aTTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAttttt                                                                                                                                <  4:193616/98‑1 (MQ=255)
                                                                    tAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTaa                                                                                                                        <  4:436453/101‑1 (MQ=255)
                                                                    tAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTaa                                                                                                                        >  3:436453/1‑101 (MQ=255)
                                                                     aGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAAACCACCTTCt                                                                       <  1:594469/149‑1 (MQ=255)
                                                                         ttGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTaa                                                                   >  2:123171/1‑149 (MQ=255)
                                                                         ttGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTaa                                                                   <  2:573839/149‑1 (MQ=255)
                                                                               tAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGaa                                                             <  2:278779/149‑1 (MQ=255)
                                                                                 ggCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGt                                                           <  4:673570/149‑1 (MQ=255)
                                                                                       tGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGa                                                     <  2:449609/149‑1 (MQ=255)
                                                                                        ggATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGAt                                                    <  2:381300/149‑1 (MQ=255)
                                                                                          aTTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATga                                                  <  2:586856/149‑2 (MQ=255)
                                                                                              ccGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACt                                                                 >  2:606050/1‑130 (MQ=255)
                                                                                              ccGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACt                                                                 <  1:606050/130‑1 (MQ=255)
                                                                                               cGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATc                                             >  1:129144/1‑149 (MQ=255)
                                                                                                       gTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTaa                                                                                                                        >  3:409162/1‑66 (MQ=255)
                                                                                                       gTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTaa                                                                                                                        <  4:409162/66‑1 (MQ=255)
                                                                                                            cGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATg                                                                                        >  4:43849/1‑93 (MQ=255)
                                                                                                            cGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATg                                                                                        <  3:43849/93‑1 (MQ=255)
                                                                                                            cGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTttct                                <  1:123171/149‑1 (MQ=255)
                                                                                                                 cTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCa                           <  2:92124/149‑1 (MQ=255)
                                                                                                                     tCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCAcc                                                                           <  2:545211/97‑1 (MQ=255)
                                                                                                                     tCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCAcc                                                                           >  1:545211/1‑97 (MQ=255)
                                                                                                                     tCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACt                                                                 >  2:344644/1‑107 (MQ=255)
                                                                                                                     tCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACt                                                                 <  1:344644/107‑1 (MQ=255)
                                                                                                                            ggTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACgg                <  2:508257/149‑1 (MQ=255)
                                                                                                                               gATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTa             >  4:113432/1‑149 (MQ=255)
                                                                                                                                     aaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGt                                                       >  3:8839/1‑101 (MQ=255)
                                                                                                                                     aaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGt                                                       <  4:8839/101‑1 (MQ=255)
                                                                                                                                     aaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGAt       <  2:129144/149‑1 (MQ=255)
                                                                                                                                          tCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCa                           <  3:196586/124‑1 (MQ=255)
                                                                                                                                          tCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCa                           >  4:196586/1‑124 (MQ=255)
                                                                                                                                          tCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTaa  >  2:65894/1‑149 (MQ=255)
                                                                                                                                             gCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGt     <  4:105656/143‑1 (MQ=255)
                                                                                                                                             gCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGt     >  3:105656/1‑143 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
GCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAA  >  NC_000913/850849‑851135

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑