Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,019,173 A→G K107E (AAG→GAG)  ubiE → bifunctional 2‑octaprenyl‑6‑methoxy‑1,4‑benzoquinone methylase/ S‑adenosylmethionine:2‑DMK methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,019,1730AG98.3% 221.4 / ‑3.6 58K107E (AAG→GAG) ubiEbifunctional 2‑octaprenyl‑6‑methoxy‑1,4‑benzoquinone methylase/ S‑adenosylmethionine:2‑DMK methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (26/31);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

GCAGCGGCGTACGCCGTGGGCAGACCGTGCTGGATCTGGCTGGTGGCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCA  >  NC_000913/4019027‑4019319
                                                                                                                                                  |                                                                                                                                                  
gCAGCGGCGTACGCCGTGGGCAGACCGTGCTGGATCTGGCTGGTGGCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCt                                                                                                                                                <  1:127106/151‑1 (MQ=255)
     ggcgTACGCCGTGGGCAGACCGTGCTGGATCTGGCTGGTGGCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGt                                                                                                                                            >  1:29527/1‑150 (MQ=255)
      gcgTACGCCGTGGGCAGACCGTGCTGGATCTGGCTGGTGGCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTaa                                                                                                                                          <  2:333188/151‑1 (MQ=255)
          aCGCCGTGGGCAGACCGTGCTGGATCTGGCTGGTGGCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATc                                                                                                                                      >  2:757102/1‑151 (MQ=255)
               gTGGGCAGACCGTGCTGGATCTGGCTGGTGGCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGgtg                                                                                                                                  >  1:301682/1‑150 (MQ=255)
                   gCAGACCGTGCTGGATCTGGCTGGTGGCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGAtt                                                                                                                             <  1:718956/151‑1 (MQ=255)
                     aGACCGTGCTGGATCTGGCTGGTGGCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGTTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTgg                                                                                                                           >  2:99497/1‑151 (MQ=255)
                       aCCGTGCTGGATCTGGCTGGTGGCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCa                                                                                                                         <  2:717865/151‑1 (MQ=255)
                        ccGTGCTGGATCTGGCTGGTGGCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCaa                                                                                                                        <  2:890187/151‑1 (MQ=255)
                        ccGTGCTGGATCTGGCTGGTGGCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCa                                                                                                                         >  1:122208/1‑150 (MQ=255)
                             cTGGATCTGGCTGGTGGCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTg                                                                                                                   >  1:307734/1‑151 (MQ=255)
                             cTGGATCTGGCTGGTGGCACCGGCGACCGGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGAAACGtt                                                                                                                    <  2:345254/150‑1 (MQ=255)
                                     ggctggTGGCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGt                                                                                                            <  2:619033/150‑1 (MQ=255)
                                      gctggTGGCACCGGCGCCCTGACAGCAAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTc                                                                                                          <  1:635838/151‑1 (MQ=255)
                                         ggtggCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAgg                                                                                                       <  2:1023120/151‑1 (MQ=255)
                                         ggtggCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAg                                                                                                        <  2:323064/150‑1 (MQ=255)
                                            ggCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGa                                                                                                    <  1:533577/151‑1 (MQ=255)
                                              cACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAAc                                                                                                  >  2:774675/1‑151 (MQ=255)
                                                ccGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGc                                                                                                <  2:297349/151‑1 (MQ=255)
                                                 cGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCt                                                                                               >  2:889141/1‑151 (MQ=255)
                                                  ggCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTg                                                                                              >  1:544241/1‑151 (MQ=255)
                                                      aCCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGgc                                                                                          >  2:37690/1‑151 (MQ=255)
                                                          gACAGCGAAATTCTCCCGCCTGGTCGGAGAAACAGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCt                                                                                       <  2:695492/150‑1 (MQ=255)
                                                                    ttCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCgg                                                                            >  2:939564/1‑151 (MQ=255)
                                                                              tGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCtt                                                                  >  2:515372/1‑151 (MQ=255)
                                                                                 tCGGCGAAACTGGCAAAGGGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGAGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGGAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGAGAACACCTTTg                                                                >  1:1063998/1‑150 (MQ=255)
                                                                                    gagaAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTg                                                            <  2:1064246/151‑1 (MQ=255)
                                                                                    gagaAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTg                                                            <  2:243316/151‑1 (MQ=255)
                                                                                      gaAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGc                                                           >  2:66224/1‑150 (MQ=255)
                                                                                        aaCTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATc                                                        >  2:1044716/1‑151 (MQ=255)
                                                                                         aCTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCa                                                       <  1:269377/151‑1 (MQ=255)
                                                                                           tGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCAcc                                                     >  2:554163/1‑151 (MQ=255)
                                                                                               aaaGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCAttt                                                 <  1:466869/151‑1 (MQ=255)
                                                                                                aaGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTc                                                >  1:374027/1‑151 (MQ=255)
                                                                                                aaGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTc                                                <  2:768912/151‑1 (MQ=255)
                                                                                                           gCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTg                                      <  2:772658/150‑1 (MQ=255)
                                                                                                           gCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTg                                      <  2:909817/150‑1 (MQ=255)
                                                                                                           gCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGc                                     >  1:971465/1‑151 (MQ=255)
                                                                                                            cTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCg                                    >  2:1020005/1‑151 (MQ=255)
                                                                                                             tGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGt                                   >  2:897256/1‑151 (MQ=255)
                                                                                                                tatCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTaa                                 >  2:612426/1‑150 (MQ=255)
                                                                                                                 atCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACg                               <  2:427094/151‑1 (MQ=255)
                                                                                                                 atCAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACg                               >  1:62134/1‑151 (MQ=255)
                                                                                                                   cAATGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTc                             <  1:187836/151‑1 (MQ=255)
                                                                                                                     aTGAATCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCa                            <  1:318090/150‑1 (MQ=255)
                                                                                                                        aaTCCATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGa                        <  1:227647/151‑1 (MQ=255)
                                                                                                                           ccATGCTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACaa                     <  1:436884/151‑1 (MQ=255)
                                                                                                                                cTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGAt                 >  2:883114/1‑150 (MQ=255)
                                                                                                                                cTCAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGAt                 >  1:981831/1‑150 (MQ=255)
                                                                                                                                  cAAAATGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATa                <  1:674353/149‑1 (MQ=255)
                                                                                                                                   aaaaTGGGCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAg             <  2:122209/151‑1 (MQ=255)
                                                                                                                                        gggCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTg        <  1:443167/151‑1 (MQ=255)
                                                                                                                                         ggCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTg        >  2:167572/1‑150 (MQ=255)
                                                                                                                                          gCCGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCg      >  2:327113/1‑151 (MQ=255)
                                                                                                                                           ccGCGAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGt     <  2:965722/151‑1 (MQ=255)
                                                                                                                                            cgcgAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGtt    >  1:865342/1‑151 (MQ=255)
                                                                                                                                              cgAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCa  <  2:28197/151‑1 (MQ=255)
                                                                                                                                              cgAGGAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCa  <  1:894825/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCAGCGGCGTACGCCGTGGGCAGACCGTGCTGGATCTGGCTGGTGGCACCGGCGACCTGACAGCGAAATTCTCCCGCCTGGTCGGAGAAACTGGCAAAGTGGTCCTTGCTGATATCAATGAATCCATGCTCAAAATGGGCCGCGAGAAGCTGCGTAATATCGGTGTGATTGGCAACGTTGAGTATGTTCAGGCGAACGCTGAGGCGCTGCCGTTCCCGGATAACACCTTTGATTGCATCACCATTTCGTTTGGTCTGCGTAACGTCACCGACAAAGATAAAGCACTGCGTTCA  >  NC_000913/4019027‑4019319

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: